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ABSTRACT 


We  report  the  first  observations  of  photon-assisted  quasiparticle 
tunneling  and  AC  Josephson  effect  in  superconducting  tunnel  junctions 
Irradiated  with  far-infrared  (FIR)  radiation.  Radiation  at  246  GHz 
(A*1.22  nin)  and  604  GHz  (496  ym)  from  an  optically- pumped  FIR  laser 
source  was  used.  Tin-tin  oxide-lead  junctions  of  ^1  (pm)  area  were 
fabricated  on  crystal  quartz  substrates  with  integral  planar  dipole  an¬ 
tennas  of  resonant  length  at  the  frequency  of  the  incident  radiation. 
The  observed  photon-assisted  tunneling  features  are  in  excellent  agree¬ 
ment  with  the  Tien-Gordon  theory,  and  the  inferred  responsivity  ap¬ 
proaches  the  quantum  limit  at  low  temperatures  for  photon  energies  less 
than  the  gap.  At  604  GHz,  with  a  176  ohm  junction,  we  have  seen  7 
Josephson  steps,  comparable  to  point  contact  performance.  The  varia¬ 
tion  of  the  step  widths  with  laser  power  is  found  to  agree  quite  well 
with  both  the  RSJ  model  and  the  Werthamer  theory.  For  low  resistance 
junctions  (e.g.  16  ohms),  we  find  the  Josephson  steps  to  be  flat  and  to 
agree  well  with  the  shape  predicted  by  the  RSJ  model  without  noise 
rounding,  while  noise  rounding  is  very  evident  with  the  higher  resis¬ 
tance  junctions. 

The  step  shape  is  in  reasonable  agreement  with  the  theory  of  P.  A. 
Lee,  using  a  noise  temperature  of  20-25  K,  attributed  to  shot  noise  at 
finite  voltage,  with  an  appropriate  small  junction  capacitance.  A 
computer  simulation  has  been  developed  which  models  the  effects  of  the 
voltage-dependent  (shot)  noise  and  the  non-linear  quasiparticle  resis¬ 
tance  of  the  junctions,  assuming  T*2K.  Reasonably  good  agreement  is 


obtained  between  the  simulated  I-V  curves  for  un-lrradiated  but  noise 
affected  junctions,  reproducing  the  observed  fall -off  In  the  IcRn 
product  with  Increasing  junction  resistance.  For  Irradiated  junctions 
we  also  obtain  good  agreement  with  the  observed  noise  rounding  of  the 
ac  Josephson  steps.  Finally  we  have  coupled  FIR  radiation  Into  a  7n 
niobium-sill con-niobium  junction  and  have  observed  three  Josephson 
steps  and  one  photon-assisted  tunneling  step. 
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FIG.  1.1  (a)  Schematic  illustration  of  the  dc  electrical 

measurement  of  an  SIS  tunnel  junction,  (b)  DC  current- 
voltage  (I-V)  curve  of  a  16  n  Sn-SnO-Pb  tunnel  junction 
taken  at  T*1.4  K. 

FIG.  1.2  (a)  Effect  of  32  GHz  microwave  radiation  on 
superconducting  weak  link,  illustrating  the  flat 
current  steps  from  the  ac  Josephson  effect,  (b)  Effect 
of  38  GHz  microwave  radiation  on  an  Al-oxide-Pb  junc¬ 
tion,  illustrating  the  photon-assisted  tunneling 
current  steps. 

FIG.  2.1  Real  and  imaginary  parts  of  the  response  func¬ 
tions  j-|(u)  and  jgtaJ.  at  T=0. 

FIG.  2.2  (a)  Critical  current  Ir(T),  normalized  to  Ic(0), 

vs.  the  reduced  temperature  T/T  for  a  Sn-Sn  tunnel 
junction,  (b)  Critical  current,  IC(T),  normalized  to 
1(0),  vs.  the  reduced  tempera turecT/T  for  a  Sn-Pb 
tonne 1  junction.  c 

FIG.  2.3  Schematic  illustration  of  the  current  biased  RSJ 
model. 

FIG.  2.4  Mechanical  analogy  of  the  RSJ  model  with  a  par¬ 
ticle  moving  In  a  tilted  periodic  potential. 

FIG.  2.5  Effect  of  noise  on  the  RSJ  model  I-V  curve  for 

Bc-0. 

FIG.  3.1  (a)  Section  of  the  commercial  mask  containing  a 

tunnel  junction  pattern  and  a  pad  pattern,  (b)  A 
close-up  of  the  central  region  of  (a)  showing  the 
tunnel  junction  pattern,  (c)  Another  section  of  the 
mask,  this  time  with  a  tunnel  junction  pattern  with 
long  leads,  (d)  Section  of  mask  containing  a  long 
microbridge  pattern. 

FIG.  3.2  Illustration  of  the  sequence  of  steps  used  to 
make  a  mask. 

FIG.  3.3  Photoresist  exposure/development  procedure 

beginning  with  the  photoresist/  A1 /photoresist  tri layer 
and  ending  with  the  suspended  photoresist  bridge  which 
Is  used  to  make  the  junction. 


FIG.  3.4  (a)  SEM  micrograph  of  a  photoresist  bridge  sus¬ 
pended  over  a  polished  sapphire  substrate,  made  with  the 
mask  In  Figure  3.1(c),  following  the  procedure  of  Fig. 
3.3.  (b)  SEM  micrograph  of  a  suspended  photoresist 

bridge  used  in  fabricating  a  small-area  junction  at  the 
center  of  a  dipole  antenna. 

FIG.  3.5  (a)  Evaporation/oxidation  steps  for  Sn-SnO-Pb 
overlap  junctions,  (b)  SEM  micrograph  of  a  junction 
made  with  this  procedure. 

FIG.  3.6  (a)  Evaporation/oxidation  procedure  for  resist- 

aligned  junctions,  (b)  SEM  micrograph  of  a  sample  made 
with  this  procedure. 

FIG.  3.7  (a)  Evaporation/oxidation  procedure  for  edge- 

aligned  junctions,  (b)  SEM  micrograph  of  junction  made 
with  this  procedure. 

FIG.  3.8  (a)  Reduction  factor,  a,  as  a  function  of  the 
effective  thickness,  te^,  for  different  values  of  er. 
(b)  Values  of  a  for  teff>0.15  for  the  values  of  er  in 

(a) ,  (c)  E-plane  and  H-plane  antenna  patterns  for  a 

dipole  antenna  on  a  thin  substrate. 

FIG.  3.9  (a)  Current  element  at  an  air-dielectric  inter¬ 
face.  (b)  E-plane  and  H-plane  power  patterns  from  this 
current  element. 

FIG.  3.10  (a)  H-plane  power  pattern  of  a  resonant  dipole 
antenna  at  an  air-dielectric  interface  for  different 
values  of  er.  (b)  E-plane  power  pattern  of  the  reso¬ 
nant  dipole  of  (a)  for  the  same  values  of  Er. 

FIG.  3.11  Detail  of  transverse  optical  access  dewar.  FIR 
laser  radiation  is  focussed  by  the  TPX  lens  onto  the 
cone,  through  the  back  side  of  the  substrate,  onto  the 
junction-antenna  structure. 

FIG.  3.12  (a)  SEM  micrograph  of  a  typical  sample  with  a 
dipole  antenna  of  ~183  urn  length  and~5  urn  line  width. 

At  the  center  of  the  dipole  is  the  tunnel  junction. 

(b)  Close-up  of  the  junction  region  in  (a). 

FIG.  3.13  Methyl  fluoride  energy-level  diagram.  The  C0« 
pump  line  coincides  with  the  v=0-l ,  J=ll,  K=2  transi¬ 
tion.  The  496  urn  FIR  lasing  transition  is  the  v=l,  K=2, 
J*12-ll  transition. 

FIG.  3.14  Optically-pumped  FIR  laser  set-up. 


FIG.  3.15  A  24  minute  recording  of  the  free-running  FIR 
laser  output  power. 

FIG.  4.1  Typical  I-V  curves  of  a  low  resistance  (~  2  n), 
larger  area  (-  4.5xl0“8  cm2),  tunnel  junction  for  de¬ 
creasing  temperatures. 

FIG.  4.2  The  square  of  2Acn,  measured  from  data  such  as 
that  presented  in  Fig.  4.1,  versus  the  temperature  T. 

FIG.  4.3  (a)  The  critical  current  at  a  temperature  T, 

normalized  to  the  critical  current  at  1.4  K,  Ic(T)/Ic(0), 
versus  the  reduced  temperature,  t=T/Tg.  (b)  The  criti¬ 
cal  current  at  a  temperature  T,  normalized  to  the 
estimated  zero  temperature  critical  current, 

IC(0)=1 .2±0.1  mV,  Uc(T)/Ic{0)}est  versus  t=T/Tc. 

FIG.  4.4  Critical  current,  Ig,  vs.  the  normal  resistance, 
Rp,  for  the  junctions  studied  in  this  work.  Data  shown 
by  solid  circles  falls  below  the  dashed  line,  IcRn=1.2mV 
for  Rn>10  fl.  The  computer  simulation  solid  line  and 
open  circles,  closely  follows  the  data.  Error  bars 
indicate  bounds  of  simulated  Ic  for  a  given  In. 

FIG.  4.5  (a)  Typical  I-V  curves  of  a  low- resistance  sample 

for  increasing  laser  power.  Here,  Vgs(Apb+Acn)/e. 

(b)  Typical  I-V  curves  of  a  high-resistence  junction  for 
increasing  laser  powers. 

FIG.  4.6  Effect  of  the  polarization  of  the  FIR  laser  out¬ 
put  relative  to  the  antenna  axis.  Parallel  polarization, 
//,  couples  well,  whereas  perpendicular  polarization,  J., 
does  not,  as  expected. 

FIG.  4.7  Normalized  step  width  as  a  function  of  the  laser 
voltacje  in  the  junction  compared  with  the  Bessel  func¬ 
tion  J Jn  J  dependence  and  the  Werthamer  theory  |Wn|,  for 
a  composite  of  data  from  two  low  resistance  junctions. 

FIG.  4.8  Normalized  step  width  as  a  function  of  the  laser 
voltage  in  the  junction,  compared  with  the  Bessel  func¬ 
tion  j  J  |  dependence  and  the  Werthamer  theory,  |  W  |. 

The  data  is  from  a  high-resistance  (176  n)  junction 
irradiated  by  604  GHz  laser  radiation. 

FIG.  4.9  Normalized  step  width  as  a  function  of  laser 
voltage  in  the  junction,  compared  with  the  Bessel  func¬ 
tion  dependence,  |0|,  and  the  Werthamer  theory,  |W|. 

The  data  shown  is  from  a  156  n  junction  irradiated  by 
246  GHz  laser  radiation. 


lx 


FIG.  4.10  I-V  curves  of  two  different  junctions,  each  with 
suitable  antenna  length,  irradiated  at  604  and  246  GHz 
at  1.4  K.  The  dashed  curves  are  obtained  from  the  Tien- 
Gordon  theory.  This  accounts  for  the  photon-assisted 
tunneling  steps  at  ±Vg±kfiu)|_/e  where  k  is  an  integer  and 
vg*^Pb+ASn^e  *  leav*n9  only  local  voltage  displacements 
nfiar  tne  Josephson  steps  labeled  by  n=0,l,2,...  . 

FIG.  <■,  (a)  Superposed  dc  I-V  curves  for  different  values 

of  trie  normalized  FIR  laser-induced  voltage  o=eV|_/fiu)|_  for 
a  176  n  junction  at  4.2  K  and  a  laser  frequency  of  604 
GHz.  (b)  Comparison  between  the  voltage  dependence  of 
dV/dl  computed  theoretically  for  various  FIR  power  levels 
and  that  obtained  by  digital  differentiation  of  the 
experimental  I-V  curves  in  (a). 

FIG.  4.12  Change  in  the  dc  current  with  laser  power  [~c?) 
at  selected  dc  bias  voltages.  Points  are  experimental 
data;  curves  are  computed  theoretically  from  the  zero- 
power  I-V  curve. 

FIG.  4.13  Plot  of  current  responsivity  normalized  to  the 
quantum  limit  (e/lta)  vs.  dc  bias  voltage.  Inset  shows 
computed  frequency  dependence  of  R/(e/t>w)  at  T=0,  1.4 
and  4.2  K,  for  a  bias  of  V^V-,  and  experimental  points 
at  246  and  604  GHz  at  4.2  K. 

FIG.  4.14  Fit  to  the  shape  of  the  2nd  Josephson  step  on  a 
16  n  junction  illuminated  with  604  GHz  laser  radiation. 

FIG.  4.15  Fits  to  the  shape  of  the  4th  Josephson  step  on 
a  176  n  sample  at  604  GHz  for  different  values  of  the 
fitting  parameters  l£,  Teff,  and  C. 

FIG.  4.16  (a)  Simulated  I-V  curve  with  noise  but  no  ac 
currents.  Switching  at  critical  current  shown  by  dashed 
lines.  Critical  currents  from  simulation  (for  5  runs) 
fall  within  bounds  indicated  by  dotted  lines,  (b)  I-V 
curve  from  the  176  n  junction  (data). 

FIG.  4.17  Zeroth  and  first  step  simulated  behavior  for  a 
junction  with  rf  and  dc  bias  but  no  noise. 

FIG.  4.18  Simulated  step  at  2a=4.0  Arrow  labeled  by  1°  is 
the  unfluctuated  critical  current  estimated  from  the 
nominal  I-Rp  product.  The  critical  current  with  noise 
is  also  shown,  as  well  as  the  estimated  errors  on  the 
simulated  critical  current,  i.e.,  the  error  bars  drawn. 


FIG.  4.19  (a)  I-V  curve  of  a  <76  n  junction  at  a*2.0; 
Josephson  steps  are  labelled  by  the  integers  n*2,  ns3, 
and  ns4.  (b)  Simulation  of  the  data  shown  in  (a),  which 
includes  noise  and  photon-assisted  tunneling  steps. 

(c)  Simulation,  dashed  lines  in  (b),  compared  with  data 
from  (a)  illustrating  the  excellent  quality  of  the  fit. 

FIG.  4.20  Fourth  Josephson  step  fits  at  two  laser  power 
levels,  2o«4.0  and  2a*5.75.  On  the  left  are  the  simu¬ 
lated  step  shapes,  measured  from  the  step  center,  solid 
circles,  compared  with  the  data,  solid  line.  On  the 
right  is  the  same  data,  but  the  simulation  now  includes 
the  photon-assisted  tunneling  effect. 

FIG.  5.1  Input  impedance  of  a  bow-tie  antenna  as  a  func¬ 
tion  of  bow  angle  for  an  antenna  on  a  fused  quartz  sub¬ 
strate.  Solid  lines  are  calculated  values,  crosses  are 
values  measured  from  a  10  GHz  microwave  modelling 
experiment. 

FIG.  A.l  Illustration  of  cross-section  of  plexiglass  con¬ 
tact  printing  sample  holder. 

FIG.  B.l  A  typical  high  current-density  sample  showing 
depressed  gap  structure  from  quasiparticle  injection  at 
6  temperatures.  At  1.4  K  a  scaled  I-V  curve  from  a 
high-resistance  junction  is  drawn  in  dashed  lines  for 
comparison.  The  three  upper  curves  are  drawn  with  a 
scale  twice  as  sensitive  as  the  lower  three  curves. 

FIG.  B.2  An  I-V  curve  for  a  Sn-SnO-Sn  junction  (Ref.  116) 
showing  the  same  sort  of  gap  depression  phenomena  ob¬ 
served  in  this  work  on  Sn-SnO-Pb  junctions. 

FIG.  B.3  Voltage  difference,  Vg,  between  the  dashed  curve 
and  solid  curve  in  Fig.  B.l  at  T*1.4  K  versus  the  cur¬ 
rent  in  the  junction. 

FIG.  B.4  Voltage  difference,  V_,  between  gap  structure  of 
low  resistance  and  high  resistance  junctions,  computed 
in  the  same  way  as  Fig.  B.3,  for  the  Sn-SnO-Sn  junc¬ 
tions  of  Ref.  116  as  a  function  of  JG,  which  is  propor¬ 
tional  to  the  current  through  the  junction. 

FIG.  B.5  I-V  curves  of  a  0.55  n  edge-aligned  point  contact 
type  sample  at  2a*0,  and  2aa0.8,  the  maximum  power 
coupled  to  the  junction,  where  2as2eVL/fiu>|_.  The  curve  at 
2a*0.8  clearly  shows  a  step  at  1.25  mV,  labelled  n=l . 
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FIG.  B.6  I-V  curves  at  several  laser  powers  for  a  14  n  187 

point  contact  edge-aligned  sample  showing  gap  structure 
at  (aPb+Acn)/e  at  2a*0,  and  up  to  five  Josephson  steps 
at  2a*6.7.  Note  the  presence  of  the  half-harmonic  step, 
labelled  n*»s,  not  previously  seen  In  point  contact  data 
in  the  FIR. 

FIG.  B.7  Step-widths  as  a  function  of  normalized  laser  188 

Induced  voltage  in  the  junction,  2a,  for  the  steps  shown 
in  Fig.  B.6.  Solid  lines  are  scaled  Werthamer  functions, 

I W  ,  and  the  dashed  lines  are  scaled  Bessel  functions, 

|J  . 

FIG.  B-8  Maximum  value  of  the  normalized  laser  induced  190 

voltage  in  the  junction,  2as2eV|_/ftu>(_ ,  versus  the  normal 
resistance,  Rj,. 

FIG.  B.9  SEM  micrographs  of  a  SNAP  junction,  (a)  Low  mag-  193 

nification  view  of  sample,  displaying  dc  leads— vertical 
parallel  lines  in  the  Figure—and  pads  to  which  gold 
wires  are  bonded,  (b)  Close-up  of  antenna  region  of  (a), 
showing  10  ym  linewidth  antenna,  ~180  ym  long.  The 
junction  is  at  the  center  of  the  antenna,  dc  leads  are 
to  the  left  of  the  center,  (c)  Close-up  of  junction 
region  at  the  center  of  (b)  taken  with  SEM  stage  tilted. 

The  junction  Is  the  smallest  area  at  the  center  of  the 
two  larger  squares,  (d)  Close-up  of  junction  region  but 
without  a  tilt  on  the  SEM  stage,  the  junction  is  the 
“1  ym  by  -1  ym  region  at  the  center  of  the  Figure. 

FIG.  B.10  I-V  curves  of  the  SNAP  junction  in  Fig.  B.9,  195 

Irradiated  with  604  GHz  (A=496  ym)  laser  radiation,  at 
normalized  laser  Induced  voltages,  2a=0,  and  the  highest 
voltage  coupled,  2a*1.7.  We  see  three  Josephson  steps 
at  Vn*n(fiw|  /2e)  where  n  is  an  integer,  and  a  single 
photon-assisted  tunneling  step  at  Vg+(fiu)L/e)  where 
Vg«2ANb,  as  labelled  in  the  Figure. 

FIG.  B.ll  Step  half-widths  as  a  function  of  normalized  197 

laser  Induced  voltage,  2a,  for  the  Josephson  steps  of 
Fig.  B.10,  at  two  temperatures  T=4.2  K  (crosses),  and 
T*1.4  K  (solid  circles).  Solid  line  denotes  Werthamer 
function  fits,  |W|,  and  dashed  lines  denote  Bessel  func¬ 
tion  fits,  | J | . 
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CHAPTER  1 


INTRODUCTION 


In  the  26  years  since  the  development  of  the  microscopic  theory  of 
superconductivity  by  Bardeen,  Cooper,  and  Schrieffer*,  research  in 
superconductivity  has  followed  two  divergent  paths.  One  path  has 
brought  the  development  of  high  transition  temperature  (Tc),  high 
critical  magnetic  field  (Hc2)  superconducting  materials  for  use  In  large 
scale  applications  such  as  beam  magnets  for  high  energy  proton  and 
electron  accelerators,  or  magnetic  confinement  "mirrors"  for  fusion 
applications  .  These  developments  depended  critically  on  the  discovery 

3 

of  A15  compounds  by  Hulm,  Matthias,  and  Kunzler  .  The  second  path  of 
research  In  superconductivity  relied  on  two  discoveries  related  to  the 
problem  of  electrons  tunneling  between  two  superconductors  separated  by 
a  thin  layer  of  Insulator.  The  first  discovery  was  of  quasiparticle  or 
single  electron  tunneling  between  superconductors  by  Giaever  ,  the 
second  was  the  discovery  of  the  tunneling  of  Cooper  pairs  of  electrons 
by  Josephson^.  The  work  reported  here  contributes  to  this  second  path 
by  making  a  systematic  study  of  electron  tunneling  between 
superconductors  Irradiated  in  the  far-infrared  (FIR)  region  of  the 
spectrum. 

We  will  start  by  giving  a  qualitative  description  of  the  dc 
electrical  behavior  of  a  superconductor- Insulator-superconductor  (SIS) 
tunnel  junction.  The  current-voltage  (I-V)  characteristics  of  the  SIS 
tunnel  junction  have  two  distinct  branches.  Fig.  1.1(a)  shows 
schematically  such  a  tunnel  junction  while  Fig.  1.1(b)  shows  Its  I-V 


curve.  In  the  current-biased  dc  behavior,  as  we  Increase  the  bias 
current,  we  find  a  portion  of  the  I-V  curve  In  which  there  Is  zero 
voltage  for  finite  currents.  This  is  the  case  until  the  current  exceeds 
a  characteristic  current  called  the  critical  current,  I  ,  at  which  point 
the  junction  switches  to  the  finite  voltage  state.  At  high  currents  the 
junction  I-V  curve  becomes  approximately  linear.  As  the  current  is 
reduced  towards  zero  the  junction  voltage  stays  approximately  constant 
at  a  characteristic  voltage  called  the  gap  voltage  or  Vgap«(A1  +  A2)/e. 
This  effect  comes  from  the  high  differential  conductance  resulting  from 
the  divergence  In  the  density  of  states  of  the  quasiparticle  excitations 
at  the  superconducting  energy  gap.  Finally  as  the  current  approaches 
zero,  the  junction  voltage  suddenly  becomes  zero.  Usually  a  large 
hysteresis  loop,  beginning  where  the  junction  Initially  goes  to  the 
finite-voltage  state  to  where  It  drops  back  into  the  zero-voltage  state. 
Is  found.  It  Is  the  ability  of  the  junction  to  switch  rapidly  between 
two  states  which  makes  SIS  devices  attractive  for  high-speed  computing 
applications6. 

Other  interesting  effects  occur  when  an  SIS  tunnel  junction  is 
illuminated  by  high  frequency  electromagnetic  radiation.  One  effect 
depends  on  the  Cooper  pairs  and  is  called  the  ac  Josephson  effect  in 
which  constant  voltage  "steps"  appear  on  the  I-V  curve  at  multiples  of  a 
characteristic  voltage,  fiu)L/2e  ,  shown  in  Fig.  1.2(a).  These  steps  are 
due  to  the  mixing  of  the  ac  supercurrent  with  the  applied  high  frequency 
electromagnetic  field  in  the  junction.  A  second  effect  comes  from  the 
interaction  of  the  quasi  particles  with  the  external  field,  and  is  called 
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FIG.  1.1  (a)  Schematic  illustration  of  the  dc  electrical  measurement 

of  an  SIS  tunnel  junction,  (b)  DC  current-voltage  (I-V)  curve  of  a  16  n 
Sn-SnO-Pb  tunnel  junction  taken  at  T*1.4  K. 


the  photon-assisted  tunneling  effect.  Due  to  this  effect,  "steps" 


appear  on  the  I-V  curve  at  characteristic  voltages  ±vgap+nCfi“i/e)  where 
n  Is  a  positive  or  negative  Integer.  This  effect  Is  shown  schematically 
In  Fig.  1.2(b).  The  photon-assisted  tunneling  (PAT)  effect  was 
discovered  In  1962  by  Dayem  and  Martin7  and  explained  in  1963  by  Tien 

o 

and  Gordon  . 

The  complete  frequency-dependent  theory  of  the  SIS  tunnel  junction. 
Including  the  PAT  effect  and  the  ac  Josephson  effect  was  developed  by 

Q 

Werthamer  In  1966.  Because  the  tunnel  junction  is  so  sensitive  to 
external  electromagnetic  fields  due  to  the  above-mentioned  effects, 
tunnel  junctions  and  other  superconducting  devices  have  been  proposed 
for  use  as  sensitive  detectors  of  millimeter  and  submillimeter 
wavelength  radiation10. 

In  this  report  we  will  be  concerned  with  the  quantitative 

measurement  of  the  Josephson  effect  and  the  photon-assisted  tunneling 

effect  In  the  sub-millimeter  wavelength  region  of  the  spectrum, 

particularly  near  and  above  the  superconducting  energy  gap.  The 

measurements  described  herein  are  the  first  reported  measurements  in  the 

far-infrared  (FIR)  on  these  effects  on  fabricated  thin-film  tunnel 

junctions11.  Prior  to  this  work  the  only  data  available  on 

12 

superconducting  devices  In  the  FIR  were  on  point-contact  devices  which 
are  mechanically  unstable  and  known  to  have  a  large  variety  of  I-V 
curves13. 
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Using  small -area  superconducting  tunnel  junctions  we  have  observed 

the  ac  Josephson  effect  and  the  photon-assisted  tunneling  effect  at 

frequencies  up  to  604  GHz,  the  highest  frequency  reported  to  date.  We 

12 

have  seen  a  maximum  of  seven  ac  Josephson  steps  ,  up  to  8.75  mV,  a 

result  comparable  to  that  of  point  contacts.  We  have  made  a  systematic 

study  of  the  ac  Josephson  effect  and  the  PAT  effect  and  found  good 

agreement  between  our  data  and  the  predictions  based  on  fundamental 

theory.  On  high  resistance  junctions,  R>100  Q  ,  we  have  observed 

rounded  Josephson  steps.  We  have  shown  that  the  rounding  is  due  to  shot 
14 

noise  ,  in  contrast  to  the  results  on  point  contacts  which  were 

consistent  with  thermal  noise  from  joule  heating  of  the  point  contact  at 

12 

high  dc  bias  voltages  . 

This  report  is  organized  into  five  chapters  and  two  appendixes, 
beginning  with  this  brief  introductory  Chapter.  Chapter  II  begins  with 
a  detailed  theoretical  review,  describing  the  fundamental 
frequency-dependent  theory  of  the  ac  Josephson  effect,  the 
photon-assisted  tunneling  effect,  and  the  fundamental  theory  of  noise 
effects  In  these  devices.  After  this  theoretical  discussion,  we  briefly 
describe  the  Important  experimental  results  obtained  prior  to  this  work. 
In  Chapter  III  we  describe  the  experimental  techniques  which  were  used 
to  make  our  measurements,  including  the  fabrication  of  the  samples  which 
were  made  at  Harvard,  the  means  of  coupling  FIR  radiation  to  the 
small -area  tunnel  junctions,  and  the  FIR  laser  system  used  to  irradiate 
the  junctions.  Chapter  IV  is  a  detailed  analysis  of  these  measurements, 
the  Josephson  effect  results,  the  photon-assisted  tunneling  effect 
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FIG.  12  (a)  Effect  of  32  GHz  microwave  radiation  on  superconducting 

weak  link,  illustrating  the  flat  current  steps  from  the  ac  Josephson 
effect,  (b)  Effect  of  38  GHz  microwave  radiation  on  an  Al-oxide-Pb 
junction,  illustrating  the  photon-assisted  tunneling  current  steps. 
[After  Ref.  3.] 
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results,  end  noise  effects.  The  main  body  of  the  report  ends  with 
Chapter  V,  which  discusses  the  results  of  the  FIR  measurements  as  a 
whole,  proposes  a  set  of  experiments  to  do  in  the  future,  and  briefly 
discusses  the  implications  of  this  work  for  FIR  mixers  and  direct 
detectors. 

The  two  appendixes  are  organized  as  follows.  In  Appendix  A  we  give 
complete  details  of  the  fabrication  process  including  the 
photolithography,  substrate  cleaning  procedures,  and  chemical  etches  for 
the  various  metals  used.  Appendix  B  begins  by  describing  measurements 
on  high-current  density  junctions  made  at  Bell  Labs  which  showed 
unusually  large  amounts  of  gap  suppression  due  to  heavy  quasiparticle 
injection.  Following  this,  data  on  point-contact  type  junctions  and 
preliminary  results  from  Nb-Si-Nb  junctions  made  at  Sperry  Research 
Labs,  are  presented. 


CHAPTER  II 


THEORETICAL  OVERVIEW 


2.1  Introduction 

This  Chapter  is  concerned  with  the  quantitative,  fundamental  theory 
of  the  phenomena  described  previously  in  Chapter  It  that  is,  why  a  SIS 
tunnel  junction  has  a  non-linear  I-V  curve  with  two  distinct  branches, 
and  why  the  1-V  curve  develops  ac  Josephson  steps  and  photon-assisted 
tunneling  steps  when  the  junction  is  irradiated  by  high-frequency 
radiation. 

We  will  sketch  the  main  features  of  the  theory  of  the  tunnel 
junction,  relevant  to  this  experiment,  beginning  in  Section  2.2.  In 
Section  2.2  the  tunneling  Hamiltonian  formalism  will  be  described  as 
well  as  its  specific  predictions  for  the  dc  and  ac  voltage  biased  tunnel 
junction  cases.  The  RSJ  model,  a  commonly  used  low-frequency  version  of 
the  theory,  is  presented  in  Section  2.3,  while  the  effects  of  noise  on 
the  dc  and  high-frequency  behavior  of  the  junction  are  outlined  in 
Section  2.4.  The  Chapter  ends  with  a  brief  discussion  of  the 
experimental  results  obtained  by  others  prior  to  this  work. 

2.2  Tunneling  Hamiltoni an/We rthamer  Theory 

The  original  calculation  of  the  frequency  dependence  of  the 
tunneling  current  in  a  superconducting  tunnel  junction  is  due  to 

Q 

Werthamer.  Werthamer's  calculation  in  turn  is  based  on  the  tunneling 
Hamiltonian  method  used  originally  by  Bardeen,*6  and  by  Cohen,  F^licov, 
and  Phillips,*6  and  the  calculations  of  Ambegaokar  and  Baratu.r,*  and 
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of  Josephson. 

The  Hamiltonian  for  a  tunnel  junction  is  given  by 

H*HL+HR+HT  (2.1) 

where  and  HR  are  the  complete  Hamiltonians  of  the  left  and  right 
superconductors,  and  Hy  is  the  term  which  couples  the  left  and  right 
superconductors. 


The  tunneling  term,  Hy,  is  given  by 


"r  "  &0  [TMckX  +  (TM>*<Ccko]  (2'21 

where  cko,  and  d^o  are  the  destruction  operators  for  the  electrons  in 
the  left  and  right  superconductors,  respectively,  and  cr  is  a  spin  index 
for  the  electrons.  The  number  operators  for  the  left  and  right  metals 
are 


\  *  £  4Cko 
\  -  s  <oV 


(2.3) 


The  tunneling  Hamiltonian  formalism  has  several  important 

1  O 

assumptions  built  into  it.  First,  the  spins  of  the  electrons  are 


unchanged  as  they  are  transferred  across  the  barrier.  Second,  the 
tunneling  of  the  electrons  occurs  instantaneously.  Another  important 


assumption  is  that  the  distribution  functions  describing  the  thermal 
distribution  of  the  quasiparticles  are  equilibrium  distribution 
functions.  Third,  the  voltage  drop  is  pictured  as  occuring  entirely 


within  the  barrier.  And  finally,  the  matrix  elements,  TV  ,  and 

K ,  q 

★ 

T  k  ,  are  equal  by  time  reversal  symmetry. 

If  a  time  dependent  voltage  V(t),  which  can  have  a  dc  term  in  it, 
is  applied  to  the  left  electrode,  sich  that  it  is  positive  relative  to 
the  right  one,  the  Fermi  levels  for  the  left  and  right  electrodes  are 
related  by  uL-yR=-eV.  Here  e=|el,  the  absolute  value  of  the  electric 
charge.  The  Hamiltonian  of  the  left  electrode  is  HL(V)=HL(0)-eVNL, 
where  NL  is  the  number  operator  for  the  left  side,  given  by  Eq.  (2.3a), 
i.e.,  the  voltage,  V(t),  only  modulates  the  energy  of  each  electron  on 
the  left  electrode. 

The  calculation  for  the  current  through  the  junction,  outlined 

below,  assumes  an  interaction  representation  in  which  the  time 

dependence  of  the  operators  is  derived  from  the  unperturbed  Hamiltonian, 

H0=Hl(V)+Hr,  whereas  the  time  evolution  of  the  eigenstates  comes  from 

* 

the  perturbation  term,  Hy,  the  tunneling  Hamiltonian.  Therefore  if.c^  (t) 
is  the  destruction  operator  for  an  electron  on  the  left  side  of  the 
barrier  with  the  voltage  applied,  then  its  time  dependence  can  be 
calculated  from  the  Heisenberg  equations  of  motion;  it  depends  only  on 
the  Hamiltonian  for  the  left  electrode,  H^(V),  and  is. 

"rt'koW  *  [Vt>-lVv)]  <*•«> 

The  solution  of  Eq.  (2.4)  is 
<L(t)  ■  e1*'11'2  c.  (t) 


(2.5) 
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al "  Tv(t) 


(2.6) 


the  effect  of  HL(V)  being  a  modulation  of  the  phase  of  the  destruction 
operator.  The  tunneling  current  operator  is  given  by  the  rate  of  change 


of  the  electron  number  operator 
I(t)  -  -eNL 


(2.7). 


As  for  the  case  of  the  destruction  operators  cko(t) ,  the  rate  of 
change  of  the  number  of  electrons  on  the  left  side  of  the  barrier  is 
calculated  from  the  Heisenberg  equations  of  motion  using  the  total 
Hamiltonian,  H,  in  Eq.  (2.1).  The  only  piece  of  H  which  does  not 
commute  with  is  the  tunneling  Hamiltonian,  Hj,  which  gives  a  non-zero 


commutator: 


\  *  r  K\] 


(2.8) 


The  expectation  value  of  the  current  operator  I(t)  is  calculated  to 
first  order  in  the  tunneling  Hamiltonian,  Hj  using  linear  response 
theory  [or  first  order  perturbation  theory),  and  is  given  by: 


Kt)  =  <Kt)>  =  fRelT  $tdx  en 

kqa  -oo 

X  <[Cka<t>dqo(t)-HT(T)] 


(2.9) 


where  <  >Q  denotes  an  expectation  value  referred  to  the  unperturbed 
Hamiltonian.  The  factor  ent,  comes  from  the  adiabatic  turn  on  of  the 
perturbation  assumed  in  the  perturbation  calculation.  Here,  n-*0+  . 


>  ^  V*  *  -  <  **  *.*  .  '  .  •  C*  V  .  .  *  -  *  *  '  ,  "  •  s.-  S  •  .  *  «  •  .  »  ,1  .  ■  ,  •  ,  •  ,  ■  O  ■*-'  v"  s'  V-  V1  s.  -  i* 


r.i’wr 


Werthamer  evaluated  Eq.  (2.9)  for  a  general  time  dependent  voltage 
V(t)  across  the  barrier,  and  his  result  is  reprinted  below: 

I(t)  =  Im  rfda,da1'{W(a3)W*(a,‘)e’i(a,’a),)t  j  (W'*L  ) 

L_oo  1  c  0 


+  w(«)w(u)' )  *  it  j2(u>.*it0o),J 


(2.10 


In  Uerthamer's  notation,  the  time  dependent  voltage  V(t)  contains  no  dc 

A 

part,  V(t)=V(t)-VQ,  where  VQ  is  a  dc  voltage.  The  functions,  W(ou),  are 

tA 

related  to  exp[(-ie/fi)  j  V(t '  )dt '  ]  by  a  Fourier  transform: 


exp{-Jf  V( t *  )dt ' }  =  5”dco  W(oj)e' 


(2.  n: 


The  phase  <f>  contains  the  time  dependence  from  the  dc  bias  voltage,  VQ, 


♦  =  W  V  =  V 


(2.12) 


and  the  phase  aQ  is  the  initial  phase  difference  between  the  two 
electrodes.  The  frequency-dependent  current  amplitudes  j^(o>)  and 
)  are  defined  by 

jj(ui)  =  2e  S  ^  d“>jdu>2  £f  (u^)-f  (^jj  |Tk^|^ 


x  A^Uj)  Aq(u)2)  Juj-Wg-w+inj  ^ 

e1“^2<“)=  2e  [f'(“l)-f'<"2>]Tk.,T-l<.-< 

x  Bl(wj)  ( ui., )  fui,  —  ui«— ui—  i ril  ^ 


(2.13) 
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where  f“( w )=[e+®w  +1]”* 

is  the  Fermi  distribution  function,  &=(kT)-1  and1i*l.  The  terms 

IQ 

Ak(«),  and  Bk(«),  and  Bk(«)  are  the  spectral  weight  functions  for 
the  superconductor  on  the  left  side 

AkU)  *  \  [(1  ♦  (ek/Ek))5(a.-ek)  +  (i  -  (ek/Ek))  6 (w+ek)j 

Bk(t»j)  a  j  ^(w-Ek)  -  6(w+Ek)j 

6k(u)  *  (^k/Ak)  Bk(w)  (2.14) 

Ek  is  the  energy  of  the  quasiparticle  excitation  and  is  given  by  the 
usual  expression: 

i 

Ek  -  (e*  +  l\)2  (2.15) 

where  Ak  is  the  superconducting  energy  gap  and  £k  is  the  kinetic  energy 
of  the  quasiparticle  relative  to  the  Fermi  energy. 

When  the  temperature  T  is  raised  above  the  superconducting 
transition  temperatures  of  the  two  junction  electrodes  the  spectral 
weight  function  Bk(w  )  -*•  0  since  Bk(u>  )%Ak>  and  Ak(w  )  -*•  6(w-ek),  the 
spectral  weight  function  for  an  electron  in  a  normal  metal. 

The  current  amplitudes  j^(u  )  and  j2(w  )  are  computed  by 
substituting  Eqs.  (2.14)  into  Eqs.  (2.13),  and  doing  the  w ^  and  u2 
integrals.  The  sums  over  wave  vectors  are  replaced  by  integrals.  The 
resulting  integrals  can  be  done  only  by  numerical  integration  at  finite 
temperatures;  however  at  T*0  they  can  be  evaluated  analytically.  Fig. 
2.1  displays  the  real  and  imaginary  parts  of  j 1 w  )  and  j2(  »  )  at  T=0 


FIG.  2.1  Real  and  imaginary  parts  of  the  response  functions  j.(w)  an< 

T=0*  The  peak  in  Retj?(  <»>  )3  at  fi u)  /2a®  1  is  the  Riedel  peak. 
rBLj2(u)  )j  is  the  pair  response  function,  while  Im[j.(u>)]  is  the 
quasiparticle  I-V  curve.  [Adapted  from  Ref.  9.]  * 
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I 

and  for  a  junctions  with  both  electrodes  made  of  the  same 

superconducting  material.  Atfiw/2A*1  there  is  a  peak  in  the  real  part 

of  j2( u)  )  and  jj( u  ).  The  peak  in  Re[j2U  )]  is  the  peak  first 

20 

described  by  Riedel,  while  the  sharp  rise  in  Im[ ( u*  )D  atfiw/2A»l  is 
the  gap  structure  first  seen  in  the  I-V  curve  of  a  tunnel  junction  by 
Giaever.^ 

2.3  Results  From  the  Werthamer  Theory 

The  main  result  of  the  Werthamer  theory,  Eq.  (2.10)  is  rather 
formidable  looking,  however,  we  can  make  sense  of  it  by  analyzing  it  in 
specific  cases  of  interest.  We  begin  with  the  case  of  v(t)=0,  i.e. 
only  a  dc  voltage  on  the  junction.  Eq.  (2.11)  impl ies  W(o>  )*  <5  ( u  ) 

A 

since  V(t)*0.  Inserting  this  value  of  W(u>  )  in  Eq.  (2.10)  we  obtain 

I ( t)  *  ImjjjOy-jj  +  C0SU+ao) 

+  Re  [j2(^]  sin(<j>+a0)  (2.17) 

where  u0*2eV0/fi,  4>=*w0t,  and  aQ  is  the  intial  phase  value,  using  Eq. 

(2.12).  The  real  and  imaginary  parts  of  jj  and  j2,  written  in  Eq. 

(2.17),  can  be  rewritten  in  standard  notation  ’  as  follows. 

W> 1  Im  [V'V*'] 

‘jl'V1  ’  Re  [j2<eV0/fl)]  (MB) 

!J2(V(I,T)  -  Inl  [l2(«o«>] 

where  VVdc*  Eq.  (2.17)  becomes 

Kt)  *  VVT>  +  IJ2(Vo*T)  cos<*  +  °o) 

+  ^(VqJ)  slnU  +  ct0) 


(2.19) 


As  VQ  -*-0,  ^qp(v0»T)  -*-0  and  Ij2(v0»T)  0,  and  we  are  left  with  the 

term  Ijj(0,T)s1n( $+aQi) ,  the  dc  Josephson  effect.  I  and  Ij2  can  be 
written  in  terms  of  conductances, 

V°o'VT>  V0 

IJ2-01(Vo-T)  Vo  <2-20 

so  that  Eq.  (2.19)  becomes 

I(V0,T;t)  *  Ij^VT)  +  {oi^v0*T)  cosfl>(t)+ 

°o(VT)1Vo  <2-21 

where  *  (t)=  <t>(t)+oQ.  Eq.  (2.21)  is  the  result  obtained  by  Josephson^ 
in  1962. 


The  factor  ^(V^T)  at  VQ*0  is  called  the  critical  current.  The 

notation  is  Ijj(0,T)»Ic(T)  and  it  can  be  evaluated^  to  find  a  simple 

form  for  the  case  of  a  tunnel  junction  with  two  identical 

superconductors  for  any  temperature  T,  with  energy  gap  a(T)  and  normal 

resistance  R  : 

n 


IC<T>  ■  Unh^) 


(2.22a; 


For  dissimilar  metals,  at  T«0,  with  gaps  ^(G)*^,  and  a2(0)=a2,  an 
analytic  form  has  been  found  ^>95 


ic(o) 


■  *1*2 

V*2 


)  K(. 


1*1 ~A2^ 

VA2 


) 


(2.22b] 


where  K  is  the  complete  elliptic  integral  of  the  first  kind.  At  finite 
temperatures,  the  expression  for  the  critial  current  of  a  tunnel 


17 


junction  is  written  in  terms  of  a  series  expansion.17  Fig.  2.2(a) 
displays  the  critical  current  as  a  function  of  T  for  a  Sn-Sn  junction, 
while  Fig.  2.2(b)  shows  the  critical  current  as  a  function  of 
temperature  for  a  Sn-Pb  junction.17  As  T-^T  ,  for  the  Sn-Sn  junction 
Ic«(l-t)  where  t«T/Tc. 


It  is  relatively  straightforward  to  find  the  effect  of  a  dc  bias 
voltage  and  an  ac  voltage.  In  this  case  V(t)=VQ+VLcos( u  Lt)  and  the 
Josephson  relation  2eV(t)  * <hd#(t)/dt  implies 
,  /2eV, 


*<t)  -  «¥•>»  -  -L* 


(2.23) 


and 


where  we  have  used  the  relation 


(2.24) 


exp(-ixsin(Q) )  =  z  J  (x)  e 

n=—  n 


-in9 


(2.25) 


where  Jn  is  the  nth  order  Bessel  function.  The  U  function  becomes 


W(w)  *  E  6(w  +  ku.  ) 

1(9-.  *  nuL  L 

Substituting  Eq.  (2.26)  into  Eq.  (2.10)  we  obtain: 


(2.26) 


I(t)  -In  E  Jk(a)Jk,(«)  j1(-k*u,L  +  Jsa,0)  e'i(k'"kHt 
kiks-« 


+  E  J.(a)Jk,(a)  j p ( “ k ' a>.  +sSt*»n )  e^k  +k)wL^+1u>0^+1a0 
kik*—  K  £  L  0 


(2.27) 


where  a*eV^/hu)^  is  the  normalized  ac  voltage  in  the  junction  and 
u0*2eV0/fi.  The  dc  I-V  curve  for  the  dc  plus  ac  voltage  sources  is 
calculated  from  Eq.  (2.27)  and  is: 

I(VJ  *  £  J?(a)  Im{ j . ( - ku, +e V  /h ) } 

0  k*-«D  ^ 

+  £  |  £  'V°^Jn-k^  Re^ j2^k‘Jsn^“L^^5^Vo^n^le< 

njo  k*— 

(2.28) 

Eq.  (2.28)  is  the  fundamental,  experimental ly  verifiable  result  of 
the  frequency-dependent  theory  of  the  Josephson  tunnel  junction.^*  In 
this  equation,  the  imaginary  part  of  ^  is  the  dc  I-V  curve-- 
Im( jx(  w  ))“IdC(‘ft<JJ  /e)*Iqp(fiu  /e) ,  while  the  real  part  of  j2  contains  the 
frequency  dependence  of  the  Josephson  effect,  and  is  rewritten  in  terms 
of  the  response  function  IJ1,  Re[j2(w  )]=I01(fi  u /e) .  Eq.  (2.28)  is 
rewritten  as: 

I(v«)  *  kL  Jk(a)  wk(iK» 

1  kLJk(o)J"-k(o)  Iji((k-,*n)r,-)i  •<»,*<&.» 

(2.29) 

Eq.  (2.29)  contains  two  terms.  The  first  term  comes  from  the 
interaction  of  the  electromagnetic  field  in  the  junction  with  the 
quasiparticles,  and  is  the  mathematical  description  of  how  the 
photon-assisted  tunneling  effect  changes  the  shape  of  the  whole  I-V 

Q 

curve.  This  term  reproduces  the  result  of  the  Tien- Gordon  theory  of 


the  photon-assisted  tunneling  effect.  The  photon- assisted  tunneling 
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steps  or  Dayem-Martin7  steps  occur  at  ±Vgap+k(<fia)  L/e)  as  can  be  seen 

from  Eq.  (2.29)  by  observing  that  a  feature  which  is  significant  on  the 

I-V  with  no  applied  radiation,  i.e.  the  gap  structure  on  the  dc  I-V 

curve  (Fig.  1.1(b))  will  appear  at  the  voltage  k(fi  uj L) /e  away  from  the 

original  structure,  but  modulated  in  amplitude  by  the  squares  of  Bessel 
o 

functions,  J£(a) . 

The  second  term,  unlike  the  first  one,  has  an  effect  only  at 
discrete  voltages,  V0*n(fi  u ^/2e) ,  producing  "spikes"  on  the  I-V  curve  at 
these  voltages.  When  the  junction  is  dc  current  biased  the  "spikes" 
become  flat  regions  on  the  I-V.  This  second  term  describes  the  ac 
Josephson  effect.  Figs.  1.2a  and  1.2b  previously  illustrated  these 
effects,  and  the  reader  is  urged  to  look  at  these  figures  again.  The 
frequency  dependence  of  the  Josephson  effect,  first  obtained  by 
Werthamer,  is  displayed  in  this  term  through  the  dependence  of  the  step 
amplitudes  on  1^,  the  pair  response  function. 

The  frequency  dependence  of  the  Josephson  effect,  observable  as  an 

enhancement  of  the  step  amplitudes  at  the  Riedel  peak,  was  first  tested 
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by  Hamilton  and  Shapiro  in  1971,  using  20-26  GHz  microwave  radiation. 

Their  results  were  however  limited  to  frequencies  only  within  x5%  of  the 

Riedel  peak.  Buckner,  Finnegan,  and  Langenberg^  studied  Sn-SnO-Sn 

junctions  at  135  GHz,  and  also  observed  the  effect  of  the  Riedel  peak, 

but  their  work  was  also  limited  to  frequencies  within  a  few  percent  of 

25 

the  Riedel  peak.  Weitz,  Skocpol ,  and  Tinkham  first  observed  the 
Riedel  peak  in  the  far-infrared  directly,  by  measuring  the  enhancement 
of  the  size  of  the  first  Josephson  step  at  the  peak  and  its  subsequent 


fall -off.  They  used  Nb-Nb  point  contacts  for  this  work.  While  their 

work  showed  an  overall  consistency  with  prior  work,  and  they  observed  a 

fall-off  in  step  amplitudes  with  roughly  the  right  frequency  dependence, 

the  observed  step  amplitudes  were  a  factor  of  two  smaller  than  expected 

from  theory,  even  with  corrections  for  reductions  in  the  step  size  from 

heating  or  noise.  It  is  at  present  an  open  question  whether  there  are 

irreducible  differences  between  tunnel  junctions  and  point  contacts,  due 

to  the  unkown  nature  of  the  contact,  although  significant  progress  has 
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been  made  recently  on  this  question. 

The  Werthamer  result,  Eq.  (2.10),  for  a  general  time-dependent 

voltage  is  very  difficult  to  solve,  and  is  appropriate  only  to  a  voltage 

27 

Diased  situation.  Harris  recast  the  Werthamer  theory  from  the 

frequency  domain  description  of  Eq.  (2.10)  to  a  time  domain 

formulation.  He  developed  general  techniques  for  digitally  computing 

the  behavior  of  a  junction  connected  to  an  arbitrary  circuit.  These 

methods,  however,  require  considerable  computer  time,  and  have  not  been 

used  extensively  in  analyzing  data  on  tunnel  junctions.  A  more  common, 

simpler  approach,  in  wide  use  is  a  model  of  the  junction,  called  the 
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resistively  shunted  junction  (RSJ)  model  developed  by  Stewart  and 
McCumber.  This  model  is  essentially  a  low  frequency  or  low  voltage 
approximation  to  Eqs.  (2.10)  and  (2.19)  and  is  discussed  in  the 
following  section. 


2.3  The  RSJ  Model 

The  RSJ  model  can  be  derived  from  the  microscopic  theory  result  of 
Eq.  (2.19)  by  assuning  first  that  the  pair  response  function  whose  imagi¬ 
nary  part  gives  an  Ij(V,T)  which  is  not  frequency  (voltage)  dependent.  So 
for  Ij(V,T),  the  maximum  dc  supercurrent  at  the  junction  temperature, 

I  (T),  is  substituted.  The  quasiparticle  current,  I„JV,T),  is  modelled 
t  qp 

as  either  a  voltage  dependent  conductance  G(V),  i.e.  Iqp*G(V)V,  or  even 
more  drastically  as  simply  a  resistor  with  the  value  of  the  normal  state 
conductance  ,  Iqp=VGn.  The  cos(<t>  )  term  is  ignored. ^  Finally,  the 
capacitance  of  the  junction  can  be  included  in  the  equations  by  a  term 
«C(dV/dt).  The  external  circuitry  which  drives  the  junction  is  modelled 
as  a  dc  current  source,  plus  a  time  dependent  (ac)  source  if  needed. 

The  equation  which  results  if  these  substitutions  are  inserted  into  Eq. 
(2.19)  is: 

IT  -  G(V)V  +  Ic  sin*  +  C  (2.30) 


where  IC*IC(T),  and  V=(1V2e)d  $ /dt.  The  total  bias  current,  IT,  through 
the  device  is  written  as. 

IT  *  Idc  +  IL  sinu>Lt  +  IN(t)  (2.31) 


where  a  noise  term,  1^,  has  been  included  to  model  thermal  fluctuations 
in  the  resistor,  Rn,  or  an  external  noise  source.  Inserting  the 
equation  for  the  junction  voltage  V=(fi/2e)<£  into  Eq.  (2.30)  ,  the 


(^£-c)  sin<j> 


+ 


(2.32) 


where  Ij  is  given  by  Eq.  (2.31).  Eq.  (2.32)  is  a  very  general  RSJ 
model,  although  not  nearly  as  complex  as  the  original  Werthamer 
equation.,  it  still  retains  solutions  of  considerable  complexity 
depending  on  specific  choices  for  G(V),  C,  and  the  rf  drive  term 
I^sin(w^t).  Fig.  2.3  schematically  illustrates  the  RSJ  model.  Fig. 
2.3(a)  is  a  schematic  of  Eq.  (2.30),  where  the  Josephson  element 
defined  by  the  symbol,  |,  in  the  Figure,  is  shunted  by  the  capacitor  C, 
and  the  conductance  G(V) ,  and  driven  by  the  current  1^.  Fig.  2.3(b) 
shows  that  the  conductance  can  be  pictured  as  the  normal  conductance  of 
the  junction  or  the  non-linear  quasiparticle  conductance,  while  Fig. 
2.3(c)  illustrates  the  total  current  can  be  viewed  as  a  dc  bias  current, 
an  ac  current  source,  and  a  noise  current. 

We  can  easily  determine  the  effect  of  an  ac  voltage  source  on  the 
junction  I-V  curve  using  the  RSJ  model.  We  take  a  voltage  source  with 
dc  and  ac  terms  V(t)=VQ+VLcos( u  ^t) .  The  phase  t  becomes 
rf(t)  =  (2eV0/f  t+(2eVL/1iu  L)sin(w  Lt)+tf0.  We  use  G(V)=Gn,  and  C=0.  Eq. 
(2.30)  becomes 

IT  =  GnV(t)  +  Ic  sin<j>(t)  (2.33) 

The  second  term  of  Eq.  (2.33),  I p= I c s i n  <.  (t),  becomes 

2  eV. 

*p  s  cos(u>0+nuL)t  (2.34) 

where  <p  Q=n/2  to  maximize  the  zero  voltage  critical  current,  and 
standard  trigonometric  identities  were  used  to  arrive  at  the  result. 
Jp(2a)  is  the  Bessel  function  of  order  n,  and  a=eVL/fiaj^  as  before.  The 


Resistively  Shunted  Junction  (RSJ)  Model 
[Current  Biased] 


FIG.  2.3  (a)  Schematic  illustration  of  the  current  biased  RSJ  model. 

I,  is  the  total  current  through  the  junction,  C  is  the  junction 
capacitance,  G( V)  is  the  quasiparticle  conductance,  the  \  denotes  the 
Josephson  element,  (b)  The  quasiparticle  conductance  can  be  linearized, 
i.e.  represented  as  (Rn)~  ,  or  left  non-linear,  (c)  The  total  current, 
IT  through  the  junction  is  the  sum  of  a  dc  term,  I  .  ,  an  ac  drive  term, 
I^sinu^t,  and  a  noise  term,  1^.  J!" 
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time-averaged  dc  current  becomes 


<IT>  =  6nVQ  +  Ic£  Jn(2ct)  6(ui0+nu)L)  (2.35) 

n=-°° 

since  the  second  term  of  Eq.  (2.34)  contributes  only  if  u)0=-nwL  where 

n  is  an  integer.  Therefore*  "steps"  appear  on  the  I-V  curve  at  dc 

voltages  which  are  integer  multiples  of  as  observed  before.  In 

this  case  the  full  width  of  the  nth  step  is  I  *2I_J„(2a) ,  instead  of  the 

n  c  n 

complicated  structure  of  Eq.  (2.29).  The  second  term  of  Eq.  (2.29) 
reduces  to  the  second  term  of  Eq.  (2.35)  in  the  limit  of  fito  L/2e  -*■  0. 

In  this  case  Ij((k-n/2)fi  u  L/e)%Ic,  and 

I  Jk(a)Jn_k(a)=Jn(2a), 

so  that  the  term  is 

I  Z  J  (2o)  s(V  ±n(*fiu> . /e)). 
c  n}0  n  0  L 

"  ? 

The  term  Z  Jk^a^qp^Vo+lc^  u  i/e)  reduces  t0  v0Gn  the  I_v 

k=-» 

curve  is  linear,  since  the  term  (V  +k(ft  u^/e))  is  compensated  for  by  a 
term  in  the  sum  with  (VQ-k(fi u  L/e) ) ,  leaving 

Jk2(“>[VoGn]  *  Vn  <2-36) 

•  o 

since  Z  Jr(a)sl.  The  resulting  equation  is: 

n= -• 

,  “  f»U), 

I(vo>  *  GnVo  +  'CS  «(*0t "<-2T» 

n>o 

and  is  equivalent  to  Eq.  (2.35). 


(2.37) 


This  discussion  demonstrates  how  the  Werthamer  theory  reduces  to 
the  RSJ  model  only  in  the  limit  of  low  frequencies.  This  point  is 
important  since  it  indicates  that  for  junctions  which  are  sensitive  to 
frequencies  approaching  the  energy  gap  frequency,  caution  must  be 
applied  in  using  results  from  the  RSJ  model  when  interpreting 
experimental  data. 

We  now  discuss  the  effects  of  noise  sources  intrinsic  to  the 
junctions  on  their  I-V  curves  for  both  the  transition  out  of  the 
zero- voltage  state  assuming  a  dc  current  bias,  and  the  rounding  of  the 
ac  Josephson  steps. 

2.4  Noise  Effects 

The  intrinsic  noise  in  electronic  devices  such  as  transistors, 

resistors  and  superconducting  tunnel  junctions  is  described  by  the 

power  spectrum  for  current  or  voltage  fluctuations.  In  superconducting 

tunnel  junctions  the  noise  could  come  from  fluctuations  in  the 

quasi  particle  (or  pair)  current,  for  example.  The  power  spectrum  for 

quasiparticle  current  fluctuations,  at  a  constant  dc  bias  voltage  VQ, 
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first  calculated  by  Dahm  et  al .  ,  is  given  by 

PI<“>  1  h  <  VVT>  «th(f<eV0«M) 

♦  ‘qp^o-T1)  coth(|<eV0-M)}  (2.38) 

where  Iqp(VQ)  is  the  quasiparticle  current  in  the  junction,  defined  in 


Eq.  (2.18).  The  power  spectrum  is  defined  as  the  Fourier  transform  of 
the  current-current  auto-correlation  function: 

QO 

Re  <1  (t+r )  I (t )>  a  _/  Pj(u>)  cos  uxdu  (2.39) 

where  the  autocorrelation  function  is 

< I ( t+T )  T[t)>  .  /  <IU+t)  !(t)>  dt  (2.40) 


In  this  notation  the  angular  brackets  imply  a  statistical  and  quantum 
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mechanical  average,  and  the  bar  indicates  a  time  average.  The 
mean-square  current  fluctuation  (noise),  <ijjj>,  can  be  calculated  from 
the  power  spectrum  by  integrating  over  all  frequencies 

<iN2>  *  f"  PjfuOdu  (2.41) 

The  power  spectrum,  Eq.  (2.38),  while  rather  complicated  can  be 
understood  simply  in  certain  limiting  cases.  If  the  bias  voltage,  VQ, 
is  small  compared  to  kT/e,  i.e.  eVQ<<kT,  and  the  frequency  is  small, 
1W<<kT,  eVQ,  the  power  spectrum  reduces  to 

PI(u)  *  £<4kT>  -3v(V°i  (2.42) 

0 

Eq.  (2.42)  reduces  to  the  standard  Johnson-Nyquis*  result 

p,(v)  .  4kT/Rn  (2.43) 

by  converting  the  frequency  measure  from  radians  to  Hz  using 

PT(\>  )=2hPt ( w  ),  and  the  fact  that  I„.,(V„)/V  =(R„)-1 ,  if  the  junction  is 


in  its  normal  state  (linear  I-V).  If  the  junction  is  below  the 
transition  temperatures  of  the  two  electrodes,  its  I-V  curve  is 
non-linear.  At  low  dc  voltages,  i.e.,  VQ  -*-0,  Iqp(V0MdI/dV)V0,  and 
*qp^o^V^Rd)~*  where  ^“(dl/dV)"1 ,  the  dynamic  resistance  of  the 
junction  at  low  voltages. 

At  high  bias  voltages,  eV0>>kT,iiuj,A,  regardless  of  whether  or  not 
the  junction  is  superconducting,  the  power  spectrum,  Eq.  (2.38), 
reduces  to 

pi'“)  *  f  V V  (2-44) 

in  agreement  .with  the  usual  expression  for  shot  noise,  Pj(v  )*2eIqp(V0) , 
after  conversion  to  frequencies  in  Hz. 

The  mean-square  voltage  noise  <vj">  is  related  to  the  current  noise 
by  the  dynamic  resistance,  Rd: 

00 

<V*>  *  Rd  <i*>  *  Rd  /  PjMdw  (2.45) 
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Eq.  (2.38)  was  used  first  by  Dahm,  et  al .  to  calculate  the  linewidth 
of  the  Josephson  radiation  emitted  by  a  tunnel  junction,  in  a 
calculation  based  on  FM  noise  theory. 

The  power  spectrum  of  Eq.  (2.38) ,  while  an  important  result, 
because  it  is  derived  with  the  assumption  of  a  constant  dc  bias  voltage, 
is  of  limited  usefullnes  in  studying  the  detailed  dynamics  of  a  tunnel 
junction.  In  the  usual  experimental  configuration  the  junction  is  dc 
current  biased  with  noise  and  without  rf,  or  dc  current  biased  and  rf 


voltage  or  current  biased.  In  practice  the  RSJ  model  and  the  limit  of 
Eq.  (2-38)  for  low  voltages  has  been  used  as  an  approximate  scheme, 
with  reasonable  success,  for  calculating  the  effects  of  noise  on  a 
junction.  In  the  low  capacitance  limit,  b.  -*•  0,  appropriate  to  shunted 
tunnel  junctions  or  metallic  weak  links,  analytic  solutions  to  the  RSJ 
model  with  a  white  noise  source  such  as  for  the  power  spectrum  of  Eq. 
(2.43),  are  available.  We  will  now  give  a  brief  description  of  the 
results  of  Ambegaokar  and  Halperin  (AH)^4,  for  ec  -*■  0,  and  Lee,^5,  for 
BC<1. 

Eq.  (2.32)  can  be  simplified  by  rewriting  it  in  terms  of 

normalized  units.  Frequencies  are  normalized  to  where 
1/2 

uij»(2eIc/fiC)  is  the  junction  plasma  frequency.  Currents  are  divided 
by  I  ,  the  critical  current,  and  voltages  by  the  I  R  product,  where, 
for  example,  IcRn»ira(0)/2e  for  a  junction  with  equal  gaps  at  T=0.  The 
resulting  equation  is 
2 

+  g(^J-)  +  sin<j>  *  p^c  +  oL  sinftx  +  pn  (2.46) 

dx 

where  x=t  Uj  is  the  normalized  time,  g(d^/dx )®Iqp( ('1^/2e)d0/dx )  )/Ic  is 

a  function  of  the  normalized  frequency  (voltage),  d<*>/dx .  The 

normalized  ac  current  and  frequency  are  p  I 1 c ,  and 

respectively.  The  noise  term  is  p  N= I N/ Ic .  If  G(V)  =  (Rn)‘1  then 

-1/2 

g(d$/dx)  becomes  g(d$/dx)=(Bc)  d^/dx  which  is  the  situation 

2  2 

described  by  AH,  and  Lee.  The  parameter  B^=(u,R„C)  =2eKR_C/‘fi  is 

con  c  n 

called  the  Stewart-McCumber  parameter,  and  is  used  to  describe  the 
hysteresis  on  the  I-V  curve  of  the  junction.  The  s  parameter  is  small, 
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8C«1,  in  the  low  capacitance,  low  hysteresis  limit  used  by  AH,  and  Lee. 

a 

The  noise  term  p  N>  is  described  by  its  autocorrelation  function  in 
the  time  domain: 


<pn(t+t')  pm(t')> 


N' 


Y8C72 


6(t) 


(2.47a) 


where  B  =(2el  R„C)/fi,  and  Y=fil  /ekT,  is  calculated  from  the 

v  w  ll  V* 


unnormalized  autocorrelation  function 
<iN(t+x)  IM(t)>  *  tP-  <s(t) 


(2.47b) 


appropriate  to  the  resistance  of  the  junction  in  the  low  voltage  Johnson 
noise  limit.  For  an  nonlinear  I-V.  curve  in  the  low  voltage  limit,  Bc  is 
replaced  by  Bc'*(2elcR^C)/fi  and  Y  remains  the  same  in  Eq.  (2.47a),  or 
equivalently  Rn  is  replaced  by  Rd  in  Eq.  (2.47b).  All  terms  on  the 
left  hand  side  of  Eq.  (2.46)  remain  the  same  since  we  have  used  the 
general  form  for  the  normalized  quasi  particle  current,  g=g(d4>/dx) . 


The  RSJ  model  can  be  rewritten  in  terms  of  the  mechanical  model  of 
a  particle  in  a  tilted,  periodic  potential.  If  we  consider  the  case 
where  the  conductance  G(V),  is  simply  Gn=(Rn)"1.  then  the  middle  term  of 
the  left  hand  side  of  Eq.  (2.46),  g(d$/dT),  is  (Bc)"ly,2[0].  We  take 

A 

for  the  moment  p  n=0,  then  Eq.  (2.46)  becomes: 


where  U(<J)s(- p  dc0-cos0)  is  the  tilted  periodic  potential  the 
particle  is  sliding  in.  A  particle  is  trapped  in  a  well  at  position  A 
in  Fig.  2.4.  Friction  or  damping  is  modelled  by  the  second  term  of  Eq. 


Particle  in  a  Tilted  Periodic  Potential 
(Junction  Analogy) 


B 


0  2^  4  *  6  t 

<t> 


FIG.  2.4  Mechanical  analogy  of  the  RSJ  model  with  a  particle  moving  in 
a  tilted  periodic  potential.  As  the  dc  current  is  increased  from  zero 
the  particle  trapped  in  a  well  at  A  can  escape  prematurely  from  the 
well  by  jumping  over  the  energy  barrier  E.  and  begin  sliding  down  the 
potential.  [After  Ref.  37.] 


(2.48).  With  no  noise  the  particle  remains  trapped  in  the  well  until 
the  tilt  is  great  enough  for  the  particle  to  become  dislodged  and  roll 
down  the  hill.  As  the  tilt  decreases,  the  particle  continues  rolling 
until  the  tilt  reaches  a  minimum  value  which  depends  on  the  frictional 
term.  If  the  friction  is  large,  i.e.  low  6C,  then  the  particle  comes 
to  rest  when  the  tilt  is  larger  than  it  would  for  the  low  friction, 
large  ec,  case.  This  means  that  the  hysteresis  of  the  junction  I-V  is 
less  for  the  low  &  junction  than  for  the  high  a  junction.  Noise  can 
cause  the  particle  to  roll  out  of  the  well  at  a  lower  tilt  than  it  would 
without  noise.  Noise  can  kick  the  particle  out  of  the  bottom  of  the 
well  at  A,  over  the  barrier  B,  causing  the  particle  to  roll  down  the 
hill.  The  phase  then  changes  with  time,  or  the  junction  has  a  voltage. 

Ambegaokar  and  Hal  peri n  used  Eq.  (2.46)  in  the  limit  of  pl=0,  and 
Bc«l.,  and  g(d«j»/dT)*(Bc)"1/,2(d(|»/dr).  AH  found  that  the  critical 
current  was  depressed  and  the  transition  out  of  the  zero  voltage  state 
was  rounded  as  the  parameter  Y*tiIc/ekT  was  decreased  from  large  values 
to  zero,  as  is  shown  in  Fig.  2.5.  The  parameter  y  measures  the 
strength  of  the  coupliny  of  the  junction  from  the  Josephson  effect, 
against  kT,  the  thermal  energy  tending  to  destroy  the  phase  coherence  of 
the  junction.  We  can  get  a  feeling  for  the  significance  of  various 
values  of  y  if  we  compute  it  for  real  parameters.  Inserting  values  of 
e,  fi,  k  into  the  equation  for  y  ,  we  obtain  Y=47.68[Ic(uA)/T(K)].  If 
Ic*  1  uA  in  the  absence  of  fluctuations,  and  T=  1  K,  then  y^.  From  Fig 
2.5  we  observe  that  the  critical  current  would  be  depressed  by  -*.20%  by 
the  thermal  fluctuations.  In  the  limit  of  Y-*0,  the  noise-rounded  I-V 


Current  (reduced 


FIG.  2.5  Effect  of  thermal  noise  on  the  RSJ  model  I-V  curves  in  the 
limit  of  0  »(2el jrC)/f» -►  o.  at  different  values  of  the  noise 
parameter.  y=fil  7ekT.  [After  Ref.  34.] 
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curve  of  the  junction  would,  in  fact,  be  linear.  The  critical  current 
(Josephson  effect)  would  be  completely  destroyed  by  the  thermal  noise. 


P.  A.  Lee's  result  is  more  general,  and  is  written  below  for  the 
situation  in  which  the  junction  is  driven  by  an  ac  signal,  for  the  nth 
Josephson  step: 


4V  .  4-  Rdl°  (l-e”*'-)  T1-1(1+Bc(T2/T1))  (2.49) 

where  the  arguments  of  this  equation  are  as  follows.  aV  is  the  dc 
voltage  measured  from  the  step  center,  o*Al/I°,  where  aI  is  the  current 
along  the  step  measured  from  the  step  center,  and  1°  is  the  step 
half-width  without  noise  and  Rd  is  the  dynamic  resistance  at  the  step 
center  without  radiation.  The  current  and  voltage  at  the  step  center 
are  IQ  and  VQ  respectively, Y  ’=Y(V0/I0Rd)  where  Y=fiI°/ekTeff ,  and 
finally  e^*(2eI°RdC)/fi.  Tj  and  T2  are  integrals  given  by 


T1 

T2 


o*d$  IQ  [v'sinU/2)]  exp(-|Va*) 

Q  d<j>  sin(<j>/2)  ' sin(4>/2)j  exp(-^-Y* ai4> ) 


(2.50) 


where  IQ  and  1^  are  modified  Bessel  functions.  Eq.  (2.49)  in  the  limit 
of  Bc*0,  agrees  with  the  AH  result  if  a=I/Ic,  y'  =Y=fiI c/ekT ,  Rd=Rn,  and 
ec*(2eIcRj;C)/ft,  etc. 
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Henkels  and  Webb,  and  Weitz  used  a  simplified  version  of  Eq. 
(2.49),  appropriate  for  8C=0,  in  their  analyses  of  the  noise  rounding  of 
Josephson  steps  at  microwave  and  FIR  frequencies,  respectively,  on  point 
contacts.  In  Chapter  IV,  we  will  describe  fits  to  step  shapes  on  high 
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resistance  tunnel  junctions,  using  Eq.  (2.49),  the  Lee  result.  We  will 
also  describe  a  more  complete  analysis,  using  a  computer  simulation  of 
the  junction  dynamics,  of  both  the  noise  rounding  of  the  Josephson  steps 
and  the  reduction  of  the  critical  current  from  noise,  based  on  the  model 
Of  Eq.  (2.46). 

Lee  also  calculated  the  effect  of  noise  rounding  in  the  limit  of 
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large  Bc.  His  result  was  used  by  Fulton  and  Dunkleberger  in  their 
studies  of  the  lifetime  of  the  zero-voltage  state  of  Sn-Sn  and  Pb-Pb 
tunnel  junctions.  The  inverse  lifetime  or  probability  per  unit  time  of 
switching  is 

t'1  *  ( uTj/2-it )  exp(-E/kT)  (2.51) 

where  E»(‘fiIc/2e)[a(2sin”*a-it)+2cos(sin"*a)]  is  the  energy  barrier  in  the 
presence  of  current, 

with  ofI/Ic,  and  u>j(l-c^)^it the  angular  attempt  frequency.  If 
the  bias  current  is  swept  at  a  rate  dl/dt,  Eq.  (2.51)  can  be  used  to 
compute  the  probability  per  unit  current  P(I)  of  the  junction  switching. 

P(I)  *  (T(I)^r1  exp{-/  dl'Ud')^)"1}  (2.52) 

This  probability,  measured  by  Fulton  and  Dunkleberger,  was  used  by 
inversion  of  Eq.  (2.52)  to  compute  the  lifetime  t(I)  for  their  tunnel 
junctions.  They  found  very  good  agreement  with  the  expression  for 
t(I)  Eq.  (2.51)  for  their  junctions  in  the  temperature  range  from 


4.2K  down  to  1.5K. 


2.5  Discussion 


Although  the  theory  presented  in  Sections  2. 2-2. 4  is  rather  well 
established  and  considerable  experimental  work  has  been  done  to  confirm 
it  in  one  aspect  or  another,  as  pointed  out  in  the  text,  much 
experimental  work  and  theoretical  work  remains  to  be  done,  particularly 
at  the  limits  of  very  low  (^1  mK)  temperatures  for  the  transition  out  of 
the  zero  voltage  state,  and  at  high  frequencies  (FIR)  for  the  Werthamer 
theory  and  the  noise  rounding  of  the  Josephson  steps  at  finite  voltages. 
We  will  begin  this  section  by  discussing  the  problem  of  the  transition 
out  of  the  zero  voltage  state. 

As  the  temperature  of  a  tunnel  junction  is  lowered  into  the  mK 
range  it  is  expected  that  the  thermally  activated  process  described  by 
Eq.  (2.51)  will  be  frozen  out.  Instead  as  the  current  is  swept  towards 
the  critical  current,  the  transition  out  of  the  zero-voltage  state  is 
expected  to  occur  by  quantum  mechanical  tunneling  through  the  energy 
barrier.  Voss  and  Webb  performed  experiments  on  high  e  (s  ^50-5000, 

V  w 

2 

or  low  current  density,  J  -v.16-160  A/cm  )  all -Nb  junctions  and  found 
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basically  good  agreement  with  the  theory  of  Caldeira  and  Leggett  ,  but 
they  found  a  value  of  5  for  a  parameter  which  should  have  been  of  order 
1.  Like  many  all  Nb  junctions  their  junctions  had  a  significant  normal 
conductance  below  the  gap,  which  should  not  occur  in  a  tunnel  junction 
which  follows  the  behavior  described  by  the  Werthamer  theory.  In  high 
current  density  junctions  (Jc^5xl05  A/cm2),  Jackel  et  al.*°  found 
evidence  for  macroscopic  quantum  tunneling  at  temperatures  from  2.2K 
down  to  1.6K.  These  experiments  suggest  that  macroscopic  quantum 


tunneling  effects  are  probably  not  important  for  the  junctions  used  in 

3  2 

this  experiment  ( 0_-v6 ,  J.-wlO  A/cm  )  in  the  temperature  range  (T\1.4K) 

W  w 

used,  but  these  effects  may  become  important  on  higher  current  density 
junctions  which  will  be  used  in  future  experiments.  It  would  be 
interesting  to  find  out  if  Sn-SnO-Sn  or  Pb-PbO-Pb  junctions  would  have  a 
closer  agreement  with  the  theory  of  Caldeira  and  Leggett,  since  these 
junctions  have  dc  I-V  curves  which  are  nearly  ideal,  i.e.,  follow  the 
behavior  expected  from  tunneling  theory. 

25 

At  high  frequencies,  Weitz  et  al .  observed  the  fall-off  of  the 

ac  Josephson  effect  above  the  Riedel  peak.  They  found  the  shape  of  the 

fall -off  was  in  agreement  with  the  theory  but  that  the  theory  predicted 

the  step  widths  to  be  a  factor  of  2  larger  than  the  observed  step  widths 

even  after  correction  for  heating  and  noise  rounding  effects.  A 

comprehensive  theory  for  the  dc  I-V  curve  of  point  contacts  has  not  been 

found,  so  the  fact  that  the  step  widths  were  smaller  than  expected  may 

have  more  to  do  with  the  nature  of  the  point  contact  than  any  problem 

with  the  Werthamer  theory,  especially  considering  the  reasonable 

22 

agreement  found  near  the  Riedel  peak  by  Hamilton  and  Shapiro,  and 
23 

Buckner  et  al.  ,  working  with  classical  tunnel  junctions. 

12  41 

Weitz  et  al.  ’  found  Josephson  steps  at  high  voltages  on  high 
resistance  junctions  to  be  rounded  by  noise.  They  were  able  to  explain 
the  observed  rounding  in  terms  of  thermal  noise  at  a  temperature 
elevated  above  the  bath  temperature,  essentially  caused  by  Joule  heating 
of  the  point  contact.  The  temperatures  found,  however,  were 
considerably  (a  factor  of  2)  lower  than  the  temperature  that  would 


correspond  to  the  shot  noise  predictions  of  the  noise  theory  of  D'hm  et 
32  42 

al .  Vernet,  Heneaux,  and  Adde  measured  Josephson  step  linewidths  as 

a  function  of  dc  voltage  for  point  contacts  and  found  good  agreement 

32  32 

with  the  theory  of  Dahm  et  al .  At  low  dc  voltages,  Dahm  et  al .  found 

good  agreement  with  their  theory  from  linewidth  measurements  on  Sn-Sn 

43 

and  Pb-Pb  tunnel  junctions.  Crozat,  Vernet,  and  Adde  found  good 

32 

agreement  with  the  theory  of  Dahm  et  al .  on  microbridges  using  a 
direct  noise  measurement  technique,  but  direct  noise  measurement 
techniques  on  point  contacts  performed  by  Claasen,  Taur,  and  Richards4^3 
showed  considerably  less  noise  than  expected  from  the  tunnel  junction 
noise  theory  of  Ref.  32. 

It  is  clear  from  this  discussion  that  considerable  experimental 
work  needs  to  be  done  to  clear  up  these  discrepancies.  It  is 
particularly  important  that  well  characterized  tunnel  junctions,  for 
which  the  theory  of  Sections  2. 2-2. 4  is  appropriate,  be  used.  The 
experiments  described  in  the  Chapters  which  follow  are  on  Sn-SnO-Pb 
junctions  having  good  gap  structure,  small-area,  and  moderate  current 
densities.  These  junctions  were  made  by  photolithographic  techniques 
which  were  not  yet  developed  when  the  experiments  described  above  were 
performed.  The  experiments  of  this  work  attempt  to  address  the 
questions  raised  about  the  size  of  the  ac  Josephson  steps  and  the 
intrinsic  shot  noise  in  the  junction. 


CHAPTER  III 


EXPERIMENTAL  TECHNIQUES 


3.1  Introduction 

During  the  course  of  this  work  three  types  of  tunnel  junction 

samples  were  made.  The  first  samples  were  small-area  junctions  designed 

for  experiments  on  gap  suppression  induced  by  high  quasiparticle  current 

injection,  and  consisted  of  Sn-SnO-Pb  overlap  tunnel  junctions 
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fabricated  by  the  method  of  Dolan  and  Dunkleberger  .  This  work  will 
be  described  in  Appendix  B.  The  second  set  of  samples  were  small -area 
junctions  fabricated  at  the  centers  of  dipole  antennas  for  use  in  the 
FIR  measurements  which  make  up  the  bulk  of  this  work.  These  junctions 
were  made  primarily  by  the  "resist-aligned"  scheme  of  Howard,46  et  al., 
although  some  samples  were  fabricated  by  the  "edge-aligned”  method.  The 
third  type  of  samples  were  fabricated  at  Sperry  Research  Labs.  The 
Sperry  samples  consisted  of  Nb-Si-Nb  SNAP47  (Selective  Niobium 
Anodization  Process)  junctions  at  the  centers  of  dipole  antennas.  Our 
results  on  the  Sperry  samples  are  reported  in  Appendix  B. 

Because  the  fabrication  of  these  samples  is  rather  involved  we  will 
devote  most  of  this  chapter  to  this  problem.  We  begin  by  describing  the 
mask  making  techniques  used  for  the  first  two  types  of  samples,  and  will 
then  discuss  the  three  methods  of  junction  making  used  at  Harvard.  In 
order  to  keep  the  discussion  to  a  reasonable  length  we  will  leave  the 
complex,  but  critically  important,  details  of  the  substrate  cleaning 
procedures,  mask  and  sample  fabrication,  and  etches  used  to  Appendix  A. 
Following  this  we  will  discuss  the  problem  of  coupling  the  FIR  laser 


radiation  to  the  samples  and  the  particular  method  we  used.  We  will  end 
the  chapter  with  a  brief  description  of  the  FIR  laser  system,  its 
operation,  and  performance. 


3.2  Samp! e  Fabrication 
3.2.1  Mask  Fabrication--Commercial  Masks 

The  first  set  of  samples  were  made  using  a  mask  fabricated  by  a 
local  firm  (Advance  Reproductions,  North  Andover,  MA).  This  mask  was 
used  in  a  contact  printing  procedure  to  produce  photoresist  stencils, 
which  in  turn  were  utilized  in  the  thin  film  evaporation/oxidation 
sequence  used  to  fabricate  the  actual  samples. 

The  mask  made  by  Advance  Reproductions  contained  the  patterns  for  a 
set  of  tunnel  junctions  at  the  centers  of  superconducting  strips  of 
varying  lengths,  designed  to  study  the  non-equilibrium  effects  in  the 
superconducting  strips  created  by  the  presence  of  the  high-current 
density  tunnel  junction  at  its  center.  The  mask  can  be  best  described 
by  imagining  viewing  it  under  a  microscope  as  the  magnification  is 
increased.  With  no  magnification  six  strips  of  chrome  0.3  mm  wide  and 
4.8  mm  long  are  seen.  There  are  four  notches  on  five  of  the  strips  and 
five  notches  on  one  of  them.  Three  of  the  long  strips  contain  tunnel 
junction  patterns,  while  the  other  three  have  long  microbridge  patterns. 
A  close-up  of  a  notch  area  of  the  strip  (the  white  areas  between  the  two 
large  dark  areas)  with  five  notches  is  shown  in  Fig.  3.1(a).  The  dark 
areas  in  Fig.  3.1  are  covered  by  chrome,  while  the  white  areas  are 


clear  on  the  mask.,  The  central  portion  of  Fig.  3.1(a)  is  magnified  in 
Fig.  3.1(b).  The  region  of  the  Figure  with  the  small  white  tab  and  the 
small  black  bridge  defines  where  the  tunnel  junction  will  be  in  the 
completed  sample.  Fig.  3.1(c)  displays  a  pattern  for  a  tunnel  junction 
at  the  center  of  a  superconducting  strip,  taken  with  the  same 
magnification  as  Fig.  3.1(b).  A  long  microbridge  pattern  is  shown  in 
Fig.  3.1(d).  This  mask  was  made  using  an  Electromask  Pattern 
Generator/ Image  Repeater  system.  The  desired  sample  geometry  was 
decomposed  into  a  set  of  rectangles  of  varying  center  positions  (x,y) 
and  widths  and  heights  (w,h).  The  widths  and  heights  of  these 
rectangles  had  to  be  greater  than  a  minimum  value  of  5  um,  and  less  than 
a  maximum  value  of  1.5  mm.  These  limits  on  the  rectangle  dimensions 
depended  on  the  particular  pattern  generator  we  used.  Large  areas  which 
were  to  be  clear  on  the  mask  were  then  decomposed  into  many  rectangles 
of  this  maximum  size.  The  results  of  this  procedure,  the  set  of  numbers 
(x,y,w,h)  for  the  whole  pattern,  were  recorded  on  9  track,  800  bpi  (bits 
per  inch)  tape  in  Extended  Binary  Coded  Decimal  Interchange  Code 
(EBCDIC).  This  tape  was  produced  at  Lincoln  Labs  using  the  pattern 
generating  software  available  there. 

The  pattern  generator  works  by  taking  a  soda-lime  glass  substrate 
coated  with  chrome  and  photoresist  and  projecting  rectangles  of  varying 
widths  and  heights  onto  it  at  locations  defined  by  the  (x,y)  coordinates 
on  the  tape  by  using  a  precision  x-y  translation  stage.  The  resulting 
mask,  after  developing  the  photoresist  and  etching  the  chrome,  was  a  lOx 
mask  (  or  reticle)  if  the  minimum  dimensions  for  the  final  mask  were  to 
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FIG.  3.1  (a)  Section  of  the  commercial  mask  containing  a  tunnel 

junction  pattern  and  a  pad  pattern.  See  text  for  details.  The  black 
regions  have  Cr  on  them  in  the  mask,  (b)  A  close-up  of  the  central 
region  of  (a)  showing  the  tunnel  junction  pattern,  the  narrow  black 
region  in  the  center  will  be  used  to  make  the  photoresist  bridge  which 
is  used  to  form  the  junction,  (c)  Another  section  of  the  mask,  this 
time  with  a  tunnel  junction  pattern  with  long  leads,  (d)  Section  of 
mask  containing  a  long  microbridge  pattern.  The  width  of  the  white  line 
is  only  \Z  uni. 
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be  below  5  ym.  Otherwise  this  mask  would  have  been  the  master  mask. 

The  minimum  final  dimension  we  used  was  2  ym,  so  we  had  the  reticle 
projected  (and  reduced  by  10  in  size)  onto  another  photomask  blank  to 
produce  the  final  master  mask.  We  made  working  copies  by  contact 
printing  from  this  master  mask.  These  copies  were  later  discarded  when 
damaged. 

The  mask  fabrication  procedure  outlined  above  is  commonly  used  in 
the  present  day  electronics  industry.  Because  of  the  need  to  produce  a 
magnetic  tape,  which  requires  considerable  software,  and  because  of  the 
high  cost  per  mask,  we  developed  a  method  for  making  our  own  masks  at 
Harvard.  Our  method  is  based  on  relatively  simple  techniques  that  were 
used  in  the  electronics  industry  in  the  1960's  but  which  were  abandoned 
after  the  invention  of  reliable  computer-controlled  pattern-generating 
equipment. 

3.2.2  Mask  Fabrication  at  Harvard 

All  of  the  fabrication  work  at  Harvard,  reported  here  below,  was 
performed  in  a  clean  room  in  which  all  personnel  using  the  room  were 
required  to  wear  disposable  talc-free  gloves,  and  lint-free  gowns.  All 
the  photolithographic  work  was  performed  under  laminar-flow  hoods. 
Surfaces  not  under  these  hoods,  which  could  collect  dust  and  grease, 
were  wiped  off  with  ethanol  before  work  was  begun  every  day. 

Masks  were  fabricated  at  Harvard  by  a  sequence  of  steps  outlinea  in 
Fig.  3.2.  We  began  with  sheets  of  clear  acetate  upon  which  pieces  of 
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black  tape  were  placed  where  we  wanted  metal  in  the  completed  sample. 
Three  sets  of  acetate  sheets  were  used  because  the  field  of  view  of  the 
microscope  diminishes  as  the  magnification  is  increased.  For  the 
largest  areas  on  the  sample,  a  ?5x  scale  pattern  was  made  for  the  pad 
mask.  A  second  acetate  sheet  plus  tape  pattern,  prepared  at  a  250x 
scale,  defined  the  FIR  antenna  and  dc  connections  to  the  pads.  The 
third  pattern  was  at  a  1250x  scale  and  contained  the  junction  and  its 
connection  to  the  antenna. 

The  patterns  prepared  in  this  way  were  taped  to  a  lightboard  which 
illuminated  them  from  the  back  in  order  to  produce  the  highest  possible 
contrast  between  the  light  and  dark  areas  on  the  pattern.  The  patterns 
were  then  photographed  with  a  high  quality  35  mm  camera  on  high-contrast 
fine-grained  Kodak  Kodal ith  Ortho  film  type  6556.  The  resulting  film 
(emulsion  mask)  had  the  features  of  the  original  acetate  patterns  but 
reduced  in  size  by  a  factor  of  25.  The  pad  mask  was  therefore  at  its 
final  size  and  was  contact  printed  onto  a  photoresist  plus  chrome  coated 
flexible  photomask  blank  (dimensions  2"x2"x0.010"  ).  This  mask  in  turn 
was  used  to  contact  print  the  pad  pattern  onto  a  thicker  chrome 
photomask  blank  (2.5"x2.5"x0.060")  which  after  a  few  more  steps  became 
the  master  mask.  The  emulsion  pad  mask  could  be  contact  printed 
directly  onto  the  master  mask  blank;  however,  the  life  of  the  emulsion 
mask  was  extended  because  it  was  only  used  occasionally,  i.e.  when  it 
was  contact  printed  onto  the  flexible  mask  blank. 

The  lOx  antenna  emulsion  mask  and  50x  junction  emulsion  mask  were 
contact  printed  onto  l"xl"x0.030"  projection  mask  blanks.  These 


projection  masks  fit  onto  the  x-y  stage  of  the  projection  attachment 
to  the  lab's  Zeiss  Universal  Microscope.  After  the  pads  were  contact 
printed  onto  the  thick  photomask  blank,  the  antenna  mask  pattern  was 
projected  onto  the  pad  mask  at  the  proper  location;  the  exposed  mask 
was  then  developed  and  etch  ad.  Because  the  widths  of  the  developed 
regions  Increased  with  continued  development,  we  found  that  stripping 
the  mask  of  the  old  photoresist  and  recoating  with  0.3  ym  thick 
photoresist  produced  the  best  results.  Using  this  freshly  recoated  mask 
with  the  pads  plus  antenna,  the  junction  was  next  projected  onto  the 
center  of  the  antenna  region.  If  the  junction  area  looked  good  after 
the  development,  the  mask  was  etched.  This  mask  was  copied  onto  thin 
flexible  photoresist  coated  chrome  photomask  blanks.  These  flexible 
masks  in  turn  v/ere  used  to  produce  the  photoresist  stencils  which  were 
used  to  make  the  tunnel  junctions. 

3.2.3  Photoresist  Stencils 

Photoresist  stencils  were  fabricated  following  the  procedure  of 

Dolan^*  and  Dunkleberger^5.  Initially  1.5  ym  of  Shipley  AZ  1350J 

photoresist51  was  spun  at  5500  RPM  onto  a  clean  I"x0.5"x0.010"  thick 

z-cut  crystal  quartz  substrate.  The  substrates  were  then  baked  in  a  hot 

plate  oven  for  25±5  minutes  at  80°C.  After  cooling  for  10-15  minutes  in 

the  laminar  flow  hood,  the  substrates  were  exposed  to  UV  light  for 

2. 5-3.0  minutes  at  a  distance  of  2’  from  the  200  watt  high  pressure 

o 

mercury  UV  light  source.  Following  this,  a  thin  500  A)  layer  of 
aluminum  was  evaporated  onto  the  photoresist.  Finally  a  layer  of  1.0  ym 
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thick  Shipley  AZ  1370  photoresist  was  spun  onto  the  photoresist/aluminum 
sandwich  and  baked  at  x70°C  for  20±5  minutes.  The  resulting 
photoresist/aluminum/photoresist  sandwiches  were  stored  In  tightly 
sealed  sample  storage  boxes  until  needed  for  the  photoresist  stencils. 

The  photoresist  coated  substrates  were  now  ready  for  the  exposure 
and  development  sequence,  summarized  in  Fig.  3.3.  The  substrates  were 
placed  in  a  vacuum  contact  printing  holder  consisting  of  a  piece  of  0.5" 
thick  plexiglass  with  a  slot  milled  to  a  depth  of  0.001"  beyond  the 
sample  thickness.  Holes  were  machined  in  the  plexiglass  to  connect  this 
slot  to  a  vacuum  port.  A  flexible  chrome  mask  was  placed  chrome  side 
down  on  the  substrate,  a  vacuum  valve  was  opened  and  the  flexible  mask 

CO 

was  then  pressed  down  against  the  photoresist  by  atmospheric  pressure  . 
The  sample/holder  combination  was  then  placed  V  away  from  the  UV  light 
source  and  exposed  for  -v  90  seconds.  The  top  layer  of  photoresist  was 
then  developed  for  30-60  seconds  in  high  resolution  photoresist 
developer  (  5  parts  water  to  1  part  Shipley  AZ  351  developer),  and 
rinsed  with  distilled  water.  Next  the  aluminum  layer  was  etched. 

Finally  the  bottom  layer  of  photoresist  was  developed  for  a  few  seconds 
beyond  the  point  at  which  it  was  clear  that  all  the  photoresist  had  been 
developed  away  from  the  substrate.  (The  bottom  layer  of  photoresist 
could  be  seen  where  it  had  been  exposed  to  the  developer  by  the  etching 
of  the  A1  layer.)  This  procedure  resulted  in  a  v5  ym  undercut  all  around 
the  edges  of  the  photoresist  and  produced  the  overhanging  photoresist 
bridge  used  in  making  the  junctions.  SEM  micrographs  of  undercut  resist 
profiles  are  shown  in  Fig.  3.4.  Fig.  3.4(a)  displays  a  photoresist 
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Photoresist  Exposure/ Development  Sequence 


.PRUO/im) 
-A?  ( 500&) 


PR  (1.5^m)- Already  exposed  to 
UV  light 

Single  Crystal  Quartz  (O.OIO" thick) 


\  5  J  UV  'W 

-Mask  (Glass  +  Cr) 


Expose  Top  PR  Layer 


Develop  Top  PR  Layer 


Etch  A?  Layer 
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PR  Bridge 


Develop  Bottom  PR  Layer 
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FIG.  3.3  Photoresist  exposure/development  procedure  beginning  wi 
photoresist/  A1 /photoresist  tri layer  and  ending  with  the  suspended 
photoresist  bridge  which  is  used  to  make  the  junction. 
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bridge  fabricated  using  the  commercial  mask  described  in  Section  3.2.1, 
and  schematically  displayed  In  Fig.  3.1(c).  A  photoresist  bridge  used 
in  fabricating  tunnel  junction  at  the  center  of  a  dipole  antenna  is 
shown  in  Fig.  3.4(b).  This  micrograph  was  taken  with  the  SEM  stage  at 
an  -w80°  tilt,  viewed  along  the  antenna  axis.  Going  off  to  the  left  are 
the  parts  of  the  stencil  used  to  make  the  dc  leads.  Note  the  vertical 
sidewalls  indicating  good  contact  was  made  between  the  mask  and  the 
photoresist. 

3.2.4  Junction  Evaporation/Oxidation  Sequence 


There  are  three  ways  to  make  tunnel  junctions  using  the  overhanging 
photoresist  bridge  shown  in  Fig.  3.4.  The  first  technique  was 
developed  at  Bell  Labs  by  Dolan44  and  Dunkleberger45.  We  call  their 
original  technique  the  “overlap"  technique.  This  sequence  begins  with 
the  deposition  of  the  base  electrode  material  (Sn)  at  an  angle  to  the 
substrate  normal  (  +55  °) ,  followed  by  a  glow  discharge  of  the  base 
electrode,  and  ends  with  the  evaporation  of  the  counter  electrode  (Pb) 
at  an  angle  of  -50°.  For  this  evaporation  sequence,  substrates  were 
placed  on  a  rotatable  sample  holder  which  was  cooled  to  liquid  nitrogen 
temperature  in  a  cryopumped  evaporator  at  a  base  pressure  of  2xl0“7 
torr.  The  base  electrode  (  always  Sn  )  was  oxidized  by  a  dc  glow 
discharge  in  99.9999%  pure  oxygen  at  a  pressure  of  30  mtorr  for  30-60 
seconds  with  12  mA  of  current  at  1.1  kV.  Fig.  3.5(a)  shows 
schematically  the  evaporation  sequence,  while  Fig.  3.5(b)  displays  an 
SEM  micrograph  of  a  sample  made  with  this  procedure.  Typically  this 


FIG.  3.4  (a)  SEM  micrograph  of  a  photoresist  bridge  suspended  over  a 

polished  sapphire  substrate,  made  with  the  mask  in  Figure  3.1(c), 
following  the  procedure  of  Fig.  3.3  .  The  micrograph  was  taken  with  a 
•u80°  tilt  on  the  SEM  stage,  (b)  SEM  micrograph  of  a  suspended 
photoresist  bridge  used  in  fabricating  a  small-area  junction  at  the 
center  of  a  dipole  antenna.  The  micrograph  was  taken  with  an  ->,80  tilt 
on  the  SEM  stage;  the  view  is  from  along  the  antenna  axis.  DC  leads 
will  be  on  the  left  in  the  completed  sample. 
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technique  produces  junctions  with  areas  from  2-10  vm. 

There  are  two  variants  on  this  method  ,  developed  by  Howard  et 
al.  which  can  produce  junctions  with  somewhat  smaller  areas.  These 
two  variants  are  called  the  "resist-aligned"  scheme,  shown  schematically 
in  Fig.  3.6(a),  and  the  "edge-aligned"  scheme,  shown  in  Fig.  3.7(a). 

In  the  resist-aligned  sequence,  the  area  of  the  junction  is  reduced  from 
that  of  the  overlap  junction  by  the  deposition  of  an  insulating  layer  of 
Ge  after  the  base  electrode  deposition  at  an  angle  (%+42°  for  example) 
which  is  slightly  less  than  the  angle  of  the  base  electrode  deposition. 

In  Fig.  3.6(b)  we  see  an  SEM  micrograph  of  a  junction  produced  this 
way.  The  final  variant  of  the  overlap  technique,  the  edge-aligned 
technique,  begins  with  the  evaporation  of  Sn  at  normal  incidence.  The  Ge 
lav'er  is  evaporated  at  an  angle  such  that  it  breaks  over  the  Sn  edge; 
the  junction  is  glow  discharge  oxidized  as  before,  and  the  Pb  counter 
electrode  is  evaporated  at  a  shallow  angle  to  insure  it  touches  the 
Sn-SnO  surface.  In  Figure  3.7(b)  we  display  an  SEM  micrograph  of  a 
sample  produced  using  this  scheme.  We  clearly  see  that  the  Ge  broke  over 
the  Sn  edge  and  was  deposited  on  the  substrate  surrace. 

The  smallest  area  junctions  produced  by  optical  lithography  are 

those  made  using  the  edge-aligned  technique.  These  could  have  areas  as 
2  _Q2 

small  as  ->.0.2  urn  or  %2xl0  cm  with  the  y-Z  urn  minimum  linewidth 

lithography  we  use.  Unfortunately,  we  were  not  able  to  make  many  good 

samples  using  this  technique.  In  fact,  almost  all  of  the  junctions 

tested  in  the  FIR  were  fabricated  using  the  resist-aligned  scheme. 

2 

Areas  ranged  from  0.5-2. 5  ym  with  resistances  from  16-380  n  for  these 


FIG  3.6  (a)  Evaporation/oxidation  procedure  for  resist-aligned 

junctions,  (b)  SEM  micrograph  of  a  sample  made  with  this  procedure. 

The  bright  metal  is  Pb,  while  the  material  between  the  two  Pb  tips  is 
Ge.  Under  the  Pb  in  the  lower  part  of  the  photo  is  the  tunnel  junction 
On  the  left  are  the  dc  leads. 
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FIG.  3.7  (a)  Evaporation/oxidation  procedure  for  edge-aligned 

junctions,  (b)  SEM  micrograph  of  junction  made  with  this  procedure.  On 
the  left  is  the  Pb  layer  contacting  the  SnO  on  the  edge  of  the  Sn  layer. 
The  Ge  layer  is  on  the  substrate,  underneath  the  Pb,  but  close  to  the  Sn 
layer,  and  the  Ge  is  also  on  the  surface  of  the  Sn. 
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resist-aligned  junctions 


Unfortunately,  good  junctions  alone  do  not  couple  significantly  to 
submillimeter  laser  radiation.  In  the  next  section  we  will  describe  the 
design  of  the  dipole  antenna  used  to  couple  the  Gaussian  free  space 
output  mode  of  the  FIR  laser  to  the  junction. 
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3.3  Antenna  Coupling  to  the  FIR  Laser  Source 

3.3.1  Introduction 

Thin  film  tunnel  junctions,  unlike  point  contacts,  can  be  reliably 
made  and  can  have  excellent  thermal  cycling  characteristics  and 
stability.  Unfortunately,  the  dielectric  substrate,  essential  to  planar 
microfabrication  technology,  makes  coupling  a  FIR  laser  gaussian  output 
mode  to  the  junction  much  more  difficult  than  coupling  it  to  a  point 
contact.  A  point  contact  makes  a  good  long  wire  antenna  by  virtue  of 
the  fact  that  it  is  a  long  (^5  mm  ),  thin  (x75  pm  )  wire  sharpened  at  the 
tip  in  contact  with  a  metal  ground  plane.  This  natural  antenna  produces 
an  antenna  pattern,  cylindrically  symmetrical  about  the  wire  axis,  with 
its  maximum  at  a  polar  angle  close  to  that  axis.53,54  The  output  mode 
of  a  FIR  laser  can  couple  quite  efficiently  to  such  an  antenna 
pattern.55 

The  simplest  antenna  structure  one  could  fabricate  on  a  dielectric 
substrate  is  a  dipole  antenna.  In  fact,  planar  dipole  antennas  were 
fabricated  for  the  337  pm  HCN  FIR  laser  line  by  Javan's  group  in  1974. 56 
Using  a  Ni-oxide-Cr  thin  film  diode  at  the  center  of  a  dipole  antenna, 
they  were  able  to  detect  the  HCN  laser  radiation  using  the  non-linearity 
of  the  diode  I-V  curve.  Gradually  shortening  the  length  of  the  antenna 
by  trimming  it  with  a  diamond  knife,  they  found  its  resonant  length  to 
be  (”v.l/3)x  instead  of  (\1/2)X  as  would  be  expected;  the  difference  is 
due  to  the  presence  of  the  dielectric  substrate.55 


Since  the  pioneering  work  of  Javan's  group  there  have  been  several 

other  attempts  at  coupling  FIR  radiation  efficiently  to  antennas  on 
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dielectric  substrates.  Besides  dipole  antennas,  other  groups  have 

constructed  slot  antennas58,  V  antennas59,  bow-tie  antennas60,  long-wire 
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antennas  ,  and  dielectric  antennas  .  For  the  purposes  of  this  work  it 
is  necessary  to  discuss  only  dipole  antennas,  since  that  was  the  only 
configuration  used.  In  the  final  chapter  we  will  describe  the  possible 
advantages  of  bow  tie  antennas  for  certain  kinds  of  experiments.  The 
discussion  of  antennas  which  follows  here  is  separated  into  three  parts. 
The  first  reviews  microwave  modelling  experiments  on  dipole  antennas  on 
thin  substrates.  In  the  second  section  the  radiated  power  pattern  of  a 
current  element  on  a  thick  substrate  is  derived  and  theoretical  results 
on  the  power  pattern  of  a  resonant  dipole  on  thick  substrates  of 
different  dielectric  constants  are  presented.  The  final  part  of  this 
section  presents  our  design  for  the  dipole  antenna  and  the  el ectrof ormed 
copper  horn  which  focusses  the  laser  radiation  onto  the  antenna. 


3.3.2  Dipole  Antennas  on  Thin  Substrates 

For  a  dielectric  substrate  of  thickness  t,  relative  dielectric 
constant  cr,  and  electromagnetic  radiation  of  wavelength  A,  the 
dimensionless  effective  thickness  of  the  substrate  te^  is  defined  by 
the  relation: 


We  use  MKS  units,  with  the  dielectric  constant  written  as 


e0  Is  the  permittivity  of  free  space,  and  er  is  the  relative 
dielectric  constant.  Mizuno,  Daiku,  and  Ono  performed  a  microwave 
modelling  experiment  at  10  GHz  for  the  case  where  tgff  <  1.  Their  main 
results  are  summarized  in  Fig.  3.8.  If  a  is  defined  to  be  the  resonant 
length  of  a  dipole  antenna  on  the  dielectric  substrate  divided  by  its 
resonant  length  in  free  space,  then  o  decreases  for  increasing  effective 
thickness,  tg^,  for  a  given  value  of  er  up  to  t^  =  0.15.  This 
result  is  shown  in  Fig.  3.8(a)  for  several  values  of  e  up  to  an 
effective  thickness  of  tg^  *0.6.  The  reduction  factor  a  was  found  to 
be  constant  for  t  ^  >  0.15,  for  a  given  er,  and  these  values  of  a  are 
plotted  against  in  Fig.  3.8(b).  The  antenna  power  pattern  was  also 
measured.  Fig.  3.8(c),  for  the  E-plane  and  H-plane.  (The  E-plane  is 
the  plane  defined  by  the  electric  field  vector,  and  the  vector  from  the 
center  of  the  dipole  to  the  maximum  of  the  power  pattern.  The  H-plane 
is  defined  analogously.®4)  Note  that  the  E-plane  power  pattern  differs 
only  somewhat  from  that  of  a  dipole  in  free  space,  having  only  x30%  more 
power  directed  into  the  substrate,  at  the  maximum  perpendicular  to  the 
substrate,  than  in  the  free  space  case.  The  H-plane  pattern,  however, 
departs  strongly  from  the  cylindrical  symmetry  of  the  free-space  dipole 
pattern,  having  considerable  power  radiated  into  the  plane  of  the 
substrate.  Here,  the  waves  are  trapped  as  surface  waves  by  total 
internal  reflection  at  the  substrate  interfaces. 

3.3.3  Power  Patterns  of  Current  Elements  and  Dipol es 

65 

Rutledge  and  Muha  have  calculated  the  power  pattern  of  a  current 
element  on  a  dielectric  half-space,  using  a  physical  argument  based  on 


FIG.  3.3  (a)  Reduction  factor,  a,  as  a  function  of  the  effective 

thickness,  t  for  different  values  of  e  .  (b)  Values  of  a  for 

t  fr>0.15  for  {he  values  of  e  in  (a),  (c)  E-plane  and  H-plane  antenna 

patterns  for  a  dipole  antenna'on  a  thin  substrate.  [Adapted  from  Ref. 
63.] 
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the  Lorentz  reciprocity  theorem.66  They  showed  that  the  current  element 
radiates  much  more  strongly  into  the  dielectric  than  into  air  (or 
vacuum).  The  discussion  which  follows  is  closely  based  on  the  analysis 
of  Rutledge  and  Muha  (Ref.  65),  and  analyses  presented  in  Refs.  66-67. 
The  reader  is  urged  to  study  these  articles  for  a  more  complete  analysis 
than  is  presented  here.  Readers  who  prefer  to  forgo  the  discussion 
presented  here  are  urged  to  skip  to  the  answers,  Eqs.  (3.3)-(3.7), 
which  are  plotted  in  Fig.  3.9(b). 

The  Lorentz  reciprocity  theorem,  as  stated  by  Ball  ant ine,  applies 

to  the  situation  of  two  antennas,  and  A^,  located  at  positions  0^  and 

02,  respectively,  and  having  arbitrary  orientations.  Consider  first  A1 

to  be  a  transmitting  antenna  and  a  receiving  antenna,  and  assume  that 

Aj  transmits  with  some  average  power  P.  If  A^  is  considered  to  be  a 

transmitting  antenna  and  A^  a  receiving  antenna,  and  A^  transmits  with 

the  same  average  power  P,  as  A^  did  previously,  then  the  phase  and 

intensity  of  the  electric  field  at  the  receiver  A^  is  the  same  as  that 

produced  before  at  A^,  when  A^  was  the  receiver  and  Aj  was  the 

transmitter.  This  result  is  independent  of  the  electrical  properties  of 

the  intervening  media,  i.e.  conductivity,  inhomogeneity,  stratification, 

and  so  on.  It  is  also  independent  of  the  specific  forms  of  the 
66 

antennas. 

This  theorem  is  particularly  useful  in  cases  where  it  is  desired  to 
calculate  the  far-field  power  pattern  of  a  transmitting  antenna,  in  a 


situation  where  a  direct  calculation  is  difficult  due  to  the  presence  of 
ground  planes  or  dielectric  discontinuities  located  near  the  antenna. 

The  antenna  could  then  be  considered  as  a  receiving  antenna,  and  the 
electric  field  at  the  antenna  produced  by  the  far-field  (essentially 
plane  wave)  power  pattern  of  another  antenna  could  be  calculated.  From 
this  electric  field  the  power  pattern  can  be  calculated. 

This  theorem  can  be  applied  to  the  situation  of  a  current  element 
at  an  air  (or  vacuum) /dielectric  interface,  shown  in  Fig.  3.9(a).  We 
wish  to  calculate  the  far-field  pattern  of  the  current  element  situated 
at  the  dielectric  interface,  1^.  We  calculate  the  pattern  for  the  plane 
perpendicular  to  the  current  element  axis,  i.e.  the  H-plane.  The 
current  element,  Ij ,  produces  an  electric  field  Ej,  in  the  far-field  at 
an  angle  0  .  In  order  to  calculate  the  power  pattern,  we  consider  the 
current  element  ^  as  a  receiver,  and  compute  the  electric  field  E2  at 
Ij  from  the  current  element  I2  located  at  Ej.The  reciprocity  theorem 
stated  above  says  that  if  the  two  currents  Ij  and  I2  are  equal  then  the 
electric  fields  Ej  and  E2  are  equal.  The  variation  of  E2  with  angle  0fl 
gives  the  power  pattern  (for  the  H-plane)  in  the  air.  If  I2  were  placed 
deep  inside  the  dielectric,  the  reciprocity  theorem  could  be  used  in  a 
like  manner  to  compute  the  H-plane  pattern  in  the  dielectric. 

In  order  to  get  an  intuitive  understanding  of  the  problem,  consider 
a  plane  wave  incident  on  the  dielectric  from  the  air  at  normal 
Incidence,  0a*0.  0^  is  the  angle  of  refraction.  If  n  is  the  index  of 

refraction  of  the  dielectric,  then  the  ratio  of  the  transmitted  electric 
field  to  that  of  the  incident  field  is  2/(n+l).  For  large  n  this  ratio 


Air 


FIG.  3.9  (a)  Current  element  at  an  air-dielectric  interface,  (b) 

Emplane  and  H-plane  power  patterns  from  this  current  element.  [Adapted 
from  Ref.  65.1 
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is  x2/n  and  is  small.  If  the  plane  wave  were  incident  from  the 
dielectric  side  of  the  interface,  the  ratio  of  the  transmitted  electric 
field  to  the  incident  electric  field  is  large,  %2n/(n+l),  or  for  large  n 
the  ratio  is  \2.  The  average  power  density  of  the  waves  is  given  by  the 

4  4  1  O 

Poynting  vector  <S>= ( 1/2 )  Re(ExH),  which  becomes  S=(1/2)Z  |E|  ,  where 

^1/0 

Z=(u/e)  '  ,  the  wave  impedance  of  the  medium.  For  most  dielectrics, 

v=vi0,  and  e=ere0,  so  that  Z=(£r)-1/2Z0=Z0/n  using  n=(er)+1/2,  and 
Z0=(ii0/^o^+1//2ss277  n  the  wave  1mPeclance  free  space. 

The  transmitted  electric  fields,  by  the  reciprocity  theorem,  are 

equal  to  the  radiated  far-fields  produced  by  the  current  element  at  the 

interface  when  it  is  considered  as  a  transmitting  antennna.  Therefore, 

the  power  density  in  the  far-field  on  the  air  side  is  Sa=(l/2ZQ) 

|E.j|  (2/n)  where  Ei  is  the  incident  electric  field,  whereas  on  the 

7  7 

dielectric  side  the  power  density  is  Sd=(n/2Z0) | Ei |2  .  The  ratiu  of  the 
far-field  power  density  in  the  dielectric  to  that  of  the  power  density 
in  the  air  is  y-n  /!•  The  antenna  at  the  interface  radiates  much 
more  power  into  the  dielectric  than  the  air. 

This  argument  can  be  generalized  to  angles  other  than  the  substrate 
normal.  Rutledge  and  Muha  have  done  this;  we  will  briefly  sketch  their 
calculation  for  the  power  pattern  on  the  air  side.  The  electric  field 
E2  at  the  substrate  interface  at  Ij  is  calculated  assuming  is  the 
radiator  (transmitter).  The  incident  electric  field  is 

E.  =  J  ““oh  e-jkor 
1  4.r 

where  h  is  the  length  of  the  current  element.  The  field  E2  is  the 


(3.1) 
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incident  electric  field  times  the  transmission  coefficient  t  ,  for 
thfc  electric  field  perpendicular  to  the  plane  of  incidence 


2  cos  ea 

fl 

cos  0,  +  n  cos  e  . 
a  d 


(3.2) 


Note  for  9a=0,  9d=0,  the  transmission  coefficient  becomes  T=2/(n+l), 
the  transmission  coefficient  for  normal  incidence  discussed  above.  As 
in  the  previous  discussion,  reciprocity  says  that  the  electric  field  E ^ 
at  Ij  due  to  Ip  is  equal  to  the  field  E^  at  Ig  is  I^I*,.  The  power 
pattern  of  ^  is  Sa(©a)=(l/2Z0) | E1 |2=(1/2ZQ) |Ei t | 2.  If  we  define 


0=(2/Zo)[u,uoIh/41rr]':  ;  then 

cos  e  2 

Sa^ea^  “  6  ^cos  e  +  n  cos  e  ^ 
a  a 


(3.3) 


is  the  H-plane  pattern  in  the  air.  For  the  H-plane  pattern  in  the 
dielectric,  9d  becomes  the  angle  of  incidence,  9a  is  the  angle 
refraction,  so  that  the  transmission  coefficient  r  is  now  given  by: 


2n  cos  0  . 

_  / _ d  \ 

T  "  vcos  0,  +  n  cos  0/ 
a  d 


(3.4) 


The  power  density  is  Sd(9d)=(n/2Z0) lE^r p  where  t  is  given  by  Eq.  (3.4) 
and  Ei  as  in  Eq.  (3.1).  The  power  pattern  is 


n  cos  @d  t 

Sd^ed^  "  6n  ^ |cos  ©a  +  n  cos  @d| ^ 


(3.5) 


in  the  dielectric. 


The  E-plane  patterns  can  be  calculated  in  a  similar  manner,  and  we 

simply  state  the  results.  Rutledge  and  Muha  found 

cos  0  cos  0  .  7 

Sa^ea^  "  6  ^n  cos  e  +  cos  q) 

a  a 


(3.6) 


for  the  E-plane  pattern  in  the  air,  and 


Sd^0d^  =  Bn 


n  cos  ea  cos  e  . 

I - § - 1_| 

'n  cos  e,  +  cos  e ■ 
a  d 


(3.7) 


for  the  E-plane  pattern  in  the  dielectric.  Equations  (3.3),  (3.5), 
(3.6),  and  (3.7)  are  plotted  in  Fig.  3.9(b)  for  n=2,  the  index  of  fused 


quartz.  Note  that  no  power  is  radiated  along  the  plane  of  the  substrate 
in  this  case  since  it  was  assumed  the  substrate  was  infinitely  thick. 


There  is  a  null  in  the  power  pattern  at  the  interface,  implying  a  weak 
coupling  between  antennas  if  they  were  made  into  arrays. 


The  power  patterns  for  resonant  dipoles  on  thick  substrates  have 

been  calculated  and  measured  by  Brewitt-Taylor  and  Rees®7,  and  are  only 

slightly  different  from  the  patterns  for  the  current  element.  Figs. 

3.10(a)  and  3.10(b)  display  their  calculated  patterns  for  the  H-plane 

and  E-plane,  respectively  for  two  different  dielectric  constants.  The 

major  difference  between  the  calculated  patterns  for  the  dipole  as 

compared  to  the  current  element  is  that  in  the  E-plane  pattern,  the 

dipole  has  a  slightly  larger  central  lobe  than  the  current  element  (  for 

the  er  =  4  case  the  only  one  we  could  use  for  direct  comparison  from 

the  figures).  Brewitt-Taylor  and  Rees  also  note  that  the  resonant 

frequency  follows  roughly  as  n-*  from  that  of  the  free  space  dipole, 

where  n=[(er+l)/2]  ^  ,  consistent  with  the  resonant  length 
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measurements  of  Mizuno  Daiku,  and  Ono.  The  peak  admittance  is 

proportional  to  n,  results  consistent  with  the  microwave  modelling 

65 

measurements  of  Rutledge  and  Muha  on  bow-tie  antennas. 
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FIG.  3.10  (a)  H-plane  power  pattern  of  a  resonant  dipole  antenna  at  an 

air-dielectric  interface  for  different  values  of  £r.  (b)  E-plane  power 

fattern  of  the  resonant  dipole  of  (a)  for  the  same  values  of  er. 

Adapted  from  Ref.  67] 
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3.3.4  Design  Decisions 

In  this  work,  radiation  from  a  FIR  laser  source  was  focussed  by  an 
f/4  TPX  lens  and  an  electroformed  copper  cone  through  the  back-side  of 
the  sample  substrate,  containing  a  small-area  tunnel  junction  at  the 
center  of  the  substrate.  Fig.  3.11  illustrates  the  complete  coupling 
arrangement,  seen  in  cross-section,  as  viewed  in  the  transverse  optical 
access  dewar. 

The  antenna  was  fabricated  following  the  design  rules  developed  by 

63  2 

Mizuno,  Daiku.  and  Ono  .  Z-cut  crystal  quartz  substrates,  2.54  x  10 

cm  thick,  were  used.  Single  crystal  quartz  has  an  index  n  *  2.34  in  the 

100  urn  to  1.2  mm  wavelength  region  of  interest  to  this  experiment.  The 

effective  thickness  was,  tgff  *  1.2  ,  for  the  496  ym  methyl  fluoride 

(CH3F)  laser  line;  it  decreases  to  tfiff  *  0.48  for  the  1.22  mm  isotopic 

methyl  fluoride  (C  H^F)  laser  line.  The  antenna  length  was  shortened 

to  77%  of  its  free  space  resonant  length  (0.48a)  to  0.37A,  or  lre$  =  183 

ym  for  the  496  ym  line  and  lrgs  =  451  ym  for  the  1.22  mm  line.  Figure 

3.12(a)  displays  a  dipole  antenna  of  182  ym  length,  and  5  ym  linewidth. 

At  its  center  is  a  small-area  junction,  at  a  magnified  scale  in  Fig. 

3.12(b).  This  particular  junction  had  a  resistance  of  24fl,  an  area  of 

-8  2  3  2 

2.4  x  10  cm  ,  and  critical  current  density,  J  ,  of  1.1  x  10  A/cnr  at 

T  *  1.4  K. 


A  conical  horn  with  a  short  cylindrical  waveguide  section  was  used 
as  a  feed  for  the  dipole  antenna,  completing  the  coupling  of  the 
junction  to  the  FIR  laser  source.  The  design  of  this  horn  was  based  on 


FIG.  3.11  Detail  of  transverse  optical  access  dewar.  FIR  laser 
radiation  is  focussed  by  the  TPX  lens  onto  the  cone,  through  the  back 
side  of  the  substrate,  onto  the  junction-antenna  structure. 
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work  of  Blaney  at  the  National  Physical  Laboratory  in  Great  Britain.69 
The  horn  had  a  conical  section  35  mm  in  length,  an  opening  11.4  mm  in 
diameter,  and  an  opening  angle  of  -».180.  The  waveguide  section  was  7  mm 
long  and  0.79  mm  in  diameter.  This  cone  was  fabricated  by 
electrodeposition  of  copper  on  a  polished  stainless  steel  mandrel,  in 
order  to  keep  the  surface  roughness  to  a  minimum.  The  cone  has  a 
calculated  gain  of  29  db,  or  a  factor  of  760,  over  that  of  an  isotropic 
radiator. The  waveguide  section  had  a  cut-off  wavelength  of  1.35  mm 
for  the  TE^  mode,  the  lowest  order  or  dominant  mode.  Cut-off 
wavelengths  for  some  higher  order  modes  are  1.03  mm  for  the  TMQ1  mode, 
0.81  mm  for  the  TE2j  mode,  0.65  mm  for  the  TM^  and  TEQ1  modes,  etc.71 
Therefore  although  the  waveguide  could  only  operate  in  the  fundamental 
mode  at  the  1.22  mm  laser  line,  unfortunately  at  shorter  wavelengths, 
such  as  the  strong  496  urn  line,  several  modes  could  be  excited  at  once. 

In  the  usual  operation  of  a  horn,  an  adjustable  tuning  stub  is 
inserted  in  the  waveguide  section  about  3Xg/4  in  front  of  the  device, 
usually  a  point  contact.  ( Ag  is  the  wavelength  of  the  mode  in  the 
guide.)  A  moveable  plunger  terminates  the  waveguide  behind  the  point 
contact.  A  variation  on  this  scheme,  in  which  only  a  moveable  plunger 
was  used  to  tune  out  the  point  contact  reactance,  performed  successfully 
at  0.66  mm  at  the  NPL,  achieving  a  reasonable  coupling  efficiency  of 
%10%.69 

In  this  work  no  attempt  was  made  to  resonate  out  the  junction 
reactance  by  either  tuning  stubs  or  plungers.  The  coupling  arrangement, 
used  in  this  experiment,  in  which  the  radiation  is  focussed  through  the 


back-side  of  the  substrate,  takes  advantage  of  the  effect  of  the 
dielectric  substrate  on  the  antenna  pattern.  During  a  run  the  cone 
could  be  rotated  about  the  center  of  the  dewar  to  optimize  the  coupling 
of  the  radiation  to  the  junction. 

The  FIR  radiation  came  from  the  optically  pumped  FIR  laser  built  by 
17 

D.  Weitz.  The  characteristics  of  this  laser  have  been  thoroughly 
described  before,  but  will  be  discussed  briefly  in  the  following 


section. 


77 


3.4  Far- infrared  Laser  Source 

3.4.1  Introduction 

In  the  thirteen  years  since  the  first  observation  of  FIR  lasing 
action  in  an  optically  pumped  molecular  gas  laser  by  Chang  and 
Bridges7'*,  over  1300  laser  lines  have  been  reported.74,75  Optically 
pumped  FIR  lasers  have  become  a  reliable  source  of  radiation  in  the  40 
urn  to  2000  ym  wavelength  region.  The  laser  used  in  this  work  has  been 
documented  in  detail  in  connection  with  earlier  work  on  point  contacts, 
so  we  will  only  briefly  discuss  its  theory,  design,  operation  and 
performance. 

3.4.2  Qualitative  Theory 

The  theory  of  the  FIR  laser  can  be  understood  most  simply  in  terms 

of  the  energy  level  diagram  of  a  specific  molecule.  Since  the  symmetric 

73  76  81 

top  molecule,  methyl  fluoride,  CH^F,  has  been  thoroughly  studied  ’  ’ 

we  will  begin  by  presenting  a  diagram  of  its  energy  levels  relevant  to 
the  FIR  laser  process,  Fig.  3.13.  The  state  of  the  molecule  is 
described  by  three  quantum  numbers,  assuming  only  the  or  C-F  bond 
stretching  vibrational  mode  is  allowed:  v,  the  vibrational  quantum 
number,  J  related  to  the  total  angular  momentum,  and  K,  the  axial 
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component  of  the  angular  momentum,  where  K  =  0,  +l,+2 . iJ. 

Methyl  fluoride  molecules  in  the  v  =  0,  J  =  12,  K  =  2  state  are  excited 
by  the  9P(20),  9.55  ym  C02  pump  laser  to  the  v  *  1,  J  *  12,  K  *  2 
excited  state.  The  lasing  transition  is  between  the  J  =  12  and  0  =  11 
rotational  states.7^  This  transition  produces  the  604  GHz,  496  ym  laser 


Energy  Level  Diagram 
for  Optically- Pumped  CH3F  FIR  Laser 


K*1  2  3 


FIG.  3.13  Methyl  fluoride  energy-level  diagram.  The  CCL  pump  line 
coincides  with  the  v=0-»l,  J=  12 ,  k=2  transition.  The  496  ym  FIR  lasing 
transition  is  the  v=l ,  K=2,  J=12-»ll  transition.  aK  denotes  rotational 
quantum  number  changing  transitions  which  tend  to  destroy  the  population 
inversion  of  the  lasing  level,  as  does  aJ.  Vibrational  mode  changing 
collisions  also  destroy  the  population  inversion.  Molecules  relax  back 
into  the  ground  state  by  collisions  with  the  cavity  walls.  [Adapted 
from  Ref.  80.] 


line,  one  of  the  strongest  FIR  laser  lines  reported  in  the  literature, 
and  the  line  used  most  in  this  work. 


The  FIR  output  power  is  limited  by  fundamental  factors,  such  as 

various  relaxation  rates  of  the  molecule  out  of  its  excited  state,  and 
82 

the  Manley-Rowe  limit  on  single  photon  energy  down  conversion 

processes.  Experimentally  controllable  factors  such  as  the  fraction  of 

pump  power  absorbed  in  the  lasing  gas,  the  fraction  of  pimp  power  lost 

due  to  absorption  in  the  walls  and  mirrors  of  the  cavity,  cavity  losses 

of  the  FIR  beam,  and  transmissivity  of  the  FIR  output  coupling  mirror, 

79-81 

also  limit  the  FIR  output  efficiency.  The  expression  for  the 

theoretical  maximum  output  power  of  the  FIR  laser  in  terms  of  the  input 
pump  power,  Pp,  pump  frequency  vp,  and  output  frequency  vpIRis: 


P 


1  ,hvFIR 
FIR  2  k  hVp 


(3.8) 


and  the  efficiency  n  is  given  by  n  =  ppiR^Pp*  T,ie  factor  3-/2  comes 


from  the  fact  that  the  degeneracies  of  the  upper  and  lower  lasing 
sub-levels  are  approximately  equal  since  the  molecules  are  in  high 
rotational  states;  this  means  that  for  every  two  molecules  excited 
only  one  can  contribute  to  the  FIR  output  by  emitting  a  photon  and 
decaying  to  the  lower  state,  thereby  destroying  the  population  inversion 
necessary  for  lasing  action.  The  factor  hVpjR/  hvp  is  the  Manley-Rowe 
condition,  limiting  the  maximum  power  to  the  ratio  of  output  photon 
energy  to  the  input  photon  energy,  as  in  a  one-photon  in,  one-photon  out 


process. 


This  maximum  theoretical  laser  efficiency  is  n  ,  =9.6xlO"J  at 

J  max 

A=496  ym.  The  measured  efficiency  for  this  laser  is  only  -*.7-10  %  of  the 
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maximum,  which  is  approximately  the  same  as  the  efficiency  measured  by 
79 

other  workers. 

Experimentally,  the  FIR  cavity  losses  can  be  limited  to  a  few 
percent  by  a  careful  design  of  the  FIR  cavity,  while  allowing  a 
reasonable  fraction  of  the  FIR  power  to  be  coupled  out.  A  highly- 

overmoded  dielectric  waveguide  cavity  with  external  mirrors  has  low 

ft? 

propagation  losses  for  the  linearly  polarized  EH^  cavity  mode,  which 

04 

transforms  to  the  TEMQo  free  space  gaussian  mode.  Losses  due  to 
absorption  of  the  pump  beam  by  molecules  in  the  excited  state  can  be 
minimized  by  keeping  the  lasing  gas  to  low  pressures.  The  fraction  of 
pump  power  absorbed  is  usually  small,  due  to  the  low  pump  absorption 
coefficients  (a=  0.01  m  for  CH^F)  f  It  would  be  possible  to  increase 
the  pump  absorption  by  increasing  the  pressure,  but  this  also  increases 
the  rate  at  which  excited  atoms  are  thermal i zed  by  rotational  quantum 
number  changing  collisions,  effectively  limiting  the  pressure  to  -*-100 
mtorr.  In  practice  as  long  a  cavity  as  practicable  is  used  to  increase 
the  pump  absorption.  The  relaxation  rate,Yy,  of  molecules  in  the 
excited  vibrational  state  to  de-excite  and  return  to  the  vibrational 
ground  state  after  lasing  is  slow  compared  to  the  rotational 
cross-relaxation  of  molecules  within  the  excited  vibrational  band.  The 
lower-laser-state  lifetime  is  effectively  (vv)"^,  and  therefore 

decreasing  this  lifetime  would  increase  the  output  power.  Practical 
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solutions  are  to  decrease  the  diameter  of  the  cavity  ,  thereby  reducing 
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the  diffusion  time  for  molecules  to  reach  the  walls  and  de-excite,  or  to 
introduce  a  buffer  gas  to  collisionally  de-excite  the  molecules.85 

3.4.3  Operation  and  Performance 

Fig.  3.14  displays  the  complete  FIR  laser  system.  A  modified 
86 

commercial  CO^  laser  provided  the  pump  beam.  The  punp  laser  had  a  12 

mm  diameter,  1.8  m  long  cavity,  with  the  active  region  enclosed  by  two 

ZnSe  Brewster  angle  windows.  This  laser  had  been  modified  by  using  an 

%20  cm  longer  cavity  than  usual  for  increased  output  power,  and  by  the 

removal  of  the  output  coupling  mirror  from  the  end  of  the  laser  cavity, 

and  its  replacement  by  a  second  Brewster  angle  window.  A  75  lines/mm 

grating  was  used  to  tune  the  laser  to  a  single  line,  and  also  served  as 

one  end  mirror.  The  output  coupling  mirror  was  mounted  on  a 

piezo-electric  transducer  (PZT)  stack  for  fine  tuning  of  the  operating 

frequency,  and  was  semi-transparent  with  a  20  m  radius  of  curvature. 

The  laser  was  mounted  on  a  granite  slab  and  enclosed  by  a  metal  cover 

for  stability  against  frequency  drift  due  to  thermally  induced  cavity 
72 

length  changes. 

In  the  typical  operation  of  the  laser,  12-25  W  of  output  power 
could  be  achieved  over  the  9.5  urn  band,  and  20-50  W  in  the  10.5  urn  band 
with  a  nearly  Gaussian  output  beam. 

The  single  most  important  condition  for  a  stable,  single-line, 
Gaussian  output  beam  was  that  the  Brewster  windows  were  clean.  After 
%100  hours  of  operation,  the  laser  would  not  operate  on  a  single  line. 
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but  would  Instead  lase  on  2-3  lines  at  once  with  a  diminished  and 
unstable  output  power.  At  this  time,  the  windows  were  removed  from  the 
laser  and  ultrasonlcally  agitated  for  * 1/2  hour  in  laboratory  detergent, 
followed  by  a  15  minute  ultrasonic  agitation  in  distilled  water,  and 
blown  dry  with  pure  nitrogen  gas.  If  the  mode  shape,  and  power,  and 
stability  were  not  drastically  Improved  by  this  treatment  the  old 
windows  were  removed  and  discarded.  New  ones  were  installed  on  the 
Brewster  mounts,  cleaned  with  the  procedure  noted  above,  and  installed 
on  the  laser. 

The  FIR  cavity  consisted  of  a  38  mm  diameter,  3  m  long  Pyrex  pipe 
with  two  flat  end  mirrors.  The  COg  laser  beam  was  focussed  into  the  FIR 
cavity  by  an  f/50  coupling  mirror  through  the  BaFg  window,  and  through 
the  3  mm  diameter  hole  in  the  input  coupling  mirror  (  flat  )  mounted  at 
the  end  of  the  dielectric  waveguide.  This  mirrror  was  mounted  on  a 
translation  stage  between  2  stainless  steel  bellows  vacuum  seals.  This 
mirror  and  also  the  output  coupling  mirror  were  on  glmbal  mountings  for 
adjustment  of  the  cavity  finesse  (essentially  the  cavity  quality  factor 
Q).  The  mirror  could  be  translated  to  adjust  the  cavity  length  to  the 
proper  resonant  FIR  cavity  mode.  The  output  coupling  mirror  consisted 
of  an  A1  thin-fllm,  capacitive  (  or  complementary  )  mesh  on  a  2  mm 
thick,  1.75  inch  diameter,  optically  polished  single  crystal  quartz 

flat.  The  capacitive  mesh  output  mirrors  were  made  by  evaporating  3000 
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A  A1  through  N1  electroformed  meshes  onto  the  quartz  flats.  *  The 

mesh  sizes  were  chosen  according  to  the  FIR  laser  line  desired,  for 

example,  a  100  lines-per-inch  Ni  mesh  was  used  to  make  the  output 
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Capacitive  mah  constants  are  given  in  tens  of  g,  the  aesh  period,  and  2a, 
the  distance  between  adjacent  A1  squares. 
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the  period  of  Hi  electroforaed 


aesh  used  to  fabricate  the  A1  capacitive 


aesh. 


coupling  mirror  for  the  496  pm,  methyl  fluoride  laser  line.  For 
completeness.  Table  3.1  displays  strong  CW  FIR  laser  lines,  the  lasing 
gas  used,  COg  pump  line,  and  FIR  output  polarization  relative  to  the  COg 
polarization,  and  capacitive  mesh  constants.  The  capacitive  mesh 
mirrors  have  a  high  reflectivity  at  the  COg  pump  frequency,  proportional 
to  the  area  coverage  of  the  A1  squares  which  make  up  the  mesh.  The  pump 
radiation  not  reflected  by  the  A1  is  absorbed  In  the  quartz  flat.  The 
reflectivity  of  the  cavity  at  the  output  frequency  could  be  altered  by 
using  output  couplers  with  different  mesh  constants. 

The  FIR  radiation  Is  transformed  from  the  EH^  cavity  mode  to  the 
TEMqo  gausslan  free  space  mode  upon  transmission  through  the  output 
coupling  mirror.  The  radiation  is  reflected  by  two  mirrors  and  then  a 
reflected  beam  from  the  single  crystal  quartz  (flat)  beamsplitter  feeds 

go 

the  thermal  detector.  The  transmitted  beam  is  focussed  by  an  f/4  lens 
of  TPX  or  polyethylene,  through  room  temperature  and  liquid  heliun 
temperature  single  crystal  quartz  windows  onto  the  electro* formed  copper 
feed  horn,  as  described  above. 

The  FIR  laser  performed  best  when  the  laser  room  was  warm,  \80°  F, 
and  when  the  COg  laser  had  been  wanned  up  for  2-3  hours  before  Injecting 
Its  beam  into  the  FIR  cavity.  The  FIR  cavity  was  always  operated  in  a 
sealed  off  mode  at  the  pressure  which  gave  the  maximum  output  power. 
Under  these  conditions  the  FIR  output  power  was  very  stable,  having  long 
term  drifts  of  at  most  *10%  in  24  minutes  of  operation,  using  only  this 
passive  stabilization  scheme  Short  term  fluctuations  were  %2-3  %, 
measured  with  a  0.3  sec  time  constant.  Fig.  3.15  displays  the  output 
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power  measured  over  a  24  minute  period.  Long  term  drifts  were 
associated  with  FIR  cavity  length  changes,  while  short  term  fluctuations 
were  connected  with  CO2  cavity  length  changes  or  gross  disturbances 
such  as  opening  or  closing  the  door  of  the  laser  room. 

Finally,  proper  alignment  of  the  C02  pump  laser  Injection  optics 
was  crucial  to  the  stable  operation  of  the  FIR  laser.  A  HeNe  laser  was 
used  for  alignment.  The  alignment  procedure  began  with  the  alignment  of 
the  HeNe  laser  and  C02  laser  beams.  A  pin-hole  made  out  of  asbestos 
sheet  was  first  aligned  with  the  pump  laser,  close  to  the  end  of  the 
laser,  and  the  position  of  the  spot  from  the  laser  beam  on  another 
asbestos  sheet  about  10'  from  the  end  of  the  C02  laser  was  noted.  A 
removeable  mirror  was  then  placed  at  the  end  of  the  C02  laser,  and  the 
HeNe  laser  beam  was  aligned  with  the  center  of  the  hole  In  the  first 
asbestos  sheet  and  the  spot  on  the  second  sheet.  The  HeNe  laser  was 
then  sent  through  the  injection  optics  of  the  FIR  laser.  The  flat 
Injection  mirror  at  the  end  of  the  cavity,  and  the  focussing  mirror, 
both  mounted  on  gimbals  and  adjusted  by  micrometer  screws,  were  adjusted 
so  the  HeNe  laser  beam  went  through  the  center  of  the  Input  coupling 
hole  through  the  center  of  the  FIR  cavity,  and  onto  the  center  of  the 
output  coupling  mirror.  The  output  coupling  mirror  was  partially 
transmitting  so  the  HeNe  beam  could  be  seen  as  a  bright,  \\  cm  diameter 
circular  spot.  There  was  never  any  need  to  remove  the  coupling  mirror 
during  the  alignment.  If  the  light  transmitted  through  the  output 
coupler  appeared  diffuse,  or  shaped  like  a  crescent  or  a  line,  then  most 
likely  the  HeNe  beam  suffered  one  or  more  reflections  with  the  walls  of 
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the  FIR  cavity  before  reaching  the  mirror.  If  there  were  transmitted 
"satellites'4  in  addition  to  the  main  spot,  then  the  cavity  mirrors  were 
not  parallel  to  each  other.  The  tilt  on  the  mirrors  was  adjusted  until 
there  was  only  one  bright  spot  at  the  center  of  the  output  mirror. 

After  this  procedure  was  completed  and  the  CO2  pump  laser  was  Injected 
Into  the  FIR  cavity,  the  FIR  cavity  length  was  adjusted  until  there  was 
a  little  FIR  output  power.  The  output  power  could  then  be  maximized  by 
adjusting  the  pressure.  Input  coupling  optics,  PZT  bias  voltage,  and 
punp  power. 

Before  I-V  curves  of  the  irradiated  junction  were  taken,  the  FIR 
laser  was  adjusted  for  maximim  power  and  stability.  I-V  curves  could 
then  be  taken  at  different  laser  powers  by  Inserting  mylar  or  plexiglass 
absorbers  at  the  output  of  the  FIR  laser. 


CHAPTER  IV 


DATA  ANALYSIS 

4.1  Introduction 

The  junctions  described  In  Chapter  III  provide  an  Ideal  test  of  the 
frequency-dependent  theory  of  superconductivity,  because  of  their  small 
area  and  nearly  ideal  dc  I-V  characteristics,  and  the  relative 
reproducibility  of  these  characteristics  from  junction  to  junction.  We 
begin  this  Chapter  with  analysis  of  the  dc  characteristics  of  low 
resistance  junctions  which  were  the  first  small-area  junctions 
successfully  fabricated  in  our  lab.  Their  excellent  dc  properties  and 
microwave  response  to  4.3  mV  with  10  GHz  radiation  demonstrated  that 
these  junctions  could  be  expected  to  work  in  the  far  infrared.  The 
results  of  the  FIR  measurements  are  presented  in  the  following  sections, 
beginning  with  the  analysis  of  the  ac  Josephson  effect  data  In  Section 
4.3.  The  photon-assisted  tunneling  data  and  analysis  is  presented  in 
Section  4.4.  Section  4.5  ends  the  Chapter  with  an  analysis  of  the 
effects  of  noise  on  the  dc  and  FIR  response  of  the  junctions,  and 
attempts  to  treat  the  whole  I-V  curve,  including  the  ac  Josephson 
effect,  the  photon-assisted  tunneling  effect,  and  noise. 

The  theoretical  discussions  of  Chapter  II  provide  the  framework 
within  which  we  will  analyze  the  data.  The  dc  characteristics  analyzed!,, 
the  Section  below  had  resistances  of  12  n  or  less,  and  had  high  IcRn 
products.  Higher  resistance  samples  (16-1600  a)  were  used  in  the  FIR 


measurements 


4.2  DC  Characteristics 

Figure  4.1  displays  the  dc  I-V  characteristics  of  a  low  resistance 
(R ,>2.4  Q) ,  high  current-density  (J  xlxlO4  A/cm^),  tunnel  junction  at  four 
selected  temperatures.  At  4.2  K  the  I-V  curve  is  non-linear  due  to  the 


presence  of  the  superconducting  Pb  electrode  and  normal  Sn  electrode, 
l.e.  the  junctions  is  a  superconductor-insulator-normal  metal  (SIN) 
junction  at  4.2  K.  Below  3.8  K,  i.e. ,  the  Tc  of  Sn,  the  junction 
develops  a  critical  current,  labeled  I£  in  the  Figure.  Additional 
structure  appears  from  the  gap  of  Sn, now  growing,  at  Apb+ASn  and  at 
APb"ASn*  as  labe^ed  *n  P’S*  4.1.  The  structure  at  the  difference  of 
the  gaps  is  due  to  thermally  excited  quasi  particles  at  the  gaps  of  Pb 
and  Sn,  which  provide  an  Increase  in  the  conductance  of  the  junction,  at 
a  voltage  In  which  the  thermally  excited  quasi  particles  of  the  Sn  can 
tunnel  Into  empty  states  in  the  Pb,  and  vice  versa.  The  two  I-V  curves, 
one  at  T»3.6  K  the  other  at  T-3.3  K,  serve  to  illustrate  the 
temperature  dependence  of  the  gap  structure.  At  the  lowest  temperature, 
T-1.4  K,  the  structure  from  the  difference  of  the  gaps  disappears  due  to 


the  freezing  out  of  the  thermally  excited  quasiparticles  as  the 
temperature  approaches  zero,  T-*0.  We  now  only  observe  the  structure 
from  the  sun  of  the  gaps  Apb+ASn,  which  arises  from  "pair  breaking,"  so 
does  not  require  any  thermal  quasi  particles. 


Experimentally ,  the  gap  of  Sn  as  a  function  of  temperature,  ASn(T), 
as  well  as  the  critical  current,  IC(T),  can  be  measured  from  the  I-V 
curves.  We  have  measured  these  structures  on  the  I-V  curves  of  low 


resistance  junctions,  whose  parameters  are  listed  in  Table  4.1.  The 
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eight  samples  shown  in  the  table  had  resistances  Rn>  ranging  from  1.5  n 

2 

to  12.5  ft,  areas.  A>  from  1.5-6. 9  (jim)  ,  and  current  densities,  J_,  from 

V 

2xl03-  2.6x10*  A/cm2.  Table  4.1  also  displays  IcRn  products  at  T=1.4  K 
which  range  from  0.40  to  1.28  mV  ,  the  lowest  (0.40  mV)  was  for  the 
highest  resistance  sample,  while  for  samples  below  4  ft,  the  IcRn 
products  were  high,  >1.0  mV,  and  did  not  correlate  with  sample 
resistance.  Finally  the  table  also  displays  estimates  of  the  junction 
capacitances  C,  in  pF,  as  well  as  the  Stewart-McCumber  parameter, 

6  *(2el .R„C)/f»,  and  the  Josephson  plasma  frequency 

v  •  1 

u.j*(2eIc/f.C)1/2=(ec)1/2/RnC. 

In  Table  4.1  the  capacitance  of  each  junction  is  estimated  from  the 

formula  Ocr(e0A/t)  where  A  is  the  junction  area,  t  is  the  oxide 

thickness,  ep  is  the  dielectric  constant,  and  eQ  is  the  dielectric 

permittivity  of  vacuum.  The  area.  A,  is  the  geometrical  area  of  the 

junction,  estimated  from  SEM  micrographs,  and  does  not  include  estimates 

of  parasitic  capacitances  such  as  the  extra  overlap  of  the  Pb  layer  on 

the  Ge  insulating  layer  visible  in  the  SEM  micrograph  of  Fig.  3.12(b). 

The  relative  dielectric  constant,  e  =6  for  SnO,  is  estimated  from  the 

ratio  of  the  wave  velocity  in  the  oxide  to  that  in  free  space,  c/c,  for 
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Sn-SnO-Sn  junctions,  as  measured  by  other  authors.  ’  For  the  oxide 

o  18 

thickness,  t,  we  use  t=20  A  as  a  nominal  value. 

The  most  poorly  known  junction  parameter  in  Table  4.1  is  the 
junction  capacitance,  and  it  is  hard  to  put  accurate  bounds  on  its 
value.  The  area  is  known  to  -v.20%,  limited  by  the  accuracy  of  the 
calibration  of  the  length  scale  on  the  SEM  (<10%),  and  the  accuracy  of 


measuring  the  area  from  the  SEM  micrographs,  probably  another  10%.  The 

relative  dielectric  constant,  er,  and  the  dielectric  oxide  thickness,  t, 

are  the  most  difficult  parameters  to  quantify  because  the  measurements 

of  c/c  upon  which  these  values  are  based  are  for  large-area,  low 

current-density  junctions  of  Sn-SnO-Sn,  and  the  measurements  were 

accurate  only  to  t300  GHz,  a  factor  of  two  lower  in  frequency  than  our 
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usual  operating  frequency  of  604  GHz.  ’  Our  junctions  are  made  of 
Sn-SnO-Pb.  The  dielectric  layer  may  contain  PbO  and  other  Pb  oxides  in 
unknown  quantities.  There  may  be  also  other  at  present  unknown 
parasitic  effects.  We  estimate  that  the  ratio  of  er/t,  is  known  only  to 
%20%,  giving  a  total  potential  error  of  v40%  on  the  values  of  the 
capacitances. 

The  Stewart-McCumber  parameter,  bc,  and  the  plasma  frequency,  ^j, 
are  quantities  derived  from  the  critical  current,  capacitance,  and 

-1/2 

resistance.  Bc-tC  and  is  therefore  known  only  to  v40%,  whereas  , 

and  is  known  to  -v.20%.  The  measured  critical  current,  I  ,  and  normal 

c* 

resistance,  Rn,  are  known  to  ’vlX,  and  limited  by  the  accuracy  of  the 
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gains  of  the  amplifiers  used  (PAR  113).  Bc  and  Uj  are  shown  with 
estimated  values  based  on  nominal  lcRp  product  of  1.2±0.1  mV,  taken  from 
the  mean  of  the  IcRn  products  of  the  five  samples,  P37-P59.  This  value 
of  the  IcRn  product,  1.2  mV,  is  close  to  the  theoretical  value,  1.30  mV, 
computed  from  Eq.  (2.2b),  assuming  Apb(0)s1.31  meV,  and  ASn(0)=0.57 
meV.  The  measured  values  of  wj  and  ec  correspond  to  the  measured  1’s  of 
the  junctions  at  >1.4  K.  The  estimated  critical  current  density, 
0c(est.),  is  based  on  the  estimated  IcRn  product,  while  the  measured 


critical  current  density,  Jc(meas.)»  was  calculated  from  the  actual 
measured  Ic  which  may  be  low  relative  to  the  nominal  value  due  to 
premature  switching  triggered  by  noise. 

Keeping  these  parameters  in  mind.  Fig.  4.2  displays  the  squares  of 
the  measured  values  of  the  gap  of  Sn  plotted  against  the  liquid  He  bath 
temperature,  determined  from  the  bath  pressure  (regulated)  using  NBS 
pressure-temperature  charts.  The  square  of  the  Sn  gap,  (2ASn)  ,  should 
be  proportional  to  Tc-T,  which  can  be  seen  from  the  BCS  expression  for  the 
superconducting  energy  gap  near  T  ,  Eq.  (4.1), 

A(T)  -  1 . 74a(0)  (l-t)*  for  T  a  T  (4.1) 

c 

where  t=T/Tc.  This  equation  says  [2ASn(T)]2=  [2(1.74)ASn(0)]2  (1-T/TJ 

so  that  the  slope  of  the  [2A$n(T)]2  versus  T  Is  -[(3.48ASn(0))2/Tc],  and 

ASn(0)  can  be  determined  from  the  slope,  given  T  ,  where  T£  is  determined 

from  the  temperature  at  which  the  gap  extrapolates  to  zero.  The  data 

from  three  of  the  four  samples  lie  on  the  same  line,  and  give 

Tc(Sn)*3.79±0.1  K,  and  ASn(0)»0.56±0.03  meV.  One  sample  shown  has 

Tc*3.81±.01  K  and  ASn(0)*0.47±.03  meV.  As  is  clear  from  Fig.  4.2,  most 

of  the  samples  have  Tc’s  and  ASn(0)*s  which  are  very  close  to  (and 
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consistent  with)  with  nominal  values  used  by  other  workers.  ’  ’ 

The  critical  current  divided  by  its  zero  temperature  value,  Ic(0) , 
when  plotted  against  t*T/Tc,  in  the  absence  of  fluctuations  should  have 
the  same  functional  form  for  all  junctions  having  electrodes  of  the  same 
materials.  This  conclusion  can  be  drawn  examining  Eq.  (2.22a)  for 
I  ( T)  for  a  junction  made  of  two  electrodes  with  equal  gaps,  and 
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FIG.  4.2  The  square  of  2Ar  ,  measured  from  data  such  as  that  presented 
in  Fig.  4.1,  versus  the  temperature  T. 


dividing  It  by  I£(0).  A  similar  result  holds  if  the  gaps  are  different. 
Fig.  4.3(a)  displays  IC(T)/IC(1.4K)  versus  t*T/Tc  for  the  four 
junctions  of  Fig.  4.2.  It  is  very  clear  that  the  data  from  all  the 
samples  lie  on  the  same  line,  i.e.  the  solid  line  in  the  Figure.  This 
curve  is  the  theoretical  curve  plotted  earlier  in  Chapter  2  as  Fig. 
2.2(b),  computed  from  the  Ambegaokar-Baratoff  results.^  Fig.  4.3(b) 
illustrates  the  same  data  as  Fig.  4.3(a)  but  the  critical  current, 

I c ( T) ,  for  each  junction  Is  normalized  to  the  estimated  zero  temperature 
value,  Ic(0),  which  is  computed  from  the  nominal  IcRn  product  of  1.2  mV 
and  the  sample  resistance,  Rn>  The  nominal  IcRn  product  is  taken  from 
the  average  of  IcRn  products  of  samples  P37-P59  of  Table  4.1.  The  error 
on  the  data  points  on  the  vertical  scale  is  %8%,  coming  from  the  estimate 
for  lc(0).  Fig.  4.3(b)  brings  out  how  the  low  resistance  samples,  i.e. 

R  <  10  $2,  have  IcRn  products  consistent  with  expectations  based  on 
theory,  but  If  R  >  10  Q,  the  IcRn  products  are  low.  This  trend 
continues  for  the  higher  resistance  junctions  actually  tested  in  the 
FIR. 

Table  4.2  displays  the  dc  characteristics  of  the  samples  measured 
in  the  FIR,  the  results  of  which  will  be  presented  in  the  next  section. 
The  analysis  of  the  errors  of  the  various  measured  quantities  of  Table 
4.1  applies  as  well  to  Table  4.2.  The  IcRn  products  vary  from  the 
highest  values  of  %0.9  mV  for  a  junction  of  resistance  of  %16  n, 
decreasing  to  a  value  of  ->.0 - 10  mV  for  junctions  with  resistance  >  300  n. 
This  fall-off  of  the  IcRn  product  with  increasing  junction  resistance  is 
shown  in  Fig.  4.4,  where  the  solid  circles  represent  the  I  and  R  data 
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FIG.  4.3  (a)  The  critical  current  at  a  temperature  T,  normalized  to 
the  critical  current  at  1.4  K,  I  (T)/I  (0),  versus  the  reduced 
temperature,  t«T/T  .  (b)  The  criticarcurrent  at  a  temperature  T, 

normalized  to  the  estimated  zero  temperature  critical  current, 
U0)-1.2±0.1  mV,  [UT)/U0)]„*  versus  t-T/T.. 


of  this  Table.  (See  text.) 


of  Table  4.2.  The  dashed  line  Indicates  the  roughly  constant  IcRn 
products  (xl.2  mV)  of  the  low  resistance  samples.  The  solid  line  and 
open  circles  represent  the  results  of  a  computer  simulation  which  Is 
presented  In  Section  4.5.2.  Critical  current  densities  (est.)  varied 
from  %80  A/ cm*  for  the  1600  Q  sample,  to  about  ^3000  A/ cm*  for  the  16  a 
sample.  Areas  ranged  from  x0.4xl0"8  cm*  to  x2.4xl0"8  cm*. 
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FIG.  4.4  Critical  current,  I_,  vs.  the  normal  resistance,  R  ,  for  the 
junctions  studied  in  this  work.  Data  shown  by  solid  circles  falls  below 
the  dashed  line,  I_R  *1 .2  mV  for  R  >10  &.  The  computer  simulation, 
solid  line  and  open  circles,  closely  follows  the  data.  Error  bars 
indicate  bounds  of  simulated  I„  for  a  given  Rw.  (See  text.) 


Ci  w 


4.3  AC  Josephson  Effect  Analysis 


The  excellent  dc  properties  of  the  junctions  discussed  in  the  last 
section  encouraged  us  to  fabricate  junctions  of  higher  resistances,  and 
to  attempt  to  couple  the  FIR  laser  radiation  to  them,  using  the  FIR 
laser  and  coupling  arrangement  described  in  Chapter  III,  Illustrated  by 
Fig.  3.11  and  Fig.  3.14. 

Typical  I-V  curves  for  increasing  laser  powers  for  two  different 
junctions  are  shown  in  Fig.  4.5.  In  the  low  resistance  large-area 
junction  (Sample  P73,  A%2.2xl0  cm  ,  Rtl6  Q) ,  we  find  flat  steps,  free  of 
noise  rounding,  which  fit  the  RSJ  step  shape  very  well.  The  low 
impedance  of  these  junctions  unfortunately  causes  a  severe  mismatch 
between  junction  and  antenna  so  only  a  small  fraction  of  the  FIR  power 
is  coupled  into  the  junction  and  only  a  few  steps  are  observed.11  In 
the  particular  junction  shown  In  Fig.  4.5(a),  we  observe  four  Josephson 
steps  and  one  photon-assisted  tunneling  step,  when  the  maximum  available 
normalized  voltage  2a*2eVL/fiu)L*1.6  is  coupled. 

The  behavior  of  a  high  resistance  small-area  junction  (Sample  P89, 
R-v.176  fi,  A%0.4xl0  cnr)  is  quite  different  and  is  shown  in  Fig.  4.5(b). 
Here  2a  as  high  as  8.4  was  obtained  from  the  same  laser  power,  and  seven 
Josephson  steps  and  six  photon-assisted  tunneling  steps  were  observed. 

In  this  case  the  steps  are  affected  by  noise,  which  causes  marked 
rounding  and  also  a  slight  tilt  (i.e.  observable  differential 
resistance)  in  the  steps.  The  coupling  of  the  radiation  to  the  junction 
is  good;  it  produces  steps  at  voltages  as  high  as  8.75  mV,  which  is 
comparable  to  the  performance  (12.5  mV)  obtained  by  Weitz  et  al.  , 
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FIG.  4.5  (a)  Typical  I-V  curves  of  a  low-resistance  sample  for 
Increasing  laser  powers.  Here,  va*UPb+Acn)/e.  (b)  Typical  1-V  curves 
of  a  high-resistance  junction  for9increasTHg  laser  powers.  Note  the 
decrease  in  the  critical  current  and  its  subsequent  increase.  The 
voltage  of  the  first  Josephson  step  is  indicated  by  a  dashed  line;  the 
I-V  curves  were  taken  at  power  levels  given  by  the  maximum  of  this  step 
width  and  at  subsequent  zeroes.  The  current  step  at  the  gap  is  also 
Identified  by  dashed  lines. 


using  point  contacts  with  long  wire  antenna  structures  and  the  same 
laser  source.  As  a  consequence  of  the  well -characterized  tunneling 
properties  of  these  junctions  It  is  possible  to  identify  other  features, 
such  as  the  photon-assisted  tunneling  steps,  not  previously  reported  in 
the  FIR  spectral  range.  We  find  six  steps  originating  from  the  gap 
structures  at  Vg«{Apb+A^n)/e  and  -Vg.  They  occur  at  voltages 
V»±Vg+n(1iio^/e) ,  where  n  is  an  integer,  as  expected  from  the  discussion 
of  Section  2.3.  The  structure  of  these  photon-assisted  tunneling  steps. 
Fig.  4.5(b),  will  be  analyzed  analyzed  in  detail  in  the  next  section. 
The  response  of  the  junction  at  4.2  K  has  been  studied  as  well.  At  this 
temperature  the  junction  is  an  SIN  junction  so  the  Josephson  steps  are 
absent,  allowing  us  to  follow  the  variation  of  the  photon-assisted  steps 
with  power. 

Another  simple  qualitative  observation  is  that  the  coupling  of  the 
junction  to  the  FIR  laser  source  is  sensitive  to  the  polarization  of  the 
Incident  electric  field,  the  coupling  being  strong  when  the  polarization 
of  the  FIR  laser  output  is  parallel  to  the  antenna  axis,  the  solid  line 
In  Fig.  4.6,  and  nearly  zero  when  the  polarization  is  perpendicular  to 
the  antenna  axis,  shown  by  the  dashed  line  of  Fig.  4.6.  This 
observation  demonstrates  the  expected  coherent  coupling  of  the  focussed 
beam  into  the  antenna  and  excludes  any  sort  of  diffuse  parasitic 
coupl ing. 

The  variation  of  the  step  widths  with  the  laser  power  can  be 
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studied.  We  have  analyzed  it  in  the  same  way  as  Weitz  et  al.  ,  by 
plotting  the  observed  step  half-widths  (normalized  to  the  critical 


current  in  the  absence  of  radiation)  versus  the  normalized  laser 
radiation  voltage  2a.  To  obtain  2a  from  the  square  root  of  the  power 
level  measured  on  the  reference  thermal  detector,  we  have  chosen  a 
single  scaling  factor  for  all  steps  for  a  given  junction  to  take  into 
account  the  coupling  efficiency  to  that  junction. 

The  frequency-dependent  theory  predicts  that  the  half-width  of  the 
nth  Josephson  step  is  given  by  the  nth  term  of  the  second  sum  in  Eq. 
(2.29),  at  the  voltage  V0*n{fitoL/2e) .  The  result  is 

I„-|£  V«>  (4.2) 

k=-« 


where  for  the  moment  we  ignore  the  photon-assisted  tunneling  effect,  and 
where  wL  is  the  laser  frequency  and  IJ1  is  the  pair  response  function 
defined  by  Eq.  (2.18).  This  response  function  calculated  by  Werthamer9 
for  a  tunnel  junction  with  electrodes  composed  of  two  different 
superconductors  at  T*0  is 
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where  xsWUj*^)  and  * ®  Ia^-^2 I / (Aj-*-^)  >  and  Rp  is  the  resistance  of  the 


tunnel  junction.  To  normalize  Eqs.  (4.3)  to  the  critical  current,  we 
observe  that  ifusQ,  then 


■  'c(0>  ■  iir  (  )  «•>  (4.4) 
n  *12 

which  is  the  expression  computed  by  Ambegaokar  and  Baratoff^,  and 
Anderson^,  Eq.  (2.22b).  For  superconducting  electrodes  made  from 
identical  superconductors, 6  =0,  and  Eq.  (4.4)  reduces  to 


Ic(°)  »  IjjCO)  =  7rA(0)/2eRn 


the  zero  temperature  limit  of  the  expression 

I  (T)R  =  (TrA(T)/2e)tanh(A(T)/2kT),  Eq.  (2.22a). 

v»  ■  • 

Inserting  Eq.  (4.4)  into  Eq.  (4.3)  we  obtain: 


I  I— 1 
'jl  '  e  ' 


L(0) 

~WT 


l 


7 — k((^4)Js) 

(l-x2)*5  l-x2 


x2  a2  k 

K((^*-n 

1-6^ 


(l-fi2)*5 


0SxS6 


<S<X<1 


,  2  vv  2  .2 
\X  -o  )  X  -6 


x>l 


(4.5), 


(4.6) 


where  x=fi m  / (a j+Ag)  and  6  =  |a2-a^ (/(a1+a2)  and  K  is  the  complete  elliptic 
Integral  of  the  first  kind.  Weitz  et  al.12  and  Hamilton22  use  a 
slightly  different  notation  for  the  pair  response  function  In(fiu>  /e). 


They  define  a  new  pair  response  function,  I  which  is  a  function  of  the 
Josephson  frequency  u  0=2eV0/fi.  Their  Ip(  &0)  can  be  obtained  from  our 
/e)  by  replacing  the  dummy  frequency  variable  u>  in  Eq.  (4.6) 
by  half  the  Josephson  frequency  ojq/2.  With  this  change  of  variables 
Eq.  (4.2)  is: 

*7  f  k?.»  Jk{a)  Jn-k(a)  Ip(l2k-nLL)|  (4.7) 
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which  is  Eq.  (7)  of  Weitz,  et  al .  There  is  a  logarithmic  singularity 
(Riedel  peak)  in  I (tvo/e=VQ]  at  the  voltage  corresponding  to  the  gap 
voltage,  i.e.  when  Vq»(a^+ AgJ/e.  This  corresponds  to  a  Josephson 
frequency  u  of  2eVQ/fi«u)0*u)L,  or  *o0 -a»L*C2(A1+A2)/fi3  • 

Eq.  (4.7)  is  easily  evaluated  on  the  computer,  and  we  use  the 
notation  Wn(2a)  for  the  step  half-width  divided  by  the  critical  current 
at  the  normalized  voltage  2a,  ( In/Ic)*Wn(2a) .  The  half-width  of  the  nth 
Josephson  step,  calculated  from  the  RSJ  model  is  Un/Ic)psj*  Jn(2°) 
which  we  obtain  from  Eq.  (2.35).  The  functions  wn(2a)  and  Jn(2a)  are 
compared  with  the  step  width  data  for  the  low  resistance  sample  of  Fig. 
4.5(a)  in  Fig.  4.7. 

We  obtain  good  agreement  between  theory  and  experiment  for  the  n=0 
and  n=2  steps,  but  the  predictions  of  the  Werthamer  theory  and  RSJ  model 
are  indistinguishable  for  the  low  values  of  2a  which  were  attainable. 

For  the  n=l  step,  the  Werthamer  theory  predicts  an  appreciably  larger 
step  than  the  RSJ  model,  but  the  data  fall  somewhat  below  both 
predictions.  This  discrepancy  may  be  related  to  the  apparent  need  to 
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FIG.  4.7  Normalized  step  width  as  a  function  of  the  laser  voltage  in 
the  junction  compared  with  the  Bessel  function  | J  |  dependence  and  the 
Werthamer  theory  |W  |,  for  a  composite  of  data  from  two  low  resistance 
junctions.  The  data  with  arrows  are  from  hysteretic  steps  and  may 
underestimate  the  intrinsic  width  (see  text). 


exceed  a  threshold  power  level  before  the  n*l  step  appears,  whereas  the 
theoretical  dependence  Is  simply  linear  In  2a.  Similar  effects  had  been 
seen  earlier  In  some  point  contact  data.  It  presumably  is  related  to 
the  fact  that  the  voltage  on  the  n»l  step  falls  In  the  hysteretic  part 
of  the  I-V  curve  below  the  energy  gap  voltage.  We  define  the  full  width 
of  the  hyster.tic  step  as  the  current  separation  between  the  jump  up 
from  the  step  in  Increasing  current  and  the  jump  down  from  the  step  in 
decreasing  current-  Premature  jumps  triggered  by  noise,  could  cause 
the  measured  width  to  underestimate  the  intrinsic  width  of  these 
hysteretic  steps.  These  statements  will  be  discussed  more  fully  in  the 
last  section  of  this  chapter.  Section  4.5. 

With  the  high-impedance  junctions,  such  as  the  one  shown  in  Fig. 
4.8,  larger  values  of  2a  were  obtained,  allowing  more  steps  to  be  seen. 
This  permitted  a  more  complete  test  of  the  comparative  quality  of  the 
fit  to  the  Werthamer  theory  and  to  the  Bessel  function  dependence  of  the 
RSJ  model.  However,  the  predictions  of  the  two  models  actually  differ 
rather  little,  because  neither  the  signal  frequency  nor  its  harmonics 
are  very  near  the  Riedel  peak  frequency  for  a  Sn-Pb  junction,  and  the 
scatter  of  the  data  in  these  measurements  limits  our  ability  to 
discriminate  between  the  two.  Still  the  Werthamer  theory  appears  to 
give  a  slightly  better  fit  to  the  observed  positions  of  the  zeroes  of 
the  step  size.  Neither  theory  fits  the  zeroth  step  very  well.  This 
step  goes  to  zero  more  quickly  than  predicted  by  either  theory;  the 
same  was  true  of  the  point-contact  data  of  Weitz  et  al.  It  may  be  due 
to  the  fact  that  the  zeroth  step  width  is  reduced  by  square-law 
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FIG.  4.8  Normalized  step  width  as  a  function  of  the  laser  voltage  in 
the  junction  compared  with  the  Bessel  function  jJ|  dependence  and  the 
Werthamer  theory,  |W|.  The  data  is  from  a  high-resistance  (176  n) 
junction.  The  theoretical  widths  have  been  multiplied  by  a  scale  factor 
of  2.0,  as  described  in  the  text. 


detection  of  noise  at  all  frequencies  In  addition  to  the  effect  of  the 
laser  radiation.  In  making  the  fits  of  the  absolute  step  widths  to  both 
models,  a  scaling  factor  of  a.2.0  was  required.  We  attribute  the  need  for 
this  scaling  factor  to  an  underestimate  of  the  normalizing  denominator 
Ic  because  of  the  above-mentioned  noise  and  fluctuation  effects  In  these 
high-resistance  junctions,  In  which  the  nominal  IcRn  product  Is  well 
below  the  theoretical  value. 

The  power  dependence  of  the  step  half-widths  was  also  measured 
using  the  1.22  mm  FIR  laser  line  on  a  high  resistance  junction  (Rntl56 
a) .  The  radiation  was  coupled  to  the  the  junction  by  an  antenna 
proportionately  longer  than  the  antenna  for  the  496  urn  laser  line.  The 
step  widths  are  plotted  against  jWn|  and  |Jn|  in  Fig.  4.9.  Constant 
scaling  factors  of  2.0  were  also  used  in  this  Figure.  The  quality  of 
the  agreement  between  the  data  and  the  Bessel  and  Werthamer  functions 
was  very  similar  to  that  obtained  with  the  496  urn  line.  Surprisingly, 
the  first  step  was  much  smaller  than  expected  from  either  theoretical 
prediction,  for  unknown  reasons.  The  second  lobe  of  the  power 
dependence  of  the  zeroth  step,  as  well  as  all  the  other  steps,  were  in 
good  agreement  with  the  two  (scaled)  predictions,  which  were  not  all 
that  different  at  the  power  levels  obtained. 

The  constant  scaling  factors,  exceeding  unity,  used  in  the  fits  for 

both  the  1.22  mm  and  the  496  urn  data,  contrasts  with  the  observations  of 
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Weitz  et  al.  ,  on  point  contacts,  in  which  the  corresponding  scale 
factors  were  less  than  unity,  becoming  more  so  with  the  higher  steps, 
presumably  due  to  increasing  heating  effects  at  the  higher  voltages.  As 


FIG.  4.9  Normalized  step  width  as  a  function  of  laser  voltage  in  the 
junction  compared  with  the  Bessel  function  dependence,  |J|,  and  the 
Werthamer  theory,  |W|.  The  data  shown  is  from  a  156  R  junction  irradiated 
by  246  GHz  laser  radiation.  The  theoretical  widths  have  been  multiplied 
by  scale  factors  of  2.0  as  described  in  the  text. 


mentioned  in  Section  2.5,  even  with  all  known  corrections,  the  magnitude 
observed  in  point  contacts  was  about  a  factor  or  2  smaller  than 
predicted  by  the  theory.  It  Is  satisfying  that  our  tunnel  junctions 
give  step  widths  which  fit  the  theory  with  less  need  for  corrections. 
This  might  be  expected  because  the  theory  was  derived  for  a  tunnel 
junction,  not  a  point  contact,  as  shown  in  Section  2.2.  Heating  effects 
are  reduced  in  the  tunnel  case  because  the  voltage  drop  occurs  across  a 
high-resistance  barrier  rather  than  in  a  metallic  constriction  with 
high-conductivity. 

Throughout  this  analysis  we  have  left  out  the  photon-assisted 
tunneling  steps  which  have  a  rather  drastic  effect  on  the  shape  of  the 
I-V  curves,  as  seen  In  Figs.  4.5(a)  and  (b).  In  the  following  section 
we  will  analyze  the  photon-assisted  tunneling  steps,  neglecting  the 
Josephson  steps  for  the  junctions  studied  at  1.4  K,  and  in  the  last 
section  of  the  chapter  we  will  try  to  put  the  ac  Josephson  effect, 
photon-assisted  tunneling  effect,  and  noise  rounding  together  in  the 
analysis. 


4.4  Photon-Assisted  Tunneling  Analysis 


The  photon-assisted  tunneling  effect  comes  from  the  first  term  of 
Eq.  (2.29),  derived  originally  by  Tien  and  Gordon  before  Werthamer,  and 
is  written  below  for  convenience  as  Eq.  (4.8): 

kv0)  - 1  Jk2(a)  W'hr”  (4-8) 

k*-°° 

As  mentioned  in  Section  2.2,  the  Tien-Gordon  result  is  that,  upon 
irradiation  of  the  junction,  the  dc  I-V  curve  will  have  the  form  of  a 


weighted  sum  of  the  dc  I-V  curves  obtained  without  the  radiation,  each 
displaced  in  voltage  by  kffiu^/e),  and  appearing  with  weight  j£(a).  In 
physical  terms,  an  electron  In  one  electrode  with  energy  kfiu>L  less  than 
the  lowest  available  state  in  the  second  electrode  will  tunnel  through 


the  barrier  upon  absorbing  k  photons  of  energy  this  will  increase 
the  current  at  a  voltage  corresponding  to  ±V  +k(fiuL/e) ,  where  V  is  the 
gap  voltage,  i.e.,  A/e  for  the  SIN  case  and  (Aj+A2)/e  for  the  S1 IS2 


methods  gave  results  consistent  to  5%,  we  are  quite  confident  of  our 
internally  calibrated  determination  of  a  from  the  power  meter  reading. 

As  can  be  seen  in  Fig.  4.10,  the  Tien-Gordon  theory  (dashed  lines) 
gives  a  good  fit  to  the  PAT  steps,  and,  although  Eq.  (4.8)  completely 
Ignores  the  Josephson  effect,  it  gives  a  good  account  of  the  entire  I-V 
curve  shape  except  for  the  Josephson  steps  themselves.96 

We  also  applied  similar  analysis  to  data  obtained  at  T«4.2  K,  where 
the  Josephson  effect  is  absent  because  T>T  (Sn).  I-V  curves  for  those 

w 

SIN  junctions,  without  and  with  radiation  are  shown  in  Fig.  4.11.  The 
agreement  between  theory  and  experiment  is  very  good,  so  that  on  the 
scale  of  Fig.  4.11(a),  one  can  not  distinguish  between  the  theoretical 
and  experimental  curves.  To  be  able  to  discern  any  differences,  in  Fig. 
4.11(b)  we  compare  the  experimentally  determined  dV/dl  with  the 
theoretical  expression.  For  a=0,  a  large  peak,  reflecting  the 
nonlinearity  at  V*Vg,  is  obtained,  and  for  a?<0,  minima  and  maxima  occur 
with  period  fiu^/e.  The  position  and  magnitude  of  the  maxima  and  minima 
agree  well  with  the  theory  [dashed  lines  in  Fig.  4.11(b)] 

An  important  parameter  characterizing  detection  efficiency  is  the 
current  responsivity  R,  defined  by  Tucker  as 

R  -  M/(>sVl!l)  (4'9) 

Here  Al  is  the  change  In  the  dc  current  upon  irradiation  of  the 
junction,  and  the  FIR  power  dissipated  in  the  junction  is  1/2VLIL  where 
IL  is  the  current  at  the  FIR  frequency.  This  responsivity  reaches  its 
quantum  limit  e/fiw  when  for  each  photon  absorbed  an  additional  electron 
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displacements  near  the 


tunnels  through  the  barrier.  Before  evaluating  R,  we  compare  the 
measured  Al  with  the  one  computed  using  Eq.  (4.8).  This  is  done  in 
Fig.  4.12  where  we  plot  the  measured  changes  in  current  at  T=4.2  K  for 
a  176  Q  junction  irradiated  at  604  GHz.  Also  shown  is  the  calculated  aI 

p 

versus  a  ,  which  is  proportional  to  the  radiation  power, 

2  2 

l/2(a  /R)(fiu>L/e)  ,  that  would  be  dissipated  in  the  junction  if  it  had 
its  normal  resistance.  It  is  evident  that  the  theory  fits  well  the 
experimental  variations  of  al  with  laser  power.  Similar  results  were 
also  obtained  for  f^=246  GHz  and  for  junctions  of  different  resistances. 

The  theoretical  responsivity  in  the  limit  of  low  power  is  then 
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calculated  from  the  dc  I-V  curve  using  the  Tucker  result 
„  e  Iqp<v*'“/e>  -  2Iqp<v>  +  W'MWe) 


Here,  the  numerator  and  the  denominator  are  proportional  to  the  quantum 

generalizations  of  d'-I/dV  and  dl/dV,  respectively.  In  Fig.  4.13,  we 

plot  R  as  a  function  of  the  bias  voltage  for  a  junction  irradiated  at 

604  GHz.  Using  either  experimental  data  or  the  theoretical  Iqp(V)  for 

T=4.2  K,  which  differ  hardly  at  all,  we  find  that  the  highest  computed 

responsivity  is  0.47  (e/fiuj).  This  maximum  occurs  at  Vg,  where  we  expect 

the  nonlinearity  to  be  greatest.  For  comparison,  we  show,  in  solid 

circles,  the  experimental ly  inferred  values  of  R  found  from  the  inital 

slopes  in  Fig.  4.12,  and  we  also  plot  (dashed  line)  R  as  a  function  of 

V  as  estimated  from  the  data  in  Fig.  4.11  by  assuming  Al  is  linear  in 
2 

a  out  to  o=1.3.  In  both  procedures  we  obtain  the  incident  power  from 
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Sample  P89  fL=604  GHz 

RN~176&  Tboth  ~4.2K 


2 

FIG.  4.12  Change  in  the  dc  current  with  laser  power  (ta  )  at  selected 
dc  bias  voltages.  Points  are  experimental  data;  curves  are  computed 
theoretically  from  the  zero-power  I-V  curve. 


p 

the  value  of  a  by  assuming  that  the  conductance  at  fL  is  given  by  the 

07 

Tucker  relation 

IL/VL“(e/2fi^L)Clqp(V+fituL/e)-Iqp(V--hu>L/e)3. 

No  large  error  should  result  from  this  assunption,  since  this 
theoretical  conductance  never  differs  from  1/Rp  by  more  than  a  factor  of 
2. 

There  is  good  agreement  among  these  various  theoretical  and 
experimental  results.  In  addition,  we  show  curves  for  T*0  and  1.4  K, 
computed  using  the  theoretical  I-V  curves  for  a  SIN  junction;  note  that 
R  reaches  the  quantum  limit  at  T=0  over  a  small  voltage  range.  Clearly 
there  is  an  advantage  in  cooling  the  SIN  junction  to  lower  temperatures ; 
for  example,  at  fL=604  GHz,  the  computed  value  of  R/(e/fiu)  increases 
from  0.47  at  4.2  K  to  0.75  at  1.4  K. 

The  frequency  dependence  of  R,  as  predicted  by  Eq.  (4.10)  for  a 
junction  biased  at  Vg,  is  plotted  in  the  inset  of  Fig.  4.13.  At  T=0,  R 
has  the  full  quantum  limited  value  up  to  the  gap  frequency,  where  fuu=2A. 
At  higher  frequencies  R  decreases  rapidly,  reflecting  the  fact  that  the 
nonlinearity  of  a  SIN  junction  is  limited  to  an  energy  range  -2a.  For 
T>0,  R  falls  below  the  quantum  limit  at  all  frequencies.  There  is  good 
agreement  between  the  theory  and  our  results  at  f^=256  and  604  GHz. 

It  is  worth  mentioning  that  the  remarkable  agreement  between  the 
Tien-Gordon  theory,  which  assumes  a  voltage-biased  junction,  and  our 


findings  is  presumably  due  to  the  presence  of  sufficient  junction 
capacitance  to  cause  a  voltage  biasing  for  ac,  although  the  junction  is 
current  biased  at  dc.  This  effective  voltage  biasing  can  also  explain 
the  agreement  between  the  values  of  a  determined  from  the  Jo§ephson 
effect  data,  and  the  values  of  a  determined  independently  from  the 
photon-assisted  tunneling  data. 


4.5  Noise  Effects 


As  noted  in  Section  4.3,  noise  rounding  effects  are  unimportant  in 
our  lower  resistance,  higher  critical  current-density  junctions.  That 
Is,  the  IcRn  product  computed  using  the  nominal  Ic  is  close  to  the 
theoretical  value  for  a  Sn-Pb  junction,  as  shown  in  Table  4.2  and  Fig. 
4.4.  Moreover,  the  Josephson  steps  are  flat  and  fitted  very  well  by  the 
simple  RSJ  form  without  noise  rounding,  as  is  illustrated  in  Fig.  4.14, 
which  shows  a  step  at  2.5  mV  in  a  16  a  junction.98 

This  situation  is  in  clear  contrast  with  that  which  prevails  in  the 
higher  resistance  junctions  which  give  the  best  coupling  to  the 
radiation.  In  these  the  nominal  I£Rn  product  falls  increasingly  below 
the  theoretical  value  as  Rn  increases,  the  discrepancy  reaching  a  factor 
of  M).l  in  a  381  n  junction ,  as'sbown  previously  in  Fig.  4.4.  The 
systematic  nature  of  the  deviation  suggests  a  fundamental  origin,  and 
we  presume  that  this  is  a  reduction  of  the  apparent  I  by  noise 
rounding. 

Given  the  large  depression  of  the  dc  critical  current  by  noise 
currents,  it  is  obvious  that  they  would  be  expected  to  affect  the  higher 
Josephson  steps  as  well.  For  example,  if  the  effects  of  noise  in 
narrowing  the  finite  voltage  steps  were  significantly  less  than  that  on 
the  zero-voltage  step,  it  could  account  for  the  need  to  scale  up  the 
ratios  In/Ic  by  a  factor  of  \Z  in  Figs.  4.8  and  4.9.  Model  calculations 
reported  below  support  this  hypothesis.  In  addition  to  narrowing  and 
rounding  the  steps,  the  noise  gives  them  appreciable  differential 
resistance,  i.e.  they  are  not  flat. 


In  this  section  we  will  examine  the  effects  of  noise  based  on  two 
distinct  approaches,  first  using  the  analytic  result  of  Lee  appropriate 
to  noise  rounding  on  low  (but  finite)  &c  junctions,  to  determine  by 
curve  fitting  an  approximate  noise  temperature  for  a  representative 
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Josephson  step  on  the  176  fi  junction.  The  remainder  of  the  section  is 
devoted  to  a  computer  simulation  of  a  capacitively  shunted  tunnel 
junction  with  a  noise  current. 

4.5.1  The  Lee  Analytic  Approximation 

OC 

The  Lee  result  is  valid  in  the  limit  of  0^<1.  In  the  limit  of 
zero  capacitance,  it  reduces  to  the  result  of  Ambegaokar  and  Hal  peri n34 
which  was  used  earlier  to  analyze  the  noise  rounding  observed  with  point 
contacts.  Lee's  result  is  displayed  in  Chapter  II  as  Eqs.  (2.49)  and 

(2.50) ,  and  will  not  be  reprinted  here.  The  integrals  and  T2  of  Eq. 

(2.50)  were  evaluated  numerically  using  Simpson's  rule  on  an  LSI-11 
microcomputer. 

Ue  remind  the  reader  that  the  Lee  result  we  use  describes  the  noise 
rounding  of  a  step  in  terms  of  RSJ  model  parameters  appropriate  to  that 
st«*p,  i.e  ,  the  unfluctuated  step  half-width,  1°,  the  dynamic 
resistance  at  the  step  center  in  the  absence  of  radiation,  Rd,  and  so  on 
(see  Section  2.4).  The  usual  ec  parameter  for  the  dc  current  bias  case, 
0c«2eIcR^C/fi,  is  replaced  by  B^*(2eI°RdC)/fi.  Since  0^2  in  our 
junctions,  on  the  4th  step,  at  the  power  level  studied  (2a=4.7),  the 
approximation  for  0^<l  is  not  strictly  applicable,  but  we  assume 
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nonetheless  that  this  approach  will  give  at  least  a  rough  indication  of 
the  effect  of  noise  on  the  steps.  In  fitting  our  data  on  the  4th  step 
at  604  GHz  (i.e.,  5.0  mV),  we  assumed  a  value  for  Rd,  the  dynamic 
resistance  near  the  step,  given  by  the  experimental  curve,  and  used  C, 
Teff*  and  1°  (the  width  of  the  nth  step  in  the  absence  of  fluctuations) 
as  parameters.  For  example,  as  shown  in  Fig.  4.15  we  obtain  reasonable 
fits  for  the  4th  step  in  a  176  n  junction  using  various  combinations  of 
(SO. 01 -0.01 5  pF,  Te^20-25  K,  and  1^2. 2-2. 5  yA,  but  wider  excursions  in 
these  parameter  values  gave  significantly  poorer  fits.  This  fitted 
value  of  is  in  reasonable  accordance  with  the  value  computed  using 
the  Werthamer  theory  (Eq.  (4.7))  and  the  nominal  I_R_  (tl.2mV) 
determined  in  low  Rn  junctions  at  1.4  K.  The  agreement  between  the 
fitted  l|  and  the  one  calculated  from  the  Werthamer  theory  is 
significant  because  the  measured  IcRn  product  is  0.26  mV,  almost  a 
factor  of  5  lower,  and  because  there  was  no  need  to  scale  the  fitted 
step  width  by  the  factor  of  xZ  used  in  Figs.  4.8  and  4.9.  This  result 
supports  our  hypothesis  that  the  scaling  factors  of  xZ  came  from  an 
underestimate  of  the  normalizing  denominator  in  the  step  width  fits,  Ic, 
which  is  reduced  by  noise  and  fluctuation  effects.  The  capacitance  used 
is  based  on  estimates  determined  from  the  sample  geometry  (Table  4.2  and 
Section  4.1). 

These  results  are  rather  interesting  because  the  noise  temperatures 

12 

found  Te^i.20-25  K  are  significantly  higher  than  found  in  earlier  work 
on  point  contacts,  in  which  the  corresponding  temperatures  were 
^eff^lO-lS  K»  a^so  on  the  Ath  Josephson  step  from  the  496ym  laser  line 
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using  the  same  laser  set-up  on  a  point  contact  of  similar  resistance. 

(The  theoretical  result  used  in  that  earlier  analysis  was  the  ec=0  limit 
of  Eq.  (2.49),  which  should  be  appropriate  for  a  point  contact.)  The 
origin  of  this  difference  is  that  in  a  tunnel  .junction,  shot  noise  gives 
a  mean-square  noise  current,  <i^>  *  2elqp  =  4kTeff//Rn’  Teff  = 

eV/2k  a  29K  for  V  =  5.0  mV  (4th  Josephson  step  at  604  GHz),  whereas  for 
a  point  contact,  metallic  heating  gave  Tgff  =  *s(Tbath  +  0  *  10K  for 

the  4th  step  at  5.0  mV  (604  GHz).  Therefore,  the  large  difference  be¬ 
tween  the  two  observed  noise  temperatures  is  due  to  the  different  physi¬ 
cal  processes  which  give  rise  to  the  noise  temperatures. 

4.5.2  Computer  Simulations 

I 

We  turn  now  to  our  computer  simulation  analysis  of  noise  effects  in 
these  junctions.  Our  approach  is  a  phenomenological  one  in  which  we 
take  the  generalized  RSJ  model  presented  in  Eq.  (2.46),  and  solve  it  1 

numerically  on  the  computer  for  the  situations  relevant  to  the 
experimental  results  discussed  in  Section  4.3.  We  reprint  this  Equation 
here  to  aid  in  the  discussion  that  follows:  ! 

2 

— J-  +  g(jjj-)  +  sin$  =  prfc  +  pl  sinn-c  +  pn  (4.11 ) 

dT 

where  g( d^/dx ) ,  the  nonlinear  tunneling  current  normalized  to  the 
critical  current,  is  a  function  of  the  normalized  frequency  d$/dt  .  On 
the  right  hand  side  of  Eq.  (4.11)  prfc  is  the  normalized  dc  bias 
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current,  the  time  t  is  normalized  to  the  inverse  plasma  frequency 
pLsinfir  is  the  ac  current  source,  and  pn  is  a  Langevin 
noise  term.  Although  it  is  possible  to  use  the  full  nonlinear 
quasiparticle  tunneling  curve  to  model  g(d4>/dx),  we  have  chosen  instead 
to  approximate  it  as  follows:  for  voltages  below  eV=ASn+Apb  we  use  an 
adjustable  linear  leakage  resistance  RL  which  models  the  usual 
experimental  observation  of  excess  leakage  current  below  the  gap  rise. 
Above  eV=A<jn+Apb  we  use  the  full  BCS  quasiparticle  I-V  curve  for  T=0, 
since  we  will  be  considering  experimental  results  far  below  Tc.  In  this 
manner  execution  times  are  made  faster,  but  the  results  of  the 
simulation,  as  we  will  show,  provide  adequate  quantitative  comparison  to 
the  experimental  results. 

The  noise  term  on  the  right  hand  side  of  Eq.  (4.11),  pN,  is 
treated  in  a  white  noise  approximation  with  the  following 
auto-correlation  function: 


<p^(t 1 )  pm(t+t' )>  = 


N 


B, 


(*/e2) 


P-  coth 
dt 


n  y  d4>> 
2(fi/e2)  dT 


6(t) 


(4.12) 


2 

where  Bc*2eIcR‘C/fi,  Y=fil c/ekTN ,  and  d$/dx  is  the  dimensionless  rate  of 
change  of  the  phase,  while  6(t)  is  a  dimensionless  5  function.  In  this 
calculation  we  are  assuming  that  the  fluctuations  are  coming  from  the 
normal  resistance..  Rn,  i.e.  we  approximate  the  quasiparticle  part  of 
the  I-V  curve,  I^V)  by  V/Rn ,  for  the  noise  term  (i.e.  on  the 
right  hand  side  of  Eq.  (4.11)).  This  is  in  contrast  to  our  use  of  a 


non-linear  quasiparticle  resistance  on  the  left  hand  side  of  Eq. 

(4.11).  Eq.  (4.12)  is  the  current-current  auto-correlation  function  in 
normalized  units,  coming  from  the  power  spectrum  of  Eq.  (2.44),  that  is 


<fN(t)  IN(t+t)>  =  2eIqp(V)  coth  (^r)  6(t) 


(4.13) 


where  V=V(t)  is  the  instantaneous  voltage,  and  V*(ii/2e) dtf/dt.  Eq. 

(4.13)  is  converted  to  Eq.  (4.12)  by  dividing  both  sides  of  the  Eq.  by 
2  -1 

l£,  and  by  dividing  the  6  function  by  (ojj)  to  obtain  a  dimensionless 
6  function,  besides  assuming  Iqp(V)=V/Rn. 


Eq.  (4.11)  is  solved  numerically  for  the  case  of  a  current  sweep 
by  using  a  4th  order  Runge-Kutta  algorithm100,  which  we  describe  only 
briefly  here.  The  quantities  we  wish  to  solve  for  are  the  phase,  0  at 
a  time  t,  which  we  denote  by  V^(t)*0(t),  and  the  rate  of  change  of  the 
phase,  V2(t)=0(t),  given  a  dc  bias  current,  pdc,  rf  drive  current 
p^sintit  ,  and  noise  current,  p^.  The  method  we  use  involves  splitting 
the  2nd  order  equation  of  Eq.  (4.11)  into  two  simultaneous  first  order 
equations,  and  then  applying  the  Runge-Kutta  algorithm  for  each.  The 


two  first  order  equations  are 

Vt  *  f!  =  f1(V1-0,t)  =  V2  (4.14) 

^2  *  f2  =  MV1  *V2s^,t)  =  °dc  +  pl  sinnx  +  PN  '  "  sinVl 
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The  equations  for  and  V2  are  solved  for  by  the  Runge-Kutta  algorithm, 
which  we  illustrate  with  V2<  We  take  r^ih,  x^^x^+h,  where  i  is  an 
integer  greater  than  zero  so  that 

^2 ( Ti  +1 )  =  ^  +  6^kl  +  ^k2  +  ^k3  +  k4^  (4.16) 

where  the  are  given  by 

kl  =  f2^Ti »  vi (Ti ) »  V2(Tn*)) 

k2  s  f2(xi+iSh,  V-j (xi+Jsh) ,  V2(x.)  +  »5hk1 ) 

(4.17) 

k3  *  f2(xi+*s h,  V^x^),  V2(x.)  +  Jshk2) 
k4  =  f2  (Ti+h*  Vl(Ti.+h),  W  +  hk3} 

A  similar  set  of  equations  applies  to  V^=0(x).  The  current  is 
incremented  in  fractions  of  the  critical  current,  usually  Apdc=0.01. 

(A  value  of  Pdc=l  is  the  critical  current  in  normalized  units.)  At  the 
initial  current  level,  i.e.  pdc=0-01»  the  algorithm  starts  with  xo=0, 
and  with  the  initial  conditions  for  the  phase  and  derivative  of  the 
phase,  V1=jJ(0)=0,  and  V2=0(O)=O.  The  time  increment,  h,  is  in  units 
of  the  inverse  of  the  plasma  frequency,  and  is  typically  0.04.  The 
algorithm  waits  for  a  period  of  time  to  permit  initial  transients  to  die 
out,  and  then  the  average  voltage  is  calculated  by  the  average  of  V2 
divided  by  (0C)^,  i.e.  V/IcRn=(Bc)“^^  dtf/dx.  After  this  the 
current  is  incremented  and  the  algorithm  calculates  0=V.|  and  0=V2  at 
the  new  current  point,  using  the  last  value  of  and  V2  from  the 
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previous  current  point.  The  total  averaging  time  for  each  current 
point  in  units  of  inverse  plasma  frequency,  (wj)""'  ,  is  given  by 

Tavs[2ir/flh][Ttota^-Twait] ,  typically  «vl04-105,  where  n=cuL/ajj  as  before, 
and  h  is  the  time  increment. 


The  noise  term,  p^,  is  calculated  at  each  point  in  time  in  the 
following  way.  First,  the  average  is  taken  of  12  numbers,  xi ,  selected 
from  the  interval  [0,1]  by  a  pseudo- random  number  generator.  This 
procedure  produces  an  approximately  gaussian  distribution  with  mean  of 
1/2  and  width  1/12.  The  mean  is  shifted  to  zero  by  subtracting  1/2  from 
this  average,  and  the  width  is  adjusted  to  give  the  required 
mean-squared  current  noise  from  £q.  (4.12)^ ,101 .  In  this  discrete 

version  then  the  noise  current  at  a  time  is  denoted  by  (p^).. : 


(eN>i  * 


2  "  R"  '  it  coth  -Z-P-) 


(cf-H 


-) 


'h6c'W)'  dT 


2(*/e2)  7ITdT 


12 

x  (  £  (xk-*s)) 
k=l  K 


(4.18) 


The  discretizing  of  the  time  variable  causes  the  6  function  6 (t )  to 
become,  6-*-(l/h)6i  Q ,  where  6..  Q=1  if  i=o  (t1-=to=0)  ,  and  is  zero 
otherwise,  where  h  is  the  size  of  the  time  step. 


We  can  now  examine  the  results  of  the  simulation  beginning  with  the 
dc  case.  The  method  summarized  by  Eqs.  (4.11)-(4.18)  is  applied  now  to 
the  situation  where  there  is  no  rf  current,  only  dc  current  and  noise, 
to  determine  the  effect  of  noise  on  the  critical  current.  For  a  noise 
tenperature  of  2  K,  Fig.  4.4  shows  that  this  computation  reproduces  the 
observed  fall-off  of  the  IcRn  product  with  increasing  Rn<  The  solid 
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circles  in  Fig.  4.4  are  the  data  from  Table  4.2,  while  the  open  circles 
are  the  results  of  the  simulation.  For  low  resistance,  low  Bc  junctions, 
the  IcRn  products  are  approximately  constant  shown  by  the  dashed  line  in 
the  Figure.  On  the  other  hand,  for  the  high-resistance  junctions 
measured  the  critical  currents  are  considerably  lower  than  expected  from 
the  behavior  observed  with  the  low  resistance  junctions,  in  agreement 
with  the  simulated  behavior  shown  by  the  open  circles  and  the  solid 
line.  The  error  bars  on  the  calculated  points,  at  the  high  R  end  of  the 
curve,  are  equal  to  the  variance  of  the  I c * s  from  the  5  runs  with  the 
simulation.  Note  that  the  agreement  is  good  for  junctions  ranging  in 
resistance  over  two  orders  of  magnitude.  The  deviation  of  the  simulated 
curve  from  the  data  may  be  due  to  the  relatively  short  averaging  time  of 
the  simulation  (->.10  ysec  to  sweep  from  zero  current  to  the  critical 
current)  relative  to  that  on  the  laboratory  scale(  ->.10  secs  to  sweep  from 
zero  current  to  the  critical  current).  Therefore,  the  fast  sweep  rate 
used  in  the  simulation  may  underestimate  the  effect  of  the  noise.  It  is 
assumed  that  the  capacitance  does  not  vary  from  junction  to  junction  in 
the  simulation,  whereas  in  fact  the  capacitances  do  vary  by  a  factor  of 
■>.5,  which  may  also  explain  the  observed  deviation  between  the  simulation 
and  the  data. 

The  simulation  produces  an  I-v  curve  with  the  shape  shown  in  Fig. 
4.16(a)  closely  reproducing  the  observed  switching  behavior  in  Fig. 
4.16(b).  The  dashed  lines  in  Fig.  4.16(a)  indicates  the  computed  I 
for  this  junction;  the  variance  on  the  computed  I  is  given  by  the 
dotted  lines  in  the  Figure.  In  the  simulation  the  switching  from  the 
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FIG.  4.16  (a)  Simulated  I-V  curve  with  noise  but  no  ac  currents. 
Switching  at  critical  current  shown  by  dashed  lines.  Critical  currents 
from  simulation  (for  5  runs)  fall  within  bounds  indicated  by  dotted 
lines  (see  text),  (b)  I-V  curve  from  the  176  n  junction  (data).  The 
noise  temperature  used  in  the  simulation  is  TN=2K,  while  the  bath 
temperature  is  1.4K.  ” 
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non-zero  voltage  state  to  the  zero- voltage  state  occurs  at  a  higher  cur¬ 
rent  level  than  the  data,  perhaps  due  to  the  assumed  leakage  resistance 
of  10  times  the  normal  resistance,  which  may  be  too  low,  as  may  the 
assumed  noise  temperature  of  2K. 

We  now  simulate  the  effects  of  the  laser  radiation  on  the 
junctions.  Fig.  4.17(a)  shows  the  simulated  critical  current  with 
laser  power  divided  by  the  zero  power  critical  current  (in  the  absence 
of  noise),  while  Fig.  4.17(b)  displays  the  full-width  of  the  first 
step,  also  in  the  absence  of  noise.  Fig.  4.17(b)  is  labelled  with  I  . 
where  the  step  width  is  measured  only  by  the  distance  in  current  from 
the  point  a  to  the  point  b,  shown  by  the  inset  of  the  Figure,  i.e.  the 
portion  of  the  step  on  the  upsweep.  Idb  labels  the  maximum  width  of  the 
step,  the  region  in  current  from  d  to  b  in  the  inset  of  the  Figure.  We 
find  for  the  finite  capacitance  case,  p  -».6-7,  studied  in  the  Figure,  that 

V 

the  step  widths,  Idb,  are  adequately  described  by  the  Bessel  function 
step  width  formula  used  before  for  a  voltage  biased  (dc  and  ac)  RSJ 
model  junction: 


V'c a  V2“> 


(4.19) 


for  the  step  half-widths.  In  order  to  make  a  comparison  with  the 
experimental  data  or  the  simpler  voltage  biased  RSJ  model,  a 
relationship  must  be  found  between  2a=2eVb/fi  ,  the  normalized  ac 
voltage  in  the  junction  (e.g.  fitted  from  the  data),  and  I|_/Ic,  the 
normalized  ac  current  used  in  the  RSJ  model  simulation.  For  the  finite 
capacitance  case  a  simple  relationship  between  the  two  has  been  found: 


FIG.  4.17  Zeroth  and  first  step  simulated  behavior  for  a  junction  with 
rf  and  dc  bias  but  no  noise,  (a)  Zeroth  step  half-width  normalized  to 
the  critical  current,  (b)  First  step  full-width  normalized  to  Ic.  Note 
that  if  the  width  of  the  upsweep  is  taken,  I  b,  the  width  is 
underestimated,  whereas  Idb,  the  full:width  $  the  step  including  the 
downsweep  reproduces  the  Bessel  function  behavior. 
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(4.20) 


o 

where  n=u)L/u>j,  and  0c*(2eIcRji|C)/'fi,  as  before.  In  Fig.  4.17  the  Bessel 

functions,  JQ(2a)  for  the  zeroth  step,  and  2^(20)  for  the  first  step 

full-width,  both  normalized  to  the  critical  current,  are 

indistinguishable  from  the  lines  drawn,  for  2a*( IL/Ic)/6.73.  This 
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result  was  predicted  by  Braiman  et  al.  for  a  resistively  shunted 
junction  model  with  G(V)*l/Rn,  for  the  limit  of  n»l.  This  latter 
condition  is  only  marginally  satisfied  by  the  sample  in  Fig.  4.17. 

Fig.  4.17(b),  together  with  the  experimental  results  of  Fig.  4.7, 
indicate  that  it  is  the  width  of  the  fully  hysteretic  step  that  is  of 
importance,  justifying  the  procedure  used  earlier  for  the  analysis  in 
Section  4.2  and  Fig.  4.6. 


If  noise  is  included  in  the  simulation  of  the  zeroth  step,  at  the 
maximum  of  the  first  hump  of  the  Bessel  function,  JQ,  that  is  at  2a=4.0, 
or  1l/Ic%27,  the  step  is  rounded  somewhat.  Fig.  4.18  plots  the 
simulated  zeroth  step  I-V  at  this  power  level.  Plotted  with  the 
simulation  is  the  computed  noise-reduced  critical  current  at  zero  laser 
power.  The  error  bars  indicate  the  limits  on  the  computed  critical 
current.  With  the  2  K  noise  temperature  chosen,  the  step  width  divided 
by  this  critical  current  is  about  0.5,  which  is  roughly  the  expectation 
for  the  step  width  ratio  without  noise.  If  the  noise  temperature  used  is 
higher  (in  fact  up  to  \10  K),  enough  to  suppress  the  fluctuated  critical 


FIG.  4.18  Simulated  step  at  2a=4.0  Arrow  labeled  by  I  is  the 
unfluctuated  critical  current  estimated  from  the  nominal  I  R  product 
(see  text).  The  critical  current  with  noise  is  also  shown,  as  well  as 
the  estimated  errors  on  the  simulated  critical  current,  i.e.,  the  error 
bars  drawn. 
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current  closer  to  the  measured  value  or  the  averaging  time  in  the 
simulation  is  longer,  the  size  of  the  zeroth  step  at  2a=4  does  not 
change  significantly,  even  though  the  critical  current  at  a=0  is  reduced 
much  more.  The  simulated  critical  current  without  laser  power  can  be 
reduced  enough  so  that  the  ratio  of  the  noise  rounded  zeroth  step-width 
at  2a*4.0  to  the  noise  rounded  critical  current  is  roughly  unity.  This 
would  justify  the  scaling  factor  of  ^2.0  for  the  Bessel  and  Werthamer 
functions  used  in  the  analysis  in  Section  2.2. 

We  can  now  attempt  to  model  the  effects  of  noise  on  the  complete 
I-V  curve  of  a  junction  at  a  particular  laser  power.  We  choose  2a=4.0, 
the  power  level  used  in  the  Fig.  4.18.  We  create  a  simulated  tunnel 
junction  I-V  curve  complete  with  photon-assisted  tunneling  steps  and 
noise  rounded  ac  Josephson  steps  in  the  following  way.  Using  the 
simulation,  we  create  a  noise  rounded  I-V  curve  with  only  Josephson 
steps.  We  subtract  off  the  smooth  background  as  a  function  of  voltage 
from  this  I-V  curve,  leaving  only  the  local  current  displacements  from 
the  ac  Josephson  steps.  A  simulated  I-V  curve  with  only  photon-assisted 
tunneling  steps  is  created  from  the  dc  quasi  particle  I-V  curve  for  this 
particular  choice  of  2a  following  the  procedure  described  in  Section 
4.4.  This  I-V  curve  is  very  similar  to  the  dashed  curve  for  the  176  n 
junction  at  a*2.1  shown  in  Fig.  4.9.  The  local  current  displacements 
from  the  ac  Josephson  effect  are  added  as  a  function  of  voltage  to  the 
I-V  curve  with  only  the  photon-assisted  tunneling  steps,  creating  an  I-V 
curve  with  both  photon-assisted  steps  and  ac  Josephson  steps.  The 
results  of  this  procedure  are  shown  in  Figs.  4.19(a)-(c).  In  Fig. 
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4.19(a)  we  show  the  experimental  I-V  curve  at  o*2  we  are  simulating. 

The  simulated  I-V  curve  is  shown  in  in  dashed  lines  next  to  the 
experimental  I-V  curve  (Fig.  4.19(b)),  while  Fig.  4.19(c)  plots  the 
simulated  I-V  curve  (dashed  lines)  on  top  of  the  experimental  I-V  curve 
shown  in  the  solid  line?.  It  is  remarkable  how  similar  the  simulation 
is  to  the  actual  data.  The  differences  between  the  simulation  and  the 
data,  particularly  at  the  3rd  step,  are  most  likely  due  to  the  fact  that 
the  simulation  does  not  include  any  effects  of  the  frequency  dependence 
of  the  step  widths.  That  the  two  effects  should  add  in  current  at  each 
voltage  is  precisely  the  expected  behavior  from  the  Werthamer  result  for 
a  voltage  biased  junction,  i.e.,  as  in  Eq.  (2.29),  further  confirming 
that  the  finite  capacitance  of  the  junction  makes  it  appear 
approximately  voltage  biased. 

We  can  apply  the  above  procedure  to  the  situation  of  the  4th  step, 
analyzed  previously  with  the  simple  result  of  Lee  (Fig.  4.14),  which 
ignores  photon-assisted  tunneling.  Fig.  4.20  displays  the  results  of 
the  procedure.  All  four  curves  show  the  4th  step,  measured  from  the 
center  of  the  step,  with  the  data  shown  as  a  solid  line,  and  the 
simulation  as  the  solid  circles.  The  power  level  is  2a=4.0  for  the 
upper  curves,  while  2a*5.75  for  the  lower  curves.  Note  that  while  the 
upper  left  I-V  curve  fits  the  data  very  well,  the  lower  left  simulated 
curve  lies  below  the  data.  The  difference  at  the  higher  power  level 
could  be  due  to  the  presence  of  a  photon-assisted  tunneling  step 
distorting  the  I-V  curve.  To  test  this  conjecture,  on  the  right  half  of 
Fig.  4.20  are  I-V  curves  measured  from  the  step  center  where 
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FIG.  4.19  (a)  I-V  curve  of  a  176  n  junction  at  a=2.0;  Josephson  steps 
are  labelled  by  the  integers  n=2 ,  n=3,  and  n=4.  (b)  Simulation  ot  the 
data  shown  in  (a),  which  includes  noise  and  photon-assisted  tunneling 
steps  (see  text),  (c)  Simulation,  dashed  lines  in  (b),  compared  with 
data  from  (a)  illustrating  the  excellent  quality  of  the  fit. 
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simulations  now  include  the  effects  of  the  photon-assisted  tunneling 
steps  in  addition  to  the  noise.  The  resulting  fit  to  the  step  data  with 
2a=4.0  is  worse,  the  simulation  lying  below  the  data,  while  for  2a=5.75 
the  simulation  lies  much  closer  to  the  data,  the  net  result  being  that 
both  fits  are  of  roughly  comparable  quality. 

It  is  rather  remarkable  that  we  have  obtained  such  good  fits  to  the 
data  for  both  the  fall-off  of  the  IcRn  with  increasing  junction  resis¬ 
tance  and  the  rounding  of  the  Josephson  steps  with  a  model  as  simple  as 
the  RSJ  model  that  we  have  used.  We  have  made  two  simplifications  which, 
on  the  face  of  it,  do  not  appear  to  be  physically  correct.  First,  we 
have  taken  the  voltage  in  the  noise  term,  p^,  of  Eq.  (4.11-4.13)  to  be 
the  instantaneous  voltage.  This  would  seem  to  be  a  problem  because  the 
noise  power  spectrum  from  which  the  noise  term,  pN,  was  calculated,  was 
derived  for  the  case  of  a  constant  dc  bias  voltage,  so  that  it  would 
seem  appropriate  to  use  an  average  voltage  to  calculate  the  noise  term. 

However,  we  replaced  the  voltage-dependent  (shot)  noise  term  in  the 

2 

computer  simulation  (<iN>  'v  2eV/Rn)  by  a  Johnson  noise  term 
2 

(<iN>  =  4kTeff/Rn)  at  an  elevated  temperature,  Tgff  =  eV/2k  -v-  30K, 
appropriate  to  the  4th  step  at  5.0  mV  (604  GHz),  and  found  that  the 
simulated  step  shape  with  this  noise  temperature  was  nearly  identical 
to  that  simulated  before  with  the  voltage-dependent  noise  term.  Also, 
V(t)  *  V  on  the  4th  step,  since  I»IC»  so  for  this  reason  the  instan¬ 
taneous  voltage  should  be  no  problem  in  the  simulation.  Second,  we  used 
the  normal  resistance  to  model  the  noise  below  the  gap  which  contrasts 
with  the  theoretical  expectation  from  Section  2.4  of  using  the  sub-gap 
leakage  resistance  to  model  the  noise  for  voltages  below  the  gap.  The 
junctions,  however,  are  well  coupled  to  external  fields  through  the 


antenna,  which  has  a  resistance  of  roughly  the  normal  resistance,  so 
the  fact  that  we  are  using  the  normal  resistance  for  the  noise  below 
the  gap  may  not  be  a  bad  approximation  to  reality.  We  have  also  tried 
using  the  sub-gap  leakage  resistance  for  the  noise,  and  for  the  typical 
sub-gap  leakage  resistance  of  10Rn  we  found  that  we  could  model  the 
reduction  in  the  critical  current  with  a  noise  temperature  of  ^30  K. 

This  noise  temperature  seems  unreasonably  high.  In  any  case  we  have 
found  rather  remarkable  agreement  between  our  data  and  the  RSJ  simula¬ 
tions  used,  which  gives  some  credibility  to  the  method. 

4.5.3  Discussion 

Before  concluding  this  Chapter  we  would  like  to  pause  briefly  to 
point  out  some  recent  analytic  approaches  to  the  problem  of  noise  in 
junctions,  which  go  beyond  the  now  classic  results  presented  in  Section 
2.4. 

103 

In  a  recent  paper,  Ben-Jacob  et  al .  ,  have  examined  the  effects 

of  shot  noise  on  the  non-zero  voltage  state  of  a  hysteretic  junction  and 
found  a  distribution  of  fluctuations  that  is  non-symmetric  in  that 
positive  fluctuations  are  more  probable  than  negative  ones.  This  is  not 
surprising  since  the  shot  noise  power  spectrum  shows  increased  noise 
power  with  increasing  voltage.  Their  results  are  confined  to  the 
hysteretic  region  of  the  junction  I-V  curve,  and  are  applicable  for 
0  >1.6,  but  not  too  large,  and  can  be  tested  by  observing  the  transition 
rate  to  the  zero  voltage  state  from  the  non-zero  voltage  state. 
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Buttiker,  Harris  and  Landauer104  have  extended  the  treatment  of 
the  noise  activated  escape  of  a  particle  out  of  a  metastable  well 
(e.g.  the  transition  out  of  the  zero  voltage  state  of  a  junction  as  in 
Fig.  2.4),  given  originally  by  Kramers,  and  applied  to  superconducting 
tunnel  junctions  by  Lee  for  the  extreme  underdamped  case  (Bc»l),  to  a 
broader  range  of  Bc's.  Their  results  are  applicable  to  shunted  junc¬ 
tions  with  &c>44,  a  condition  satisfied  by  only  three  of  the  seven 
junctions  studied. 

The  theoretical  developments  cited  above  are  applicable  to  shunted 
junctions  and  only  in  certain  limits  of  parameters.  These  analyses  have 
not  been  extended  to  the  case  of  a  junction  with  the  full  non! inearity 
of  the  quasiparticle  I-V  curve,  nor  a  junction  which  is  driven  by  a  FIR 
or  microwave  source.  The  most  encouraging  approaches  are  those  recently 
developed  by  Caldeira  and  Leggett105,  Chakravarty106 ,  and  Ambegaokar 
Eckern,  and  Schon10^.  These  approaches  were  developed  for  the  problem 
of  the  transition  out  of  the  zero  voltage  state  where  thermally 
activated  processes  are  frozen  out,  and  depend  on  complicated  functional 
integral  formalisms.  These  treatments,  however,  can  take  the  full 
nonlinearity  of  the  quasiparticle  I-V  curve  into  account,  but  physically 
intuitive  results,  or  results  for  junctions  driven  by  high  frequency 
sources  are  not  yet  available. 

We  conclude  this  section  and  the  chapter  by  noting  that  we  have 
found  rather  remarkable  agreement  between  the  theory  of  the  tunnel  junc¬ 
tion  described  in  Chapter  II,  and  our  data,  at  the  highest  frequency  yet 
reported  for  a  thin-film,  fabricated  superconducting  tunnel  junction. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

In  this  chapter,  we  will  summarize  the  main  results  of  this  work, 
as  well  as  briefly  discuss  experiments  which  could  be  done  in  the  future 
to  extend  this  work  or  to  resolve  some  of  the  questions  which  are  at 
present  still  open. 

We  have  photol ithographically  fabricated  small-area  Sn-SnO-Pb 
superconducting  tunnel  junctions  which  have  nearly  ideal  dc  character¬ 
istics.  We  havt  coupled  submillimeter  laser  radiation  to  these  junc¬ 
tions  by  fabricating  them  at  the  centers  of  dipole  antennas  which  are 
resonant  at  the  wavelength  used  (X=496  pm  or  x=1.22  mm). 

In  a  high  resistance  junction  (Rn~176  fi),  with  a  submicron  area 

p 

(A-0.4  (pm)  ),  we  have  coupled  enough  radiation  to  observe  seven  ac 
Josephson  steps,  up  to  a  voltage  of  8.75  mV,  a  result  comparable  to  that 
observed  in  point  contacts  (-12.75  mV).  In  addition  to  this,  six  photon- 
assisted  tunneling  steps  have  been  observed  on  this  same  junction,  an 
effect  not  observed  in  previous  work  with  point  contacts,  although  it  is 
expected  from  the  microscopic  (Werthamer)  theory  of  the  tunnel  junction. 
We  have  systematically  studied  the  power  dependence  of  the  widths  of  the 
Josephson  steps  at  X=496  pm  (604  GHz)  and  at  X=1 .22  mm  (246  GHz)  and 
find  reasonable  agreement  with  both  the  RSJ  model  and  the  Werthamer 
(voltage-biased  frequency-dependent)  theory,  which  are  not  all  that 
different  at  the  laser  frequencies  tested  thus  far. 

In  addition,  we  have  systematically  studied  the  photon-assisted 
tunneling  effect  at  246  and  604  GHz  at  1.4  K  (T<T.(Sn))  where  the 
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Josephson  effect  is  present,  and  also  at  4.2  K  where  it  is  absent.  We 
find  good  agreement  with  the  Tien-Gordon  (Werthamer)  theory.  We  infer 
that  the  responsivity  of  these  junctions  approaches  the  quantum  limit  at 
low  temperatures  for  photon  energies  less  than  2a. 

We  have  studied  the  effects  of  noise  on  the  behavior  of  these  junc¬ 
tions.  We  have  observed  that  noise  systematically  reduces  the  measured 
IcRn  product  (at  1.4  K)  with  increasing  junction  resistance.  Noise  also 
rounds  and  produces  a  slight  tilt  on  the  ac  Josephson  steps  of  high 
resistance  tunnel  junctions.  The  Josephson  step  rounding  has  been 

interpreted  in  two  ways:  First,  using  a  simple  analytic  result  of  P.  A. 
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Lee,  applicable  to  low  a  junctions,  based  on  the  RSJ  model  of  the 
tunnel  junction  with  reasonable  estimates  of  junction  parameters,  we 
have  found  good  fits  to  the  Josephson  step  rounding  with  noise  tempera¬ 
tures  of  roughly  20-25  K,  which  is  thought  to  reflect  the  shot  noise 
which  is  not  included  in  this  model.  Second,  we  have  developed  a  com¬ 
puter  simulation  to  study  the  effects  of  noise  on  the  tunnel  junctions. 
This  simulation  uses  an  RSJ  model  which  has  been  modified  to  include 
the  nonlinear  quasiparticle  I-V  in  the  junction  dynamics.  The  simula¬ 
tion  includes  ac  current  sources,  and  a  noise  current  source  which,  for 
simplicity,  is  that  expected  theoretically  for  the  junction  in  the 
normal  state,  and  models  shot  noise  behavior  at  high  voltages  and  ther¬ 
mal  noise  at  low  voltages.  The  simulation  produces  dc  I-V  curves  which 
look  quite  close  to  the  measured  I-V  curves.  Using  a  noise  temperature 
of  ^2  K,  we  reproduced  the  observed  fall -off  of  the  IcRn  products  with 
increasing  junction  resistance.  The  noise  rounding  behavior  of  the 
steps  is  also  modelled  well  by  the  simulation,  if  the  effects  of  the 


photon-assisted  tunneling  steps  are  included. 

The  general  agreement  between  our  data  and  the  theoretical  predic¬ 
tions  of  the  frequency-dependent  (Werthamer)  theory  of  superconductivity 
is  excellent.  None  the  less,  there  is  a  need  to  extend  this  work  beyond 
its  present  limits.  An  original  goal  of  the  work  was  to  measure  the 
frequency  dependence  of  the  ac  Josephson  effect  above  the  superconducting 
energy  gap.  We  have  had  partial  success  towards  this  goal,  since  we  have 
successfully  coupled  FIR  laser  radiation  into  the  junction  atdio)L/2e= 

0.51  mV,  and  T>uiL/2e=l  .25  mV  from  the  246  GHz  and  604  GHz  laser  lines, 
respectively.  The  higher  frequency  is  still  only  -  2/3  of  the  gap.  Two 
significant  questions  still  remain.  What  is  the  frequency  dependence  of 
the  Josephson  effect  above  the  gap?  Secondly,  will  the  simple  Tien- 
Gordon  (Werthamer)  photon-assisted  tunneling  result  still  work  as  well  at 
1-2  THz  as  it  does  at  0.6  THz? 

Now  that  we  have  shown  that  it  is  possible  to  couple  submillimeter 
laser  radiation  into  a  thin-film,  photolithographically  fabricated 
junction,  one  should  attempt  to  couple  radiation  into  a  single  junction 
over  a  wide  range  of  FIR  laser  lines.  To  do  this  we  must  have  an 
antenna  which  has  a  response  which  is  substantially  frequency- independent. 

The  bow-tie  antenna  is  such  an  antenna.  This  antenna  has  been  shown 
to  work  well  in  the  frequency  range  of  interest  to  us  (0.246-2.5  THz)6.0,65 
The  bow-tie  antenna  has  a  constant  input  impedance  over  the  desired 
frequency  range,  and  being  a  planar  antenna  is  easily  fabricated  photo¬ 
lithographically.  The  impedance  of  the  bow-tie  antenna  depends  on  bow 


angle;  it  is  lower  for  larger  bow  angles,  as  is  shown  in  Fig.  5.1.  The 
theory  for  the  antenna  impedance  (solid  line)  agrees  well  with  values 
measured  (crosses)  in  a  10  GHz  microwave  modelling  experiment.®®  A  bow 
angle  of  -90°  gives  an  almost  purely  resistive  impedance  of  -100  a  if  it 
were  an  a  fused  quartz  substrate.  As  shown  in  this  work,  it  is  feasible 
to  fabricate  a  small -area  junction  of  roughly  this  impedance,  for  good 
matching. 

A  second  improvement  would  be  to  successfully  fabricate  junctions 
of  smaller  area  and  higher  current-density  than  fabricated  thus  far  in 
our  lab.  All  but  one  of  the  junctions  measured  in  the  FIR  were  fabri¬ 
cated  using  the  resist-aligned  technique  discussed  in  Chapter  3.  We 
attempted  to  fabricate  smaller  area,  edge-aligned  junctions  with  the  Sn- 
Pb  metal  system  which  worked  so  successfully  with  the  resist-aligned 
junctions,  but  we  had  relatively  little  success.  These  fabrication 
o-oblems  may  be  due  to  some  metallurgical  problem  which  is  not  understood 
by  us  at  present.  The  simplest  course  of  action  would  be  to  try  another 

metal  system,  such  as  Pbln-Pb  (20  weight  percent  In)  which  was  used  with 
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good  success  at  Bell  Labs.  The  smallest  area  junction  which  could  be 

made  using  the  edge-aligned  technique  and  our  present  day  2  ym  linewidth 

contact  lithog.aphy  is  -0.2  (ym)  .  A  junction  of  reasonable  resistance 

(Rn~20-50  a)  could  then  be  expected  to  have  current  densities  in  the 
4  2 

range  of  -5-2x10  A/cm  ,  an  order  of  magnitude  higher  than  the  current 
densities  of  the  junctions  studied  in  this  work.  A  further  increase  in 
current  density  (or  decrease  in  junction  area)  could  be  expected  to  come 
from  a  decrease  in  the  minimum  linewidth  by  the  use  of  either  deep  UV 
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FIG.  5.1  Input  impedance  of  a  bow-tie  antenna  as  a  function  of  bow 
angle  for  an  antenna  on  a  fused  quartz  substrate.  Solid  lines  are 
calculated  values,  crosses  are  values  measured  from  a  10  GHz  microwave 
modelling  experiment.  [After  Ref.  65.] 
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contact  lithography  (minimum  linewidth  Ml.2-0. 3  ym)  or  the  use  of 
electron  beam  exposure  of  a  two-layer  electron-beam  resist  system 
(minimum  linewidth  <0.1  ym).  A  group  at  Bell  Labs  has  fabricated 
superconducting  tunnel  junctions  with  areas  below  0.01  (ym)^  ^°,  a 
factor  of  roughly  50  smaller  in  area  than  the  smallest  junction  we  have 
studied.  Clearly  there  is  considerable  room  for  improvement  in  junction 
fabrication. 

A  third  area  where  the  experiment  can  be  improved  is  concerned  with 
the  question  of  whether  external  noise  gives  rise  to  the  somewhat  above 
bath  (~2  K)  noise  temperatures  which  were  shown  to  be  able  to  cause  the 
low  IcRn  products  of  the  high  resistance  junctions.  The  experiments 
were  not  carried  out  in  an  electromagnetically  shielded  room.  It  would 
be  beneficial  to  construct  an  electromagnetic  shield  around  the  optical 
dewar. 

Of  practical  significance,  this  work  has  shown  that  the  infe-red 
responsivity  of  these  junctions,  when  viewed  as  a  direct  detector, 
approaches  the  quantum  limit  of  one  electron  transferred  across  the 
barrier  for  each  photon  absorbed.  It  is  desirable  to  study  the  perfor¬ 
mance  of  SIS  and  SIN  mixers  and  direct  detectors  in  the  submillimeter 
wavelength  region.  At  this  time  at  least  two  groups  are  preparing  to 
study  SIS  mixers  operating  in  the  submillimeter  wavelength  region  with 
quantum  limited  performance,  as  described  below. 

At  Caltech,  an  SIS  (Pbln-Pb)  mixer  designed  for  250  GHz  has  a 
measured  noise  temperature  of  350  K  at  250  GHz  and  has  been  studied  up 
to  400  GHz  where  it  has  a  noise  temperature  of  -2000  K.^  The  350  K 


noise  temperature  is  still  considerably  higher  than  the  quantum  limit  of 
TN*fiti»^/e~l 2  K  for  250  GHz,  but  it  is  the  lowest  reported  mixer  noise 
temperature  at  this  frequency.  The  mixer  block  is  machined  from  a  single 
piece  of  copper.  It  uses  a  circular  waveguide  (0.97  mm  diameter)  fed  by  a 

corrugated  horn,  and  an  adjustable  backshort  for  tuning. The 
junctions  were  fabricated  on  0.1  nm  thick  quartz  substrates  which  fit 
into  a  groove  in  the  mixer  block  which  spanned  the  waveguide.  This 
scheme  could  be  extended  to  higher  frequencies  by  using  a  smaller  dia¬ 
meter  waveguide  and  so  on.  A  group  at  U.  C.  Berkeley  is  using  a 
different  approach  for  a  600  GHz  mixer  in  which  a  Nb-Si-Nb  SNAP  junction 
is  coupled  to  the  signal  and  local  oscillator  sources  through  a  long- 
wire  antenna  mounted  in  a  comer  cube.  The  long-wire  antenna  is  designed 
to  touch  a  gold  pad  on  top  or  near  the  tunnel  junction  to  couple  the 

radiation  to  the  junction  while  the  junction  reactance  is  tuned  out  with 
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an  open  circuited  transmission  line.  Both  approaches  have  good 
potential  for  yielding  quantum-limited  heterodyne  mixing  performance. 

In  sunmary,  we  have  performed  the  first  measurements  of  thin-film 
superconducting  tunnel  junctions  in  the  submillimeter  wavelength  region 
of  the  spectrum.  Our  data  agrees  well  with  the  fundamental  theory  of  the 
tunnel  junction.  Considerable  work  remains  to  be  done  to  extend  this 
work  to  still  higher  frequencies  and  to  develop  practical  SIS  mixers 
which  can  operate  in  this  frequency  range. 


APPENDIX  A 


DETAILED  SAMPLE  FABRICATION  PROCEDURES 


A.l  Introduction 

This  appendix  describes  in  detail  the  complex  procedures 
needed  to  make  small -area  superconducting  tunnel  junctions. 

The  first  section  provides  complete  details  of  the  photolithographic 
procedures  used  in  fabricating  the  samples.  The  second  section 
is  a  thorough  description  of  the  thin  film  evaporation  techniques, 
while  the  third  section  details  the  substrate  cleaning  procedures. 
The  Appendix  ends  with  formulas  for  the  etches  used  throughout  this 
work. 


A. 2  Photo! ithographic  Procedures 

A.  Photography 

1.  Prepare  black-tape  patterns  on  transparent 
acetate  sheets. 

—  250x  scale  for  +10  microscope  reduction. 

—  1250x  scale  for  +50  micrscope  reduction. 

—  Note:  a  +10  reduction  produces  -v5  ym  minimum 

linewidths,  while  +50  reduction  produces  1-2  ym 
minimum  linewidths.  No  microscope  reduction  (i.e.  +1  ) 
produces  20  ym  minimum  linewidths,  due  to  the 
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grain  size  of  the  film. 

-  The  black  tape  regions  on  the  pattern  become 
clear  on  the  film.  Clear  regions  on  the  emulsion 
will  expose  photoresist  producing  clear  regions 
on  a  substrate.  Therefore  black  regions  on  the 
film  will  determine  where  metal  is  in  the  sample 
because  the  metallization  is  done  using  a  lift-off 
technique. 

Details  of  photography. 

-  Film  type  -  Kodalith  Ortho  film  type  6556 
with  the  red  estar  base,  bulk  roll. 

-  f/22  aperture  on  camera. 

-  1-7  second  exposure  time. 

-  Developer  concentrations,  for  1600  ml  total; 

make  up  800  ml  Kodalith  Fineline  developer  A  using 
107.4  gn  of  powder,  make  up  800  ml  Kodalith 
Fineline  developer  B  from  53.7  gm  powder  by 
mixing  with  distilled  water.  Temperature 
of  solutions  should  be  20  0  C.  Mix  solutions 
only  just  before  film  development  since  the 
mixed  developer  is  good  only  for  1-2  hours. 

-  Shoot  film  63"  from  the  lightboard. 

-  Film  process. 

Develop  2  45 

Stop  0:10  -  Kodak  Indicator  Stop  Bath 


Fix  2:00  -  Kodak  Rapid  Fixer 
Rinse  10:00  in  running  water 
Dip  in  Photo-flo  and  hang  dry. 

—  Result--  *25  reduction  of  original  pattern  with 
black  and  white  reversed. 

—  Ultimate  resolution  of  film  ^5  ym  or  200  lines/mm. 

B.  Intermediate  Masking  Step 

1.  Use  glass  microscope  slides  cut  into  1"  x  1" 
pi eces. 

2.  Clean  with  standard  cleaning  procedure. 

o 

3.  Coat  with  1000  A  Cr,  using  Cr  plated  Tungsten 
rod  evaporation  source. 

4.  Coat  with  1.5  ym  AZ  1350J  photoresist,  spun  on 
at  5500  RPM,  baked  at  80°  C  for  20-25  min. 

5.  Contact  print  by  putting  1"  x  1"  substrate 

into  contact  with  emulsion.  The  emulsion  and  photoresist 
must  be  in  contact.  Sandwich  glass  and  emulsion 
betweeen  two  glass  plates.  Expose  for  90  seconds 
at  a  distance  of  2'  from  the  UV  lamp. 

6.  Develop  photoresist  in  high-resolution  make-up 
developer  —  5  parts  distilled  water  to  1  part 
developer  (AZ  351)  --  for  a  total  of  30  -  60  secs. 

There  should  be  no  photoresist  scum  on  the 
surface  of  the  Cr,  i.e.  it  should  look  quite  shiny. 

7.  Clean  in  photoresist  stripper  solution. 
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8.  Store  in  tightly  sealing  containers. 


Projection  steps. 

1.  Contact  print  pad  mask  onto  2"  x  2"  x  0.060" 

low  reflectivity  Cr,  master  mask  blank,  coated  with 
0.5  pm  AZ  1350B  photoresist. 

2.  Expose,  develop,  rinse,  and  blow  dry  with  pure  nitrogen. 

3.  Projection: 

—  Attach  light  source  to  projection 
attachment  of  microscope. 

—  Adjust  microscope  optics  so  light 
can  travel  thru  projection  attachment 
onto  substrate. 

—  Set  lamp  power  supply  to  10  V. 

—  Use  8x  objective. 

—  Expose  35-40  secs. 

—  Attachment  dial  12-14  mm  setting. 

--  Put  substrate  into  perfect  focus, 
then  lower  microscope  stage  slightly 
so  that  projected  pattern  is  slightly 
out  of  focus  to  attempt  to  focus 
the  UV  radiation  on  top  of  the 
photoresist. 

--  Always  use  red  filter  when  inspecting 
photoresist. 

4.  Develop  30-60  secs.,  rinse,  dry. 


5.  Inspect  in  microscope,  if  OK  go  on  to  the 

next  stage,  otherwise  strip  photoresist,  and  recoat 
mask  blank  with  0.5  ym  AZ  1350B  photoresist,  then 
go  back  to  step  3. 

6.  If  OK  etch  Cr. 

7.  Strip,  clean,  recoat  with  0.3  ym  AZ  1350B  photo¬ 
resist,  bake  20-25  min.,  let  it  sit  in  air  for  20  mi 
to  regain  lost  humidity. 

8.  *50  projection  thru  microscope 

—  40x  objective  on  microscope 
--  Lamp  on  5  V. 

--  Expose  for  10  secs. 

—  Focus  as  before. 

9.  Develop,  inspect,  and  etch  if  good. 

10.  Expose  large  dots,  using  8x  objective,  to  indicate 
proper  sample  orientation  during  evaporation 
sequence. 

11.  Etch  Cr,  clean  mask,  and  store  in  mask 
storage  box. 

Working  Copies 

1.  Using  cleaned  Master  Mask,  contact  print 
Master  onto  a  2"  x  2"  x  0.010"  flexible 

0.5  ym  thick  photoresist  coated  low  reflectivity 
Cr  mask  blank. 

2.  Expose  for  60  secs.,  develop,  rinse,  dry. 


3.  Inspect,  clean  if  OK,  and  this  mask  will 
become  a  working  copy,  otherwise  recoat 
and  repeat  steps  1  and  2. 

4.  Store  in  clean,  tightly  sealing  containers. 

Preparation  of  Sample  Substrates 

1.  1"  x  0.5"  x  0.010"  z-cut,  optically  polished, 
single  crystal  quartz  substrates. 

2.  Clean  using  standard  cleaning  procedure, 
described  in  Section  IV.A  below. 

3.  Coat  with  1.5  urn,  A2  1350J  photoresist,  bake 
at  80°  C  for  20-25  min. 

4.  Expose  for  2. 5-3.0  min.  at  a  distance  of  2' 

from  the  UV  lamp, 
o 

5.  Deposit  500  A  A1  from  Tungsten  boat  source, 

using  diffusion-pumped  evaporator  at  a  base 
-6 

pressure  of  3  x  10  torr,  with  liquid 
nitrogen  trap  filled. 

6.  Coat  with  1.0  ym  AZ  1370  photoresist,  bake 
at  70°  C  for  20  min, 

7.  Substrates  are  ready  now  for  contact  printing. 

8.  Store  in  tightly  sealing  sample  storage  boxes. 

Contact  Printing  Procedure 
1.  Insert  sample  into  depression  in  plexiglass 

contact  printing  sample  holder  so  that  substrate 


surface  is  nearly  at  the  same  height  as  the 
surrounding  plastic.  See  Figure  A.l. 

2.  Position  flexible  working  mask  over  substrate, 

Cr  side  down,  center  pattern  over  substrate. 

3.  Open  vacuum  valve  so  substrate  and  mask 
are  tightly  held  together.  Check  for  good 
contact--! ook  for  dirt  particles  lifting 
mask  away  from  the  substrate. 

4.  Expose  for  90  seconds  at  2'  from  the  lamp. 

5.  Close  the  vacuum  valve,  remove  the  substrate  from 
the  holder,  repeat  for  other  prepared  substrates. 

Developing  Procedure 

1.  Develop  the  top  layer  until  the  photoresist  is 
entirely  removed  and  the  A1  layer  looks  shiny 
in  such  areas,  rinse,  and  dry. 

2.  Inspect  to  make  sure  the  photoresist  is  completely 
removed.  Optional. 

3.  Etch  A1  layer  until  A1  is  entirely  removed 
from  the  places  uncovered  by  step  1,  rinse, 
and  dry. 

4.  Develop  the  bottom  layer  until  a  few  seconds  after 
the  quartz  substrate  looks  shiny,  i.e. 

in  areas  that  only  had  the  bottom  layer  of 
photoresist  remaining,  watch  the  photoresist 
slowly  develop  away  until  it  appears  that 


Plexiglass  Contact  Printing  Sample  Holder 
(not  to  scale) 


1  mil  between 


FIG.  A.l  Illustration  of  cross-section  of  plexiglass  contact  printing 
sample  holder.  Substrate  sits  in  slot  in  holder  with  substrate  surface 
lu  u  025  inches  below  the  top  of  the  holder.  A  vacuum  is  applied  so 
that  the  mask  is  pulled  into  hard  contact  with  the  substrate. 


the  substrate  is  clean,  hold  the  substrate 
in  the  developer  a  few  seconds  longer,  rinse 
quickly,  dry. 

Etch  the  A1  layer  again.  This  is  done  to  make 
it  possible  to  see  the  undercut  when  the 
sample  is  placed  under  the  microscope  for 
inspection. 

Inspect,  if  sample  needs  more  of  an  undercut 
then  repeat  steps  4  and  5. 

If  OK,  store  the  sample  in  clean  room 
storage  boxes  until  needed  for  the 
evaporation  sequence. 

Note: 

--  Developer  is  the  high  resolution 
make-up  mentioned  above. 

--  It  is  important  that  _aV[  the 
photoresist  is  removed  from  the 
surface  of  the  quartz. 

—  Usually  9  out  of  10  substrates  could  be  used 


in  the  evaporation  sequence. 
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A. 3  Evaporation  Procedures 

A.  Use  of  Cryopunp  evaporator  --  before  evaporation. 

1.  Bake  out  molecular  sieve  the  night  before 
using  the  evaporator. 

2.  Base  pressure  should  be  2-4  x  10"7  torr 
before  beginning. 

3.  Check  glow  discharge  apparatus,  look  for  shorts. 

4.  Center  of  glow  discharge  wire  should  be 
directly  under  the  sample  holder,  but 
should  be  2-3“  wider  than  the  sample  holder 
all  the  way  around.  The  wire  should  be 

at  a  negative  potential  with  respect  to  ground. 

The  base  plate  of  the  evaporator  should  be  grounded. 

5.  The  digital  thickness  monitor  crystal  should 
be  tested  for  thickness  of  accumulated 
metal,  replaced  if  necessary. 

6.  Check  metals  --  make  sure  there  is  enough. 

a.  Use  only  etched  high-purity, 

99.999-99.9999  %,  Pb.  The  Pb  is 
etched  in  an  acetic  aci d/peroxide/ 
distilled  water  solution  detailed  in 
etch  section  of  this  Appendix.  Dip 
Pb  in  etch,  rinse  quickly  in  distilled 
water,  then  in  methanol.  Keep  in 
methanol  until  ready  to  put  in  evaporator. 


Dry  with  clean  room  wipers  or  dry  nitrogen 
gas. 

b.  Sn  does  not  need  to  be  etched,  but  should 
also  of  very  high  purity. 

c.  Ge  does  not  need  to  be  etched  either, 
also  high  purity. 

7.  Put  clean  metals  in  clean  tungsten  boat  sources. 

Boat  sources  should  be  cleaned  by  ultrasonic 
agitation  in  Methanol  before  inserting  into  evaporator. 

8.  Sample  holder  should  be  reasonably  clean;  if  too 
much  metal  has  been  evaporated  onto  it,  it  should  be 
disassembled,  etched  in  copper  etch,  rinsed 

with  distilled  water,  and  methanol,  and  blown 
dry  with  nitrogen  gas. 


Evaporation 

1.  Put  sample  in  sample  holder,  screw  onto 
bottom  of  liquid  nitrogen  pot,  rotate  into 
position  for  Sn  evaporation. 

2.  Rough  pump  evaporator  to  30-50  mtorr. 

3.  Close  rough  pump,  begin  pumping  with  cryopump. 

4.  In  xl  hour  pressure  should  be  in  mid  10"7 
range. 

5.  When  pressure  is  3-4  x  10’7,  begin  liquid 
nitrogen  flow  to  the  liquid  nitrogen  pot. 


In  1/2  hour,  the  pot  should  be  nearly 
at  liquid  nitrogen  temperature. 


6.  Keep  liquid  nitrogen  flowing  while 
depositing  metal  films. 

o 

7.  Evaporate  Sn  at  100-200  A/sec  for 

o 

a  total  thickness  of  1000-1500  A  when  corrected 
for  the  angle  of  the  substrate  with  respect  to 
the  source.  Base  pressure  of  the  evaporator 
should  be  2-4  x  10”  7  torr.  The  boat  should  be 
slowly  heated  until  the  pressure  begins  to  rise 
in  the  chamber;  at  this  time  the  shutter  should 
be  opened  and  the  rate  should  be  adjusted  if 
necessary. 

8.  Rotate  substrate  to  proper  angle  for  Ge  evaporation. 

o 

9.  Evaporate  Ge,  only  need  300-500  A  when  corrected 
for  the  angle  of  the  evaporation. 

10.  Rotate  the  substrate  to  the  proper  angle 
for  the  Pb  evaporation. 

11.  Close  gate  valve,  bleed  in  pure  oxygen,  until 
the  pressure  is  30  mtorr.  Wait  1  minute. 

12.  Pump  out  chamber  with  cryopump  for  2-3  minutes, 
until  at  the  proper  base  pressure. 

13.  Close  gate  valve,  bleed  in  30  mtorr  of  oxygen, 


wait  1  minute,  turn  on  glow  discharge.  Voltage 
should  be  1.1  kV,  current  12-15  mA,  and  the 
color  should  be  gray-green.  Continue  discharge 


for  45-75  secs.,  depending  on  junction  resistance 
desired.  Pressure  should  drop  to  about  20  mtorr, 
and  the  current  should  fall  2-4  mA,  as  the  oxygen 
is  consumed. 

14.  Open  gate  valve,  chamber  pressure  should  go  down 

rather  quickly  to  base  pressure,  2-4  x  10”7  torr. 

o 

Evaporate  Pb  100-200  A/sec  to  make  a  total  thickness 
o 

of  3000-4000  A  when  corrected  for  the  evaporation 
angle. 

15.  Turn  off  liquid  nitrogen  flow,  blow  nitrogen  gas 
thru  the  liquid  nitrogen  lines  for  about  1/2  hour 
to  warm  up  sample  holder. 

16.  Wait  until  sample  holder  has  warmed  up. 

17.  Take  sample  out  of  holder,  put  holder  back 
in  the  evaporator,  and  pump  down  to  a  high 
vacuum.  Evaporator  should  always  be  left 
under  high  vacuum  condition,  except  when 
samples  are  inserted  or  removed  from  the 
evaporator. 

C .  Sampl e 

1.  Photoresist  should  be  removed  by  soaking 

the  sample  in  acetone  for  1/2  hour,  followed  by 
rinsing  it  with  acetone,  then  methanol,  and  blowing 
it  dry  with  nitrogen. 


A. 4  Cleaning  Procedures 

A.  Simple  Procedure 

If  the  substrates  are  basically  clean,  i.e.  as  they  come 
from  the  supplier,  then  the  simple  cleaning  procedure  is 
all  that  is  called  for: 

1.  Insert  substrates  into  Teflon  holder. 

2.  Ultrasonically  agitate  in  Trichloroethylene  (TCE) 
for  10  mins. 

3.  Ultrasonically  agitate  in  Acetone  for  10  mins. 

4.  Ultrasonically  agitate  in  Methanol  for  10  mins. 

5.  Blow  dry  with  nitrogen  gas. 

B.  Procedure  for  Dirty,  But  Not  Coated  Substrates 

Some  times  it  is  helpful  initially  soak  the  samples 
in  a  t2%  solution  of  Micro  cleaner  and  distilled  water 
heated  to  near  boiling  or  boiling,  for  1/2  hour, 
ultrasonically  agitated  for  10  mins.,  and  blown  dry  with 
nitrogen  gas.  Steps  1-5  above  are  then  followed. 

C.  Procedure  for  Coated  Substrates 

If  the  samples  had  been  used  before,  there  are  two 


possibilities: 


1.  If  covered  by  photoresist,  the  substrates 
should  be  soaked  in  photoresist  remover, 
which  is  nearly  boiling,  for  20  mins., 
rinsed  with  distilled  water,  and  blown  dry. 

2.  If  the  substrates  have  metal  films  on  them 
they  should  be  soaked  in  Aqua  Regia  for 

1/2  hour  until  the  films  are  removed,  rinsed, 
and  dried. 

3.  After  completing  Cl  or  C2  above  then  steps  A1-A5 
above  should  be  followed. 

The  photoresist  stripper  is  1  part  Shipley  Photoresist 
Remover  1112A  to  2  parts  distilled  water. 


A. 5  Etches 


In  the  formulas  which  follow  all  proportions  are  by  volume 
unless  otherwise  noted. 

Aluminum. 


84 

H3°4 

Phosphoric  Acid 

6 

h2o 

Distilled  Water 

5 

ch3cooh 

Acetic  Acid 

5 

hno3 

Nitric  Acid 

Aqua  Regia: 

1 

hno3 

Nitric  Acid 

3 

HC1 

Hydrochloric  Acid 

3 

h2o 

Distilled  Water 

Copper: 

1 

hno3 

Nitric  Acid 

3 

h2o 

Distilled  Water 

Chrome: 

165  gm 

NH4CeN03 

Ceric  Ammonium  Ni 

43  ml 

hcio4 

Perchloric  Acid  ( 

Ho0 


Ceric  Ammonium  Nitrate 
Perchloric  Acid  (70%) 

Distilled  Water  to  make  1  liter  solution 


APPENDIX  B 


OTHER  DATA  AND  ANALYSIS 


B.l  Introduction 

In  this  Appendix  we  will  discuss  experimental  results  which  were 
obtained  either  prior  to  (and  in  preparation  for)  the  experiments  which 
comprise  the  bulk  of  this  work,  or  which  were  completed  after  the  major 
part  of  the  work  was  completed.  Section  B.2  discusses  experiments  per- 

O  O 

formed  on  very  low  resistance  (Rn“0.1  a),  larger-area  (A-5-20xl0  cm  ) 

junctions  fabricated  at  Bell  Laboratories  in  Holmdel,  NJ.  Results  of 

far-infrared  (FIR)  measurements  on  edge-aligned  tunnel  junctions 
-8  2 

(A<10  cm  ),  with  point-contact  or  microbridge  type  I-V  curves,  are 
reported  in  Section  B.3.  The  Appendix  ends  with  Section  B.4  which 
describes  results  from  experiments  with  small-area  SNAP47  junctions 
fabricated  at  Sperry  Research  Laboratories. 

B.2  Low  Resistance  Junctions  Fabricated  at  Bell  Labs 

When  the  construction  of  the  cryopumped  evaporation  system  at 
Harvard  was  nearing  completion,  the  author  of  this  report  made  a  set  of 
tunnel  junctions  at  Bell  Labs  using  tunnel  junction  evaporation/oxidation 
facilities  available  there.  Photoresist  stencils  for  this  work  were 
made,  however,  at  Harvard. 

The  Bell  Labs  (Holmdel)  evaporation  system  was  constructed  out  of  a 
bakeable  stainless  steel  chamber  pumped  by  a  mechanical  roughing  pump 


from  atmospheric  pressure  to  the  millitorr  range  and  an  ion  pump  at  high 

vacuum.  Contamination  from  the  mechanical  roughing  pump  was  prevented 

from  migrating  to  the  chamber  by  the  presence  of  a  molecular  sieve  trap 

and  a  liquid  nitrogen  cold  trap  in  the  roughing  pump  lines.  A  mass 

spectrometer  was  attached  to  the  chamber  to  monitor  the  gases  in  the 

chamber  and  to  detect  contaminants.  To  minimize  the  adsorption  to  the 

chamber  of  water  vapor  and  other  contaminants  from  the  atmosphere,  the 

chamber  was  always  kept  under  high  vacuum,  except  when  samples  were 

transferred  into  or  from  the  chamber.  A  liquid  nitrogen  baffle  was 
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mounted  in  the  chamber  to  help  pump  the  water  vapor. 

The  stencils  used  were  fabricated  using  the  commercial  mask  des¬ 
cribed  in  Section  3.2.1  of  Chapter  3,  and  the  photoresist  exposure/ 
development  sequence  summarized  in  Fig.  3.3.  These  Sn-Pb  junctions  were 

made  using  the  simple  rverlap  technique  (Fig.  3.5)  with  areas  ranging 
-8  2 

from  -5-20x10  cm  ,  and  resistances  from  -0.10-2.0  ft.  The  Sn  film  was 
typically  -1000  A  thick,  while  the  Pb  film  was  somewhat  thicker,  -2000  A 
These  films  were  evaporated  from  metal  shot  in  resistively  heated  boat 
sources  onto  room  temperature  sapphire  substrates  onto  which  the  photo¬ 
resist  stencils  had  been  previously  prepared.  The  Sn  film  was  oxidized 

by  glow  discharge  in  pure  oxygen  at  7  millitorr  for  times  ranging  from 

4  5  2 

7-60  seconds,  producing  junctions  with  current  densities  -10  -10  A/cm  , 
and  resistances  and  areas  noted  above. 

The  I-V  characteristics  at  selected  temperatures  of  a  typical 
sample  prepared  in  the  above  manner  is  shown  in  Fig.  B.l.  This 


6 
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particular  junction  has  an  area  of  -12x10  cm  ,  Rn~0.35  ft,  and  Ic ( 1 .4 
K)=2. 3  mA.  The  capacitance,  C,  is  estimated  to  be  00.32  pF, 
ec=(2eIcR^C/fi)=  0.30  at  1.4  K,  and  the  toJ=(6c)i'2/Rn05.2xl012(sec'1 ) . 

6  and  w,  are  computed  using  the  actual  measured  I  ,  not  from  an 

C  u  C 

inferred  zero  temperature  Ic.  See  the  discussion  in  Chapter  4,  Section 
4.2  for  more  details  on  this  question. 

Fig.  B.l  shows  I-V  curves  at  six  temperatures,  T=4.2,  3.61,  3.31, 
3.07,  and  1.4  K.  We  observe  that  at  3.61  K,  the  junction  switches  at  I 
from  the  zero  voltage  state  to  the  structure  at  ^pb*Asn*  and  then  at  a 
higher  current  level,  to  the  structure  at  Apt+ASn.  At  higher  currents 
still,  there  is  a  decrease  in  voltage  as  a  function  of  current,  until  a 
point  is  reached  where  the  voltage  again  begins  to  increase  with 
increasing  current,  i.e.  the  normal  resistance  is  being  approached. 

This  trend  is  continued  at  lower  temperatures,  T<3  K,  shown  in  the 
second  set  of  I-V  curves.  Note  that  the  second  set  is  plotted  with  a 
scale  that  is  1/2  as  sensitive  as  the  first  set  of  I-V  curves.  Similar 
data  was  obtained  by  Yeh  and  Langenberg"116  on  Sn-SnO-Sn  junctions.  Fig 
B.2  shows  an  I-V  curve  from  Yeh  and  Langenberg  at  1.3  K,  displaying 
qualitatively  the  same  type  of  gap  depression  observed  in  our  data.  On 
the  I-V  curve  at  1.4  K,  the  lowest  temperature  reached  in  this  experi¬ 
ment,  the  I-V  curve  of  a  high  resistance  (Rn~1600  ft),  low  current- 
density  junction  is  plotted  in  dashed  lines  for  comparison.  The  I-V 
curve  from  the  high  resistance  junction  is  scaled  by  multiplying  the 
measured  current  of  the  high  resistance  at  a  particular  dc  voltage  t. 
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FIG.  B.l  A  typical  high  current-density  sample  showing  depressed  gap 
structure  from  quasi  particle  injection  at  6  temperatures.  At  1.4  K  a 
scaled  I-V  curve  from  a  high-resistance  junction  is  drawn  in  dashed 
lines  for  comparison  (see  text).  The  three  upper  curves  are  drawn  with 
a  scale  twice  as  sensitive  as  the  lower  three  curves. 
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FIG.  B.2  An  I-V  curve  for  a  Sn-SnO-Sn  junction  (Ref.  116)  showing  the 
same  sort  of  gap  depression  phenomena  observed  in  this  work  on  Sn-SnO-Pb 
junctions. 


the  ratio  of  the  resistances  at  voltages  just  above  the  gap  region,  i.e., 
*scaled“*high  rx(^00  fl/0.29  a).  As  can  be 'clearly  seen  the  voltage  at 
the  sum  of  the  gaps  is  rather  low,  -1.75  mV,  near  the  point  at  which  the 
junction  switches  from  the  non-zero  voltage  state  to  the  zero  voltage 
state,  and  decreases  to  a  minimum  of  -1.68  mV  before  increasing  again. 
This  is  in  contrast  to  the  behavior  seen  on  the  high  resistance  junction 
which  has  approximately  the  same  voltage  at  the  point  the  low  resistance 
junction  has  gap  structure  at  -1.75  mV,  point  A  in  the  Figure.  As  the 
current  is  Increased  the  voltage  Is  Increased  slightly  to  -1.90  mV,  at 
point  B,  giving  a  maximum  difference  of  -280  yV.  At  higher  current 
levels  the  curves  match  onto  each  other.  The  voltage  difference  between 
the  two  curves,  6Vg,  is  plotted  as  a  function  of  current  in  Fig.  B.3. 

The  behavior  described  above  has  been  studied  extensively  by  Yeh 

and  Langenberg.^  The  difference,  -6Vg,  is  interpreted  by  them  as 

coming  from  gap  depression  due  to  the  heavy  quasi particle  Injection  in 

the  low-resistance  junction  at  the  gap  edge,  and  quantitatively  described 

in  terms  of  phenomenological  equations  developed  by  Rothwarf  and 

Taylor.  ^  These  equations  relate  the  excess  quasiparticle  density 

above  the  gap  compared  to  the  number  of  thermally  excited  quasi particles 

as  a  function  of  the  current  Injected.  The  gap  depression  is  related  to 

*  IIS 

the  excess  number  of  quasi particles  by  the  T  model  of  Parker  and 

*  119 

the  y  model  of  Owen  and  Scalaplno  ,  which  give  equivalent  results  for 
the  low  temperature,  1.4  K,  of  the  sixth  I-V.  Fig.  B.3  (for  l£3mA)  looks 
rather  like  Fig.  5  of  Ref.  B.4  reproduced  as  Fig.  B.4  here  in  which  -6Vg 
is  as  we  defined  it  above,  and  JG  is  proportional  to  the  current.  The 


I(mA) 


B.3  Voltage  difference,  V  ,  between  the  dashed  curve  and  solid 
in  Fig.  B.l  at  T»1.4  K  versus  the  current  in  the  junction. 
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FIG.  B.4  Voltage  difference,-6V  ,  between  gap  structure  of  low 
resistance  and  high  resistance  junctions,  computed  in  the  same  way  as 
Fig.  B.3  for  the  Sn-SnO-Sn  junctions  of  Ref.  116  as  a  function  of  JG 
which  Is  proportional  to  the  current  through  the  junction.  There  is 
clearly  a  great  deal  of  similarity  between  this  curve  and  Fig.  B.3. 


solid  lines  In  Fig.  B.4  are  fits  from  the  non-equilibrium  model  described 
above.  Unlike  Yeh  and  Langenberg,  we  will  not  carry  the  analysis  any 
further. 

Qualitatively  their  I-V  curves  and  ours  are  similar,  although  the 
gap  depression  seen  by  them  Is  roughly  a  factor  of  ~5  less  than  seen  In 

our  data.  This  might  be  due  to  the  higher  current-density  of  our  junc- 

4  2 

tions,  -10  A/cm  compared  to  the  current-density  of  their  junctions, 

-102  A/cm2. 

The  analysis  above  was  not  pursued  further  because  this  sort  of 
behavior,  while  Interesting,  was  not  expected  in  the  high  resistance, 
smaller-area  junctions  which  were  to  be  used  In  the  FIR  work.  Non¬ 
equilibrium  behavior,  such  as  seen  here,  may  become  Important  in  very 

-9  2 

high  current-density,  extremely  small-area  junctions  (A<10  cm  ) 
junctions  which  may  be  studied  In  the  FIR  In  the  future. 

B. 3  Point- Contact  Edge-Aligned  Junctions 

The  successes  with  the  resist-aligned  junction  technique  reported 

in  Chapter  IV  motivated  us  to  attempt  to  use  the  edge-aligned  junction 

technique  to  produce  junctions  with  areas  smaller  than  the  smallest 

-8  2 

junctions  fabricated  thus  far,  FwOAxlO  cm  .  The  edge-aligned  junc- 

-8  2 

tions  held  the  promise  of  junctions  with  areas,  A~0.2xl0  cm  ,  a 
factor  of  two  or  more  smaller  in  area,  which  would  significantly 
Increase  the  power  coupled  at  the  frequencies  used  thus  far,  and  would 
enable  us  to  more  easily  couple  to  higher  frequencies,  'ftu>L>2Upb+ASn)/e. 


We  fabricated  a  series  of  junctions  using  the  edge-aligned  junction 
technique  discussed  In  Section  3.2.4.  Of  six  junctions  fabricated,  only 
one  was  a  proper  tunnel  junction,  displaying  the  kind  of  I-V  curve 
expected  from  tunneling  theory.  That  junction  had  a  resistance  of  -1600 
Q  and  was  reported  In  Section  4.2.  Parameters  for  the  remaining  junc¬ 
tions  are  listed  In  Table  B.l.  These  junctions  had  resistances  ranging 
from  -0.55  ft  to  -75  ft.  All  but  the  0.55  n  junction  had  some  sort  of  gap 
structure  at  (Apb+Ajn)/e,  and  all  had  excess  currents  *  of  -0.8  Ic. 

In  simple  terms  an  "excess  current"  exists  if  the  linear  high-voltage 
I-V  extrapolates  to  a  non-zero  current  at  zero  voltage,  instead  of 
extrapolating  to  zero  current  at  zero  voltage  like  a  resistor  or  a  tunnel 

junction  would.  This  behavior  Is  characteristic  of  a  metallic  weak  link 
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like  the  point  contacts  of  Weltz  et  al.  The  IcRn  products  were  mostly 
much  higher  than  IcRn*l-3  mV,  expected  from  tunneling  theory.  The 
average  IcRn  was  -1.7  mV. 

Figure  B.5  displays  the  I-V  characteristics  of  the  lowest  resistance 
sample  Rn~0.55  ft,  at  two  power  levels,  2a*0,  and  the  highest  power 
coupled,  2o*0.8,  where  2a=2eVL/fiu^,  as  usual.  This  I-V  curve  is  linear, 
with  little.  If  any,  gap  structure.  The  low  sample  resistance  causes  a 
large  Impedance  mismatch  between  the  sample  and  antenna,  which  explains 
why  the  coupled  power  is  so  low.  However,  the  critical  current  is 
depressed  to  about  75%  of  its  zero  power  value  and  a  clear  Josephson 
step  Is  seen  at  1.25  mV. 

I-V  curves  for  a  series  of  laser  power  levels  for  a  junction  with 
I  R  product  close  to  the  theoretical  value  of  1.3  mV,  are  shown  in  Fig. 


Sample  B-7 
Rjj  **»  0.55  £1 
Ic~  3.4  mA 
T  *  1.4K 


2a«0 


2a  «  0.8 


7  mA 


FIG.  B.5  I-V  curves  of  a  0.55  ft  edge-aligned  point  contact  type  sample 
at  2a*0,  and  2a«0.8,  the  maximum  power  coupled  to  the  junction,  where 
2o*2eV.  /ftw.  .  The  curve  at  2o*0.8  clearly  shows  a  step  at  1.25  mV, 
label ifed  n*l. 


B.6.  This  junction  has  strong  gap  structure  at  Vg*(Apb+A$n)/e*  and  has 
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structure  at  lower  voltages,  possibly  subharmonic  gap  structure. 

More  power  Is  coupled  Into  this  low  g  ,  non-hysteretic  junction,  and  a 
total  of  five  steps  are  seen  at  the  highest  power  coupled,  2a=6.7.  A 
new  feature  on  these  I-V  curves,  not  seen  previously  with  point  contacts 
in  the  FIR,  is  a  half-harmonic  step  at  -0.61  mV,  labelled  n=l/2  in  the 
Figure.  With  the  exception  of  the  half-harmonic  step  these  I-V  curves 
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are  remarkably  similar  to  the  point  contact  I-V  curves  of  Weitz  et  al. 

The  large  amount  of  power  coupled  to  the  junction  makes  it  possible 
to  make  further  comparisons  with  point  contact  data  by  plotting  the  step 
half-widths  versus  the  normalized  ac  voltage,  2a,  as  had  been  done  pre¬ 
viously  In  the  analysis  in  Section  4.3.  Fig.  B.7  shows  the  step  widths 
for  this  sample  compared  with  the  Bessel  function,  | J  | ,  and  Werthamer 
function  predictions,  |Wj.  In  order  to  obtain  good  fits  with  the  step 
width  data,  these  theoretical  step  width  dependences  had  to  be  multiplied 
by  scaling  factors  less  than  unity  which  decreased  in  magnitude  with 
increasing  voltage.  The  scaling  factors  decreased  from  -0.4-0. 6  for  the 

zeroth  and  first  steps  to  -0.1 -0.2  for  the  fourth  and  fifth  steps.  Weitz 
12 

et  al.  observed  somewhat  larger  scaling  factors  which  varied  from 
-^0. 5-1.0  for  the  zeroth  and  first  steps  and  decreased  to  -0.2  for  the 
fifth  step.  These  results  can  be  contrasted  with  the  low  resistance, 
R^O  ft,  tunnel  junctions  of  Chapter  IV  which  showed  no  need  for  any 
scaling  factors,  or  the  higher  resistance  junctions  in  which  the  scaling 
factors  were  greater  than  unity  and  had  little  dependence  on  voltage. 


We  leave  the  n*l/2  step  unfitted  because  the  Werthamer  theory  and 
the  simple  voltage-biased  RSJ  model  do  not  provide  for  a  subharmonic 
step.  The  presence  of  this  step  may  be  due  to  the  much  higher  capaci¬ 
tance  of  the  tunnel  junction  as  compared  to  the  point  contact.  It  may 
also  be  due  to  the  impedance  seen  by  the  junction  at  the  laser  frequency 

which  may  be  much  higher  than  the  junction  impedance,  giving  a  bias 
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which  may  approach  a  current  bias.  Braiman  et  al.  have  predicted 
the  existence  of  these  subharmonic  steps  under  similar  conditions  but 
did  not  provide  any  analytic  expressions  with  which  to  fit  data. 

An  approximate  value  of  (2<*)|nax  for  each  junction  can  be  determined 
from  the  power  dependence  of  the  step  widths.  These  values  of  (2a) 
are  plotted  in  Fig.  B.8,  which  shows  an  approximately  linear  relation¬ 
ship  between  (2al„  and  R  ,  as  indicated  by  the  solid  line  in  the 
max  n 

Figure.  This  relationship  implies  that  the  source  impedance  at  the  laser 
frequency  must  be  much  larger  than  the  junction  impedance,  on  the  order 
of  200  n. 

The  data  described  above  are  interesting,  and  further  analysis  might 

be  beneficial,  but  we  did  not  pursue  it  any  further  because  we  were 

trying  to  get  away  from  the  problems  of  point  contacts,  rather  than 

inventing  a  planar- fabricated  point  contact.  A  planar  point  contact 

would  be  less  well -cooled  than  the  three  dimensional  structure  of  a 

regular  point  contact,  and  would  be  expected  to  have  a  lower  high- 
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frequency  limit.  *  The  reason  why  the  Sn-Pb  edge-aligned  junctions 
were  seldom  real  junctions  must  have  something  to  do  with  the  soft 
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FIG.  B.8  Maximum  value  of  the  normalized  laser  induced  voltage  in  the 
junction,  2a*2eV./f>w.  ,  versus  the  normal  resistance,  R  .  The  solid 
curve  is  drawn  to  show  the  approximate  linear  dependence  of  2a  on  normal 
resi stance ;  R„  (see  text). 
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materials  used.  Howard  et  al.  reported  resist-aligned  junctions  of 
Sn-Pb,  but  the  edge-aligned  junctions  were  of  Pbln-Pb.  Further  work  is 
clearly  needed  to  fabricate  smaller  area  junction  which  have  low- 
hysteresis  tunneling  I-V  curves. 

B.4  SNAP  Junctions 

The  soft-metal  junctions  fabricated  at  Harvard  and  used  in  all  the 
work  reported  thus  far  are  unable  to  cycle  between  room  temperature  and 
liquid  He  temperature  more  than  a  few  times.  Technology  based  on 
refractory  metals  such  as  Nb  provide  the  best  hope  for  the  robust, 
thermally-cyclable  type  of  junctions  needed  for  practical  applications. 
In  this  vein  we  have  collaborated  with  Sperry  Research  Labs  in  order  to 
test  the  suitability  of  Nb  junctions  for  high  frequency  (FIR)  applica¬ 
tions,  such  as  for  SIS  heterodyne  mixers. 

The  results  at  604  GHz  which  will  be  described  below  are  for  junc¬ 
tions  fabricated  by  the  selective  niobium  anodization  process  (SNAP)47 
developed  at  Sperry  Research  Labs.  The  SNAP  process  begins  with  a  clean 
polished  Si  wafer,  upon  which  are  sputtered  successive  layers  of  Nb, 
amorphous  Si  (a-Si),  and  Nb.  The  wafer  now  is  essentially  one  gigantic 
tunnel  junction  which  can  be  stored  for  several  months  without  any 
deterioration.  The  wafer  is  then  coated  with  photoresist,  which  is 
exposed  and  developed,  leaving  areas  uncovered  which  are  to  be  anodized. 
Areas  which  are  covered  by  the  photoresist  become  the  tunnel  junctions 
which  will  be  either  tested  singly,  or  through  later  wiring  steps  will 
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be  combined  to  form  logic  circuits.  The  anodization  process  proceeds 
until  the  entire  layer  not  covered  by  the  photoresist  is  converted  to 
niobium  oxide.  Next,  a  contact  is  made  to  the  non-anodized  niobium 
upper  electrode  in  which  the  anodized  Nb  Isolates  the  upper  and  lower 
contacts.  A  second  insulating  layer  can  also  be  deposited  before  the 
upper  Nb  contact  is  made. 

Fig.  B.9  displays  a  SNAP  junction  at  the  center  of  a  ~180  urn  long 

dipole  antenna  for  coupling  to  the  FIR  laser  source.  Fig.  B.9(a)  shows 

one-half  of  a  2  mm  by  4  mm  long  chip  upon  which  the  SNAP  junction,  its 

antenna,  and  dc  leads  are  fabricated.  The  antenna  is  the  top  of  the  "T" 

in  the  top  part  of  the  Figure.  At  the  bottom  of  the  Figure  pads  can  be 

seen  to  which  gold  wires  have  been  attached,  which  connect  to  larger 

pads  which  are  connected  to  wires  in  the  dewar  insert  (described  in 

Section  3.4.1).  Fig.  B.9(b)  is  a  close-up  of  the  antenna  of  10  ym 

linewidth.  To  the  left  are  the  dc  leads.  At  the  center  of  the  antenna 

Is  the  junction,  shown  from  two  views  in  Figs.  B.9(c)  and  (d).  Fig. 

B.9(c)  Is  taken  with  the  SEM  stage  tilted,  while  Fig.  B.9(d)  shows  a 

view  with  the  SEM  stage  in  the  normal  position.  The  small  square  at  the 

center  of  the  region  is  the  tunnel  junction.  Outside  the  tunnel  junction 

there  is  a  larger  square  (of  *3x3ym  area)  which  has  Nb,  Si,  and 

anodized  Nb,  while  outside  this  square,  in  addition  to  the  Nb,  Si,  and 

anodized  Nb,  there  Is  a  layer  of  Si 0.  Therefore,  looking  from  the  inside 

of  the  smallest  square  and  going  outward,  the  upper  Nb  layer  (contact 

layer)  is  in  contact  with  Nb,  then  NbgOg,  then  SiO.  The  tunnel  junction 

2 

has  an  area  of  ^1  ym  .  The  specific  capacitance  of  the  a-Si  barrier  has 


FIG.  B.9  SEM  micrographs  of  a  SNAP  junction,  (a)  Low  magnification 
view  of  sample,  displaying  dc  leads— vertical  parallel  lines  in  the 
Figure— and  pads  to  which  gold  wires  are  bonded,  (b)  Close-up  of 
antenna  region  of  (a),  showing  10  yin  linewidth  antenna,  tl80  um  long. 
The  junction  is  at  the  center  of  the  antenna,  dc  leads  are  to  the  left 
of  the  center,  (c)  Close-up  of  junction  region  at  the  center  of  (b) 
taken  with  SEM  stage  tilted.  The  junction  is  the  smallest  area  at  the 
center  of  the  two  larger  squares,  (d)  Close-up  of  junction  region  but 
without  a  tilt  on  the  SEM  stage,  the  junction  is  the  -v.1  um  by<\,l  ymregi 
at  the  center  of  the  Figure  (see  text). 


FIG.  B.10  I-V  curves  of  the  SNAP  junction  in  Fig.  B.9  irradiated 
with  604  GHz  ( x*496  ytn)  laser  radiation,  at  normalized  laser  induced 
voltages,  2a»0,  and  the  highest  voltage  coupled,  2aa1.7.  We  see  three 
Josephson  steps  at  V  «n(fi<»>,  /2e)  where  n  is  an  integer,  and  a  single 
photon-assisted  tunneling  step  at  V +(1W  /e)  where  V  =2amk,  as  labeled 
in  the  Figure.  9  L  9  ND 


been  measured  and  is  2.5  pF/ym  ,  giving  the  junction  a  capacitance  of 
~2.5  pF. 

The  junction  shown  In  Fig.  B.9  has  a  resistance  of  «7  ft,  and  a 
critical  current,  Ic,  of  215  yA  at  1.4  K.  Fig.  B.10  displays  I-V  curves 
at  two  power  levels  of  the  496  ym  laser  line:  2a*0,  and  at  the  highest 
power  reached,  2o«1.7.  When  2a*0,  we  clearly  see  the  expected  gap  struc¬ 
ture  at  Vgs2ANb*2.9  mV.  The  leakage  current  below  the  gap  Is  much 
larger  than  seen  with  the  Sn-Pb  junctions  used  previously.  At  2a*1.7, 
we  observe  that  the  critical  current  Is  depressed  to  about  35%  of  Its 
zero-laser-power  value.  Two  steps  are  easily  seen,  labelled  n*l  and 
n«2  In  the  Figure,  while  a  third  small  step  Is  barely  visible,  labelled 
by  n«3.  A  small,  photon-assisted  tunneling  step  Is  seen  at 
Vg+(Tui)L/e)=i5.4  mV.  This  step  may  perhaps  be  so  small  because  of  the 
large  leakage  current  seen  on  this  particular  sample. 

The  variation  of  the  step  widths  with  power  has  been  measured  and 
Is  shown  In  Fig.  B.ll.  None  of  the  steps  fit  the  expected  Bessel  func¬ 
tion,  1 Jn | ,  or  Werthamer  function,  |Wn|,  power  dependences  very  well. 

This  plot  was  made  by  scaling  the  laser  power  meter  reading  to  the  value 
of  2a  which  would  give  the  observed  zeroth  step  size  at  this  maximum 
power.  The  same  scaling  factor  was  used  for  the  first  and  second  steps 
as  had  been  done  In  Section  4.3.  Data  from  two  temperatures  Is  shown 
where  the  crosses  denote  the  data  at  4.2  K,  and  the  solid  circles  show 
the  data  at  1.4  K.  Unfortunately  the  first  and  second  step  widths  are 
very  much  smaller  than  expected  from  either  the  Bessel  functions  or 
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FIG.  B.ll  Step  half-widths  as  a  function  of  normalized  laser  Induced 
voltage,  2a,  for  the  Josephson  steps  of  Fig.  B.10,  at  two  temperatures 
T*4.2  K  (crosses),  and  T«1.4  K  (solid  circles).  Solid  line  denotes 
Werthamer  function  fits,  |W|,  and  dashed  lines  denote  Bessel  function 
fits,  |0 { .  The  data  for  the  first  and  second  Josephson  steps  lie 
considerably  below  the  theoretical  curves  but  for  reasons  which  are 
unkown  at  present. 


Werthamer  functions.  Approximate  scaling  factors  for  these  functions 
are  ~0. 1-0. 2  for  the  first  step,  and~0.5  for  the  second  step.  This 
situation  is  much  different  than  seen  in  Section  4.3  with  Sn-Pb  junc¬ 
tions  of  roughly  comparable  capacitance  and  resistance  in  which  the 
Bessel  or  Werthamer  functions  were  not  scaled.  The  need  for  these 
scaling  factors  may  be  related  to  the  rather  large  leakage  current 
observed  which  is  not  expected  from  tunneling  theory.  It  may  be  also 
due  in  part  to  the  fact  that  both  steps  are  in  the  hysteretic  region  of 
the  I-V.  There  is  more  of  a  temperature  dependence  In  the  step  widths 
than  expected  with  nearly  a  factor  of  two  difference  in  the  step  size  at 
the  highest  power  level  for  both  the  first  and  second  step.  We  know  of 
no  simple  reason  for  this  at  present. 

We  conclude  by  noting  that  we  have  been  able  to  couple  604  GHz  laser 
radiation  into  a  SNAP  junction,  and  that  we  have  observed  up  to  three 
Josephson  steps  and  one  photon-assisted  tunneling  step.  The  step  widths 
for  the  first  and  second  steps  are  considerably  smaller  than  expected. 
Further  work  is  needed  to  explain  why  this  is  true.  It  would  also  be 
desirable  to  work  with  junctions  with  areas,  A<1  ym  ,  and  with  higher 
resistance  to  decrease  the  Impedance  mismatch  between  the  junction  and 
the  antenna,  in  order  to  Increase  the  power  coupled  to  the  junction. 
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